求助——积分求面积遇到问题

丫头丫头丫头

如图，x2,y2为转子方程，我想求x2,y2曲线围成的面积，用matlab编写程序如下：
z1=6;
z2=z1-1;
R=48.78;
r=8.13;
Q=R-r;
e=7.05;
syms theta
k=12.6;
f=r/e;
a=sin(z1*theta);
b=f+cos(z1*theta);
m=atan(a./b);
p=m-theta;       %theta
c=sin(theta+p);
n=asin(f*c);
o=n-p;          %phi
g=(z1/z2)*o;    %psi

x2t=Q*sin(theta)+e*sin(z2*theta)-k*sin(p);       %原始齿形
y2t=Q*cos(theta)-e*cos(z2*theta)+k*cos(p);
x3t=x2t.*cos(o-g)-y2t.*sin(o-g)-e*sin(g);        %共轭齿形
y3t=x2t.*sin(o-g)+y2t.*cos(o-g)-e*cos(g);


d2=diff(x3t);
ds2=inline(y3t.*d2);
s2=5*quadl(ds2,0,2*pi/5)
但是运行出来有问题，请高手指点，转子的面积该如何求，先谢谢啦{:soso_e163:}