function f=fun(E)
char E;
a=0.8;b=0.7;v0=8.1;e=2.71828;m=1;h=1;
k=(-2*m*E/(h^2))^(1/2);
l=(-2*m*(E+v0)/(h^2))^(1/2);
A=((l-k)*(e^((l+k)*a)-e^((l-k)*a)))/((l+k)*((((l-k)/(l+k))^2)*e^((-l-k)*a)-e^((l+k)*a)));
f=((l+k)*e^((l-k)*(2*a+b)))/((l-k)*e^((-l-k)*(2*a+b)))-((k-l)*e^((k+l)*(a+b))-A*(k+l)*e^((k+l)*(a+b)))/((k-l)*e^((k-l)*(a+b))-A*(k+l)*e^((k-l)*(a+b)));
end
求E在[-v0,0]范围内的数值解
——用solve解出来是
solve(@fun,E)
Warning: Explicit solution could not be found.
> In solve at 81
ans =
[ empty sym ]
——用fsolve没办法得到多解
——想用画图看，但出来的图是{空白的}！！曲线不见了
所以不知道该怎么办了，求大神们支招
Find the treasures in MATLAB Central and discover how the community can help you!
