προβλημα σε script του matlab

[theo2008]

καλησπερα παιδια.εχο ενα προβλημα σε σκριπτ στο ματλαμπ και θα ηθελα την βοηθεια σας(οχι απαραιτητα την λυση)

Εχω τον παρακατω κωδικα στον οποιο θελω για v=-12:0.05:2 να υπολογισω το ολοκληρωμα teliki για καθε τιμη του v και μετα να υπολογισω το total παλι για καθε τιμη του v και τελος να κανω σε μια plot το v και το total.Το ολοκληρωμα μου θελω να ειναι απο 0 εως +απειρο.Παρακατω παραθετω τον κωδικα καθως και την σκεψη μου για το ολοκληρωμα.

 

A=12.56e-4;

L=0.44;

vdot=v/L;

G=2;

T=300;

q=1.6e-19;

k=8.62e-5;

%Edot μεταξυ 0.01-0.7

Edot=0.01;

%Ndot μεταξύ 10^8-10^15

Ndot=10^8;

%diaf μεταξύ 0.01-0.5

diaf=0.05;

 

paronomasths=1+exp((E-q*vdot)/(k*T));

sinartisi1=1./paronomasths;

sta8eroioroi=(G*Ndot)/(sqrt(pi/2)*diaf);

ek8eths=((E+Edot+q*vdot)/diaf).^2;

sinartisi2=exp(-2*ek8eths);

sinartisi3=sinartisi1.*sta8eroioroi;

teliki=sinartisi3.*sinartisi2;

 

Η σκεψη μου για το ολοκληρωμα ειναι η εξης:αρχικα πριν το paronomasths να βαλω μια syms Ε και μετα την teliki να βαλω το

%ipologismos olokliromatos

subtotal=double(int(teliki,E,0,inf))

total=q*A*L*subtotal

 

Δηλαδη ο κωδικας μου να γινει ως εξης:

 

syms E

paronomasths=1+exp((E-q*vdot)/(k*T));

sinartisi1=1./paronomasths;

sta8eroioroi=(G*Ndot)/(sqrt(pi/2)*diaf);

ek8eths=((E+Edot+q*vdot)/diaf).^2;

sinartisi2=exp(-2*ek8eths);

sinartisi3=sinartisi1.*sta8eroioroi;

teliki=sinartisi3.*sinartisi2;

 

%ipologismos olokliromatos

subtotal=double(int(teliki,E,0,inf))

total=q*A*L*subtotal

 

Αυτο παρατηρησα πως για μια τιμη του v πχ-2 δινει καποιο αποτελεσμα το οποιο το βρισκω λογικο νουμερο.Οταν αλλαζω ομως τιμη στο v πχ-3 τοτε η τιμη του ολοκληρωματος παραμενει η ιδια και συνεπως και του total.Κανω κατι λάθος?Υπαρχει μηπως κανας καλυτερος τροπος υπολογισμου ολοκληρωματος?

Ευχαριστω.

 

 

-----Προστέθηκε 29/11/2008 στις 09 : 43 : 24-----

 

 

Ισως αμα χρησιμοποιησω την αριθμητικη ολοκληρωση του gauss?Αλλα δεν γνωριζω πολλα σχετικα με αυτην την μεθοδο... Καποιος που να μπορει να βοηθησει?