Анализ функции

Доброго времени суток.

Имеется R-функция вида:
a^2-x^2+b^2-y^2-sqrt((a^2-x^2)^2+(b^2-y^2)^2)=0

a=20, b=40.

В матлабе пишу программу вида:

[x, y] = meshgrid(-400:1:400, -400:1:400);
a=20; b=40
Z=a.^2-x.^2+b.^2-y.^2-sqrt((a.^2-x.^2).^2+(b.^2-y.^2).^2);
[C,h]=contour(x,y,Z); grid
clabel(C,h)
title ('contour'); xlabel('x'); ylabel('y')

Получаю следующее построение.
http://xmages.net/i/3330388

Подскажите, пожалуйста, В чем причина того, что получившиеся линии уровня кроме первой имеют не прямоугольный вид. В чем причина их скругления? И как будет выглядеть фигура, если например "a" выразить через "y".

Буду признателен за помощь.



Редактировалось 1 раз(а). Последний 25.01.2012 14:43.

Спасибо за анализ.