matlab解方程求助

matlab小白求教

主程序
clc
clear
ys=-1800:5:1800;
u=0.17;
alpha=0;
g=7;
E0=[-0.001,-1000];
Es=zeros(size(ys));
for i=1:numel(ys)
y=ys(i);
fun = @(E)root2d(E,y,u,alpha,g);
[z,fval]=solve(fun,E0);
end
plot(y,E)
子程序
function myfun = root2d(E,y,u,alpha,g)
myfun(1) = 1./(u-g^2./(y-E))-(3pi/8/sqrt(2))(sqrt(-E)-alpha/sqrt(2)atanh(alpha/sqrt(-2E)));
end
求大神帮我看看哪里有问题