% sind(x) * cosd(y)

x=0:30:180 
[x1 x2]=size(x);

y=0:30:360
[y1 y2]=size(y);

xx = x(ones(y2,1),:)

yy = y(ones(x2,1),:)'

xxv = sind(xx)
yyv = cosd(yy)
xy = xxv .* yyv


