[交流] 【求助】帮忙看个小程序。

Originally posted by lenienter at 2011-03-02 22:23:53:
我写了小程序，用来做数值积分的。老是出错。看不出是怎么回事。请各位高手帮忙看看。
p1=0.04;
p0=0.96;
D=64.03;
n=@(x) 0.00578*512.36*x.^0.5;
N=quadl(n,0.0,128.06);
m=@(x) x.^2.5*0.00578*512.36*( ...

p1 = 0.04;
p0 = 0.96;
D = 64.03;
n = @( x ) 0.00578 * 512.36 * x .^ 0.5;
N = quadl( n, 0.0, 128.06 );
m = @( x ) x .^ 2.5 * 0.00578 * 512.36 .* ( p1 * p0 - 0.5 * p1 *p0 * ( exp( -x ./ 44.58 ) + exp( -x ./ 1.98 ) ) );
M = quadl( m, 0.0, 128.06 );
n2 = p1 - M / N;