求助：matlab编程实现椭圆齿轮,请大家帮忙，谢谢了。

参考：基本参数确定后，请问先要确定每个齿的位置吗？
弧长的变上限积分应该怎么用程序实现
   

椭圆方程x=a*cos（x），y=sin（x）；弧长=a*（1-e^2*sin(x)^2）在（0，x）上的积分，
齿形方程：XL=x-s*cos（au）*cos（beita+au）
          YL=y-s*cos（au）*sin（beita+au）   (左齿形)
               XR=x+s*cos（au）*cos（beita-au）
          YR=y+s*cos（au）*sin（beita-au）   (右齿形)
其中：tan（beita）= - a/b*(ctan(x)).

还可以参考这个pdf文件

最近也在做这些呢

期待有人回答。