Original Post
hi my friends. I have a problem and i want your help.I have the next code.I want to find the integral(0,inf) of final for v=-12:0.05:2.for one value of v,for example -11.95,the code gives me a result.When i try to change the value of v,-10.95 for example the value of subtotal is again the same and i dont know why.Why this happens?is there a problem in my code?Do you know what can i change to have for each value of v the equivalent subtotal and total?Because i want to make a plot(vdot,total) and with this code i think that the result will be a straight line which is wrong. A=12.56e-4; L=0.44; vdot=v/L G=2; T=300; q=1.6e-19; k=8.62e-5; %Edot between 0.01-0.7 Edot=0.01; %Ndot between 10^8-10^15 Ndot=10^8; %DE between 0.01-0.5 DE=0.05; syms E; denom=1+exp((E-q*vdot)/(k*T)); function1=1./denom; static=(G*Ndot)/(sqrt(pi/2)*DE); ek8eths=((E+Edot+q*vdot)/DE).^2; function2=exp(-2*ek8eths); function3=function1.*static; final=function3.*function2; %or if i try final=function1.*function2.*static is correct? subtotal=double(int(final,E,0,inf)) total=q*A*L*subtotal Maybe if i try to solve it with gaussian integral?But i dont know how to use this method... Thank you.