matlab求定积分

求f在区间[a,b]上对t的定积分。
符号解：例如f=t^2*exp(-t);a=0,b=1
>> clear
>> syms t%定义符号变量
>> f=t^2*exp(-t);
>> int(f,t,0,1)

ans =

2 - 5*exp(-1)
数值解：例如例如f=t^2*exp(-t^3);a=0,b=1
>> clear
>> t=0:0.0001:1;%中间的0.0001是步长
>> f=t.^2.*exp(-t.^3);%求节点处的函数值。这一句的"^"和"*"前面的"."表示对应元素乘方、相乘。因为t是向量，只有方阵才可以乘方，矩阵乘法也必须满足相应的维数关系
>> sum(f*0.0001)%节点处的函数值相当于矩形的高，步长相当于矩形的宽，乘积再求和就得到总面积的近似值。

ans =

0.2107

作者：TianSL