I want to solve the equation (iv) as a function of E0 as shown in the attached picture. I tried to calculate it using the following code. I am stuck in this problem! please help
3 次查看(过去 30 天)
显示 更早的评论
if true % code end r=5*10^-9; V=4/3*pi*r^3; Nd=10^17; E0=19; Ld=18.5*10^-9; E=E0+1/6*(r/Ld)^2; Nr=Nd*exp(-(E0+(1/6)*(r/Ld)^2)); eqn=Nd-Nd*exp(-(E0+(1/6)*(r/Ld)^2)); syms E0 eqn=Nd-Nd*exp(-(E0+(1/6)*(r/Ld)^2)); int(eqn,[0,4/3*pi*r^3])
![](/matlabcentral/answers/uploaded_files/137684/IMG_20181023_113111.jpg)
1 个评论
Torsten
2018-10-23
I wonder why your integration is over E0 and not over R:
syms R
eqn = (Nd-Nd*exp(-(E0+(1/6)*(R/Ld)^2)))*4*pi*R^2;
sol = int(eqn,R,0,r)
回答(1 个)
madhan ravi
2018-10-23
r=5*10^-9;
V=4/3*pi*r^3;
Nd=10^17;
E0=19;
Ld=18.5*10^-9;
E=E0+1/6*(r/Ld)^2;
Nr=Nd*exp(-(E0+(1/6)*(r/Ld)^2));
eqn=Nd-Nd*exp(-(E0+(1/6)*(r/Ld)^2));
syms E0
eqn=Nd-Nd*exp(-(E0+(1/6)*(r/Ld)^2));
int(eqn,[0,4/3*pi*r^3]) %no errors occurred
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!