椭圆参数方程：x=acost,y=bsint
则曲率z=(d^2y/dx^2)/[1+(dy/dx)^2]^(3/2)
将dy/dx,及d^2y/dx^2算出来代进去即可