清空筛选条件
清空筛选条件

Numerical Integration problem

1 次查看（过去 30 天）
Keith
Keith 2011-11-30
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

请先登录，再回答此问题。

回答（1 个）

Andrei Bobrov
Andrei Bobrov 2011-11-30
use quadgk
for quadgk(F,d,z), if z = 4 then F(4) = NaN, division by zero.

类别

在 Help Center 和 File Exchange 中查找有关 Develop Apps Using App Designer 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by