format longD=3400;f=4267;L=6350;a=@(x)fsolve(@(a) D/(2*f)-sin(a+x),zeros(1,length(x)),optimset('display','off'));y=@(x) L*sin(a(x))-D*0.5+0.5*D*cos(a(x))*tan(x)/sin(a(x)+x)+0.5*D*sin(a(x))/sin(a(x)+x);fplot(y,x);ylabel('h')grid on
里面还有一条语句函数x=linspace(0,60);