怎样用matlab求定积分

来源:百度知道 编辑:UC知道 时间:2024/09/22 15:35:41
738.5*x^2*exp(-0.00105*x^2) 积分范围为0到1.9

%by dynamic
%see also http://www.matlabsky.com
%contact me matlabsky@gmail.com
%2009.2.
%
积分Matlab中提供了数值和解析解法两种

1.数值方法
matlab中所有数值积分函数的汇总比较和使用说明参见我写的教程http://www.matlabsky.com/thread-359-1-3.html
>> fun=@(x)738.5*x.^2.*exp(-0.00105*x.^2);
>> quadl(fun,0,1.9)

ans =

1.6846e+003

2.符号计算
syms x
>> syms x
>> int(738.5*x^2*exp(-0.00105*x^2),x,0,1.9)

ans =

-2004500/3*exp(-7581/2000000)+52750000/63*pi^(1/2)*42^(1/2)*erf(19/2000*42^(1/2))

>> vpa(ans)

ans =

1684.62230190651358892102680947

int('738.5*x^2*exp(-0.00105*x^2)',0,1.9)
结果为:1684.6223 。

int('738.5*x^2*exp(-0.00105*x^2)',x,0,19)