перевести програмку из дельфи в си

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, ExtCtrls;
 
type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Button1: TButton;
    PaintBox1: TPaintBox;
    Button2: TButton;
    Button3: TButton;
    Button4: TButton;
    Edit2: TEdit;
    procedure Button1Click(Sender: TObject);
    procedure FormCreate(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 {$R *.dfm}
 
function max(a,b: extended):extended;
begin
  if a>b then max:=a
         else max:=b;
end;
 
function min(a,b: extended):extended;
begin
  if a<b then min:=a
         else min:=b;
end;
 
procedure TForm1.Button1Click(Sender: TObject); {????? ???-???????}
  const N=300;
  type Tmas = array [0..N,0..trunc(N/2)] of double;
  var i,m,code,K:integer;
      u, v : Tmas;
      c, tao, h: extended;
begin
  Val(Edit1.Text,c,code);
  tao:=1/(3*c*N);
  h:=1/N;
  K:= trunc(N/2);
  for i:=0 to K-50 do
    u[i,0]:=1;
  for i:=K-49 to K+49 do
    u[i,0]:=0.5;
  for i:=K+50 to N do
    u[i,0]:=0;
  for m:=0 to K do
    u[0,m]:=1;
  for m:=0 to K do
    u[N,m]:=0;
  for m:=0 to K-1 do
    begin
    v[0,m+1]:=u[0,m] - (c*tao/h)*(u[1,m] - u[0,m]);
    for i:=1 to N-1 do
      begin
      v[i,m+1]:=u[i,m] - (c*tao/h)*(u[i+1,m] - u[i,m]);
      u[i,m+1]:=1/2*(u[i,m] + v[i,m+1] - (c*tao/h)*(v[i,m+1]-v[i-1,m+1]));
      end;
    end;
  with PaintBox1.Canvas do
    begin
      MoveTo(0,200);
      for i:=0 to N do
       begin
             Pen.Width:=1;
             Pen.Style:=psSolid;
             Pen.Color:=clBlue;
             LineTo(i*trunc(600/N),4*trunc(150 - 100*u[i,K]));
             Pen.Color:=clRed;
            if (u[i,K]<0.99) and (u[i,K]>0.03)
            then
            begin
            Pen.Width:=3;
            MoveTo(i*trunc(600/N)-1,4*trunc(150 -100*u[i,K]+1));
            LineTo(i*trunc(600/N)+1,4*trunc(150 -100*u[i,K]-1));
            MoveTo(i*trunc(600/N)+1,4*trunc(150 -100*u[i,K]+1));
            LineTo(i*trunc(600/N)-1,4*trunc(150 -100*u[i,K]-1));
            end;
       end
    end;
end;