how to solve the given error

Rajeev Ray

2016 4 4

0 个回答

4 次查看（30 天）

0 个投票

me=9.109*10^-31;
h=6.626070040*10^-34;
hr=h/(2*pi);
c1=((2*0.28*me).^1.5)/(2*pi*(hr).^3);
r=2.3*10^-9;
V0=1.33*pi*(r^3);
y1=V0*c1;
E=315;
delta=50;
deltaG1=[-0.6 -0.61 -0.71 -0.79 -0.84];
H=0.011;
rho=@(E_s)y1*E.^.5;
fun=@(E,E_s)rho(E_s)/(delta*sqrt(2*pi)).*exp(-(E-E_s).^2/(2*delta^2));
values=integral(@(E_s)fun(E,E_s),0,Inf);
values
fun3=@(hr,values,E1)(2*pi/hr)*values*(0.011.^2).*(1/((4*pi*400*25.7).^0.5)).*exp(((400+(-0.61)+E1).^2)/((4*400*25.7).^0.5));
values_s=integral(@(E1)fun3(hr,values,E1),-Inf,Inf);
values =
     4.4609e+31
Warning: Infinite or Not-a-Number value encountered. 
> In integralCalc/iterateScalarValued (line 349)
  In integralCalc/vadapt (line 132)
  In integralCalc (line 103)
  In integral (line 88)