[laosam280]

应该是这样的吧？
b=2*atan((y.^2-x^2*1.483^2)./(x^2*1.49496^2-y.^2))./sqrt(x^2*1.49496^2-y.^2);
x=2*pi/1.550918;
y=0:0.002:8;
那你干脆用循环算了，将b的结果保存在一个矩阵里。
b=ones(length(y),1);
for i=1:400
y=0:0.02:8
x=2*pi/1.550918;
b(i)=2*atan((y(i)^2-x^2*1.483^2)./(x^2*1.49496^2-y(i)^2))./sqrt(x^2*1.49496^2-y(i)^2);
end
plot(b,y,'.-')

呵呵，大概就是这样吧