Optimización de código en Matlab

Implementar una función que calcule los cuadrados de los 1000 y 10000 primeros números naturales de 3 formas distintas:
• Utilizando un bucle for sin declarar las variables.
• Utilizando un bucle for declarando previamente la variable.
• Calculando el cuadrado del vector correspondiente (operación componente a
componente).
Calcular, con la ayuda de la instrucción de Matlab cputime el tiempo que se tarda en realizar el cálculo con cada uno de los tres métodos. La función a implementar debe tener como variable de entrada el número N, y como salida el tiempo que se tarda con cada uno de los tres métodos (t1, t2, t3).


fprintf('Para %d\n',N)
 
tic;
for i=1:1:N
    cuadrados(i)=i^2;
end
toc;


tic;
cuadrados=1:1:N;
for i=1:1:N
    cuadrados(i)=i^2;
end
toc;

tic;
numeros=1:1:N;
cuadrados=numeros.^2;
toc;

end