Is it possible make a for double loop simpler.
x=linspace(-2,2,100); s=1/2;[X,Y]=meshgrid(x,x);A = X.^2; B = Y;for k=1:size(x,1) for j=1:size(y,2) if s<X(k,j) L(k,j)=A(k,j).*sin(s); M(k,j)=B(k,j).*sin(s); else %(s>X(k,j))
L(k,j)=A(k,j).*cos(s/2); M(k,j)=B(k,j).*cos(s/2); end endend
Best Answer