Matlab程序(一)多峰函数在曲线上的积分
求取多峰函数在x^2 +(y/3)^2=1所描绘的平面曲线上的积分
思路:将曲线积分转化成一元积分
代码如下:
function=@(t) (3*(1-cos(t)).^2.*exp(-(cos(t).^2) - (3*sin(t)+1).^2) ...
- 10*(cos(t)/5 - cos(t).^3 - 3*sin(t).^5).*exp(-cos(t).^2-3*sin(t).^2) ...
- 1/3*exp(-(cos(t)+1).^2 - 3*sin(t).^2)).*sqrt((sin(t)).^2+(3*cos(t)).^2);
result=integral(function,0, 2*pi);
result即为最终求取的积分值