何度もありがとうございます。
z は理想的な楕円の場合からのずれになりますよね？
z = Ax^2+By^2+Cxy+Dx+Ey+1　
で、A〜Eとx,yが適切ならz=0になると。

f(x,y)=Ax^2+By^2+Cxy+Dx+Ey+1
っておいて
s=倍f(x,y)}^2 になりませんか？