%
% Plots a nodal cubic
%


t = [-1.3: 0.01: 1.3];

x = t.^2 - 1;
y = t.* (t.^2 - 1);

plot(x,y)


%
% Lemniscate of Bernouilli
%
%
% half
%

t = [-1: 0.01: 1];

x = t.*(1 + t.^2)./(1 + t.^4);
y = t.*(1 - t.^2)./(1 + t.^4);

figure(1)
plot(x,y)

%
% second half
%

x1= t.*(1 + t.^2)./(1 + t.^4);
y1 = t.*(t.^2 - 1)./(1 + t.^4);


figure(2)
plot(x, y, x1, y1)
