data Butterfly;
do t=-400 to 4 by 0.01;
a = cos(t);
b = cos(4*t);
x = sin(t)*(exp(a)-2*b-(sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)));
y = cos(t)*(exp(a)-2*b-(sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)));
output;
end;
do k=-2 to 2 by 0.001;
x1=k;
y1=abs(k)+2.5;
output;
end;
proc plot ;
plot y*x='*' y1*x1='*'/overlay;
run;