又一个牛B函数：阴阳图极坐标函数

刚才在网上发现了上面这张猛图。急着想验证一下，但不知道Mathematica如何画极坐标的隐函数，于是写了一个Free Pascal的小程序。大家也可以试着把这个小程序粘贴到Free Pascal里运行一下看看。

{$ASSERTIONS+}
uses graph,math;
 
function f(t,r:extended):extended;
begin
   exit((cos(t-r)-sin(t)) * (r**4-2*r**2*cos(2*t+2.4)+0.9)+(0.62*r)**1000);
end;
 
var
   i,j,gd,gm:integer;
   t,r:extended;
begin
   gd:=D8bit;
   gm:=m640x480;
   InitGraph(gd,gm,'');
   Assert(graphResult=grOk);
 
   for i:=0 to 1000 do
   for j:=0 to 1000 do
   begin
      t:=(2*pi)*i/1000;
      r:=2*j/1000;
      if f(t,r)<0 then
PutPixel(round(cos(t)*r*100)+300,round(sin(t)*r*100)+200,15)
end;
readln;
CloseGraph;
end.

果然不出所料，下面就是程序运行的结果：

每一次提到这类的东西，我都有一种把原来的那些牛B东西一起列出来的冲动。大家不要说我赚PV哦。