怎么求出数值解并画出图形

zdd13717421

程序为：
syms w1 m
xj=9; q=1.6*10^(-19); delta=10^(-4);
d1=1.5; d2=0.19; d3=0.4; d4=0.4;
N1=2.5*10^15; N2=2*10^16; N3=6*10^16; N4=2*10^15;
v=0.12; x=9:0.01:11.45;
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj
g1=int(e01*xj/m-q*N4*(m^2-xj^2)/(m*2*delta),m,xj,xj+w1)+v
w1=solve(g1,'w1');
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj;
y=(e01*xj./x-q*N4.*(x.^2-xj^2)./(x.*2*delta)).*(xj<=x&x<=xj+w1);
plot(x,y)
错误提示：

Warning: Explicit solution could not be found.
> In solve at 81
  In f1 at 17
??? Error using ==> mupadmex
Error in MuPAD command: Either base or exponent must be a scalar.
Error in ==> sym.sym>sym.mpower at 198
            B = mupadmex('mllib::mpower',A.s,p.s);
Error in ==> f1 at 18
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj;

[ 本帖最后由 ChaChing 于 2010-7-13 23:51 编辑 ]

ChaChing

个人水平专业有限, LZ没交代下相关资讯, 真的有点难懂!
但好奇, 既然要求数值解, 为何还使用sym?