function heart(a,b,c)
% good param choices:  heart(1/3,3/2,1/8)  =  regular heart
%                      heart(1/3,3/2,1/3)  =  narrower pt
for t=1:360
% note this is just a cardioid with an angular scaling function
% that has a peak near 0 degrees, with =1 everywhere else
r(t)=(1+cos((t-180)*pi/180))*(1+a*b^(-abs(c*(t-180))));
end;
t=1:360;
t=t-90;  % rotate heart so it's vertical
t=(t-180)*pi/180;
polar(t,r,'r');
title('r  =  ( 1 + cos(t-90) )( 1 + 0.3*1.5^-^|^0^.^3^*^(^t^-^9^0^)^| )');