function do1()
fid=fopen('fig1-2.txt','w');
syms x y f n o w 
for x=0.5:0.1:10
    o=.1053722210e-1*x*((5.233778845-.1053722210e-1*i*x)/log(1.+337.5000000/(-1.944694560+2.449489743*y^.5)^2))^(1/2);
    w=(7161972.441+7161972.441*i)*besselj(1,(.2513274122e-1+.2513274122e-1*i)*x^.5)/(x^1.5*besselj(0,(.2513274122e-1+.2513274122e-1*i)*x^.5)); 
    eqn=.1666666667e-3*w/y/(1+.1580246913e-4*w*log(1+1350.000000/(1.341640787+y^.5)^2)*cos(o))+(1.8-y)/(1.8+y);
    f=solve(eqn,y);
    fprintf(fid,'%f',x);
    fprintf(fid,'%s','      ');
    fprintf(fid,'%f',double(real(f)));
    fprintf(fid,'%s','      ');
    fprintf(fid,'%f',double(imag(f)));
    fprintf(fid,'%s\n','');  
end
sta=fclose(fid);

谁能帮我看看这个符号函数方程有什么错误吗，总是运行出问题，请好心人帮帮忙！