怎么用MATLAB画阿基米德螺线

怎么用MATLAB画阿基米德螺线
给出程式

阿基米德螺线 是所有形式为（极坐标方程）r = aθ

假设（你可以自己改成其他的）在程序中a=3，用b表示θ[0,3*pi]

一个命令    ezpolar('2*b',[0,3*pi])

图形如 

曲线的长度计算∫sqrt（[x'[t]^2 + y'[t]^2 + z'[t]^2] ）dt 

极坐标是r = aθ，直角坐标式x=aθcos（θ），y=aθsin（θ)

dl=sqrt（[x'[θ]^2 + y'[θ]^2 ）dθ=|a|*sqrt(θ^2+1)dθ

长度L=|a|*∫sqrt(θ^2+1)dθ

matlab程序 

sym x;

L=int('3*sqrt(x^2+1)',0,3*pi)%符号解

L1=vpa(L)%符号转成数值解

结果L =9/2*pi*(9*pi^2+1)^(1/2)-3/2*log(-3*pi+(9*pi^2+1)^(1/2))

 L1 =138.39649839205735200693400612079

可以自己改a和θ的范围