求各位大神帮忙看看这个程序错在哪儿
clear;
n=1.526;
k=6.851;
for x=0:90
for y=0:180
if n==sin(x)/sin(y)
b=exp(-k*10/cos(y));
r=0.5*(sin(x-y)^2/sin(x+y)^2+tan(x-y)^2/tan(x+y)^2);
A=(1-r)^2*b/(1-r^2*b^2);
end
end
end
fplot(x,A,'r')
set(gcf,'Name','透光率与入射角关系','Color','r')
xlabel('入射角'),ylabel('透光率')
