Question

matlab求数值积分

大家帮我看一下这个求积分的式子怎么错了？？

Answer 1

使用quad有问题
>> quad('0.0732*exp(-(x+0.02287).^2/64.90)',0,inf)
Warning: Infinite or Not-a-Number function value encountered.
> In quad at 109

ans =

NaN
使用int就可以了
>> syms x
>> int(0.0732*exp(-(x+0.02287).^2/64.90),0,inf)

ans =

-(183*10^(1/2)*649^(1/2)*pi^(1/2)*(erf((2287*10^(1/2)*649^(1/2))/64900000) - 1))/50000

>> double(ans)

ans =

0.520936708170281