plot {x = 16 pow(sin(phi), 3), y = 13 cos(phi)-5 cos(2 phi)-2 cos(3 phi)-cos(4 phi)}