Интересная функция

unit Unit1;
 
{$mode objfpc}{$H+}
 
interface
 
uses
  Classes, SysUtils, Forms, Controls, Graphics, Dialogs, StdCtrls, ComCtrls,
  ExtCtrls, TAGraph, TASeries, Math;
 
type
 
  { TForm1 }
 
  TForm1 = class(TForm)
    Chart1: TChart;
    Chart1LineSeries1: TLineSeries;
    Image1: TImage;
    Label1: TLabel;
    TrackBar1: TTrackBar;
    procedure FormShow(Sender: TObject);
    procedure TrackBar1Change(Sender: TObject);
  private
  ak:integer;
  procedure DoGraph;
  public
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.lfm}
 
{ TForm1 }
 
//-------Решение куб. уравнения по формуле Кардано
function Kardano(a,b,c,d:Extended):Extended;
var p,q,qq,alfa,beta:Extended;
begin
p:=(3*a*c-b*b)/(3*a*a);
q:=(2*b*b*b-9*a*b*c+27*a*a*d)/(27*a*a*a);
qq:=p*p*p/27+q*q/4;
if qq>0 then
   begin
   alfa:=power(-q/2+Sqrt(qq),1/3);
   beta:=power(-q/2-Sqrt(qq),1/3);
   Result:=alfa+beta;
   end;
if qq=0 then Result:=0;
if qq<0 then Result:=0;
end;
 
//------Функция в полярных коррдинатах сводится к решению куб. уравнения
function MyFunc(ap,fi:Extended):Extended;
var cf,sf,x:Extended;
a,b,c,d:Extended;
begin
a:=1;
cf:=cos(fi*PI/180);sf:=sin(fi*PI/180);
b:=ap*cf*sf*sf;
c:=0;
d:=-850*sf;
Result:=Kardano(a,b,c,d);
end;
 
//-----------График функции (x^2+y^2)^2+-a*x*y^2=850y    -43<=a<=43
//-----------Полярн.коорд r^3+-a*cos(ф)*sin^2(ф)*r^2-850*sin(ф)=0;     r=f(ф)
procedure TForm1.DoGraph;
var i:integer;
r,x,y:Extended;
begin
Chart1LineSeries1.Clear;
if ak<0 then Chart1LineSeries1.LinePen.Color:=clRed else Chart1LineSeries1.LinePen.Color:=clBlue;
for i:=0 to 180 do
   begin
   r:=MyFunc(Form1.ak,i);
   x:=r*cos(i*PI/180);
   y:=r*sin(i*PI/180);
   if i=0 then begin x:=0;y:=0;r:=0;end;
   Chart1LineSeries1.AddXY(x,y);
   end;
end;
 
procedure TForm1.TrackBar1Change(Sender: TObject);
begin
ak:=TrackBar1.Position;
Label1.Caption:='a='+IntToStr(ak);
DoGraph;
end;
 
procedure TForm1.FormShow(Sender: TObject);
begin
ak:=0;DoGraph;
end;
 
end.