clear all; close all; clc; %define curve phi = linspace(0,2*pi,512); h = @(phi)2-2*sin(phi)+(sin(phi).*sqrt(abs(cos(phi))))./(sin(phi)+1.4); %transform to cartesian coordinates x = h(phi).*cos(phi); y = h(phi).*sin(phi); %plot result plot(x,y) axis equal; title('A lovely gift.')