problem with integral - Programmers Heaven

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Categories

Welcome to the new platform of Programmer's Heaven! We apologize for the inconvenience caused, if you visited us from a broken link of the previous version. The main reason to move to a new platform is to provide more effective and collaborative experience to you all. Please feel free to experience the new platform and use its exciting features. Contact us for any issue that you need to get clarified. We are more than happy to help you.

problem with integral

theo2008theo2008 Posts: 2Member
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...
Sign In or Register to comment.