各位高手帮忙看看这个程序的问题在哪？ 解b1和c1的方程，出来的应该是一条曲线




 t=0:0.01:1;
h=0.5+0.05.*t.^2;
m=0.65;
n=0;
y=m-h-(1/4);
z=n-(1/5);
b1=0:pi/100:pi;
f1='0.03+0.25*cos(b1-1.0472)+0.34*sin(b1+0.5236-c1)=sqrt(0.04+y^2)';
f2='0.052-0.25*sin(b1-1.0472)+0.34*cos(b1+0.5236-c1)=-z';
plot(t,b1)