Как изобразить кривые Ламе в octave или matlab?

Question

[ZaurK]

Здравствуйте! Пытаюсь изобразить кривые Ламе в octave как тут, но получается лишь одна ветвь. Подскажите, что делаю не так?

% Кривые Ламе
n = 0.5
a = 1.0
b = 1.0
f = 0:pi/180:pi/2
x = a*cos(f).^(2/n)
y = b*sin(f).^(2/n)
figure()
plot(x, y)
grid on
axis equal
xlab = xlabel('Ось абсцисс Х');
ylab = ylabel('Ось ординат Y');
title('Кривые Ламе');

Answer 1

[U235U235]

f = 0:pi/180:2*pi

Answer 2

[D_dMer]

%Суперэллипс
%clear all;
m=1;
n=1/2;
a = 10;
b = 10;
f = 0:pi/360:2*pi;
%f = 0:pi/180:pi/2;
x = (cos(f).^(2/n)).*(a*sign(cos(f)));
y = (sin(f).^(2/n)).*(b*sign(sin(f)));
%figure()
plot(x1, y1, x, y)
grid on
axis equal
xlab = xlabel('Ось абсцисс Х');
ylab = ylabel('Ось ординат Y');
title('Кривые Ламе');