用MATLAB用定积分求椭圆周长

利用MATLAB用 定积分计算椭圆的周长
怎么利用MATLAB然后用定积分计算椭圆周长
实在是不会了

画个椭圆：

img


椭圆的方程为

x = asinφ，y = bcosφ其微分为

dx = acosφdφ，y = -bsinφdφ弧元为

img

img

所以可以编写.m文件:

%x轴长
a=4;
%y轴长
b=3;
syms theta;
f=4*a*sqrt(1-abs(a.^2-b.^2)/a.^2*sin(theta).^2);
C=int(f,0,0.5*pi);
vpa(C)

在命令行会输出值

img