Numerical Integration problem
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Hi can anyone advise: I am trying to integrate numerically an equation from z_star_=25 down to d=4. However, the solution derives NaN at 6 down to 4.... The code/equation:
d=4;
z_star=25;
u_star_z_star=0.5;
k=0.4;
F = @(x)(((sin((pi()/2)*((x-d)/(z_star-d))).^1.28).^0.333)./(x-d))*(u_star_z_star/k);
%%%--- TABULATING SOLUTION BETWEEN LIMITS --- %%%
A=zeros(85,2);
a=1;
b=1;
z=z_star
for z=z:-0.25:d;
Q(a) = quadv(F,d,z);
D(b)=z
a=a+1;
b=b+1
end
for i=1:length(A);
A(i,1)=D(i);
A(i,2)=Q(i);
end
0 Kommentare
Antworten (1)
Andrei Bobrov
am 30 Nov. 2011
for quadgk(F,d,z), if z = 4 then F(4) = NaN, division by zero.
0 Kommentare
Siehe auch
Kategorien
Mehr zu Numerical Integration and Differentiation finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!