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.