689:
690:
691:
692:
693:
694:
695:
696:
697:
698:
699:
700:
701:
702:
703:
704:
705:
706:
707:
708:
709:
710:
711:
712:
713:
714:
715:
716:
717:
718:
719:
720:
721:
722:
723:

{ ... }
function ArcCos(const X: Extended): Extended;
begin
  Result := ArcTan2(Sqrt(1 - X * X), X);
end;

function ArcSin(const X: Extended): Extended;
begin
  Result := ArcTan2(X, Sqrt(1 - X * X))
end;

function ArcTan2(const Y, X: Extended): Extended;
asm
        FLD     Y
        FLD     X
        FPATAN
        FWAIT
end;

function Tan(const X: Extended): Extended;
{  Tan := Sin(X) / Cos(X) }
asm
        FLD    X
        FPTAN
        FSTP   ST(0)      { FPTAN pushes 1.0 after result }
        FWAIT
end;

function CoTan(const X: Extended): Extended;
{ CoTan := Cos(X) / Sin(X) = 1 / Tan(X) }
asm
        FLD   X
        FPTAN
        FDIVRP
        FWAIT
end;					

1:

angle:= arccos((shadow.ll.Position.Y-shadow.ul.Position.Y)/sqrt((shadow.ul.Position.X-shadow.ll.Position.X)*(shadow.ul.Position.X-shadow.ll.Position.X)+(shadow.ul.Position.Y-shadow.ll.Position.Y)*(shadow.ul.Position.Y-shadow.ll.Position.Y));					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:

function fac(n: Byte): Integer;
var
  i, j: Integer;
begin
  j := 1;
  for i := 1 to n do
    j := j * i;
  Result := j;
end;

{ Arkuscosinus }
function _ArcCos(r: Extended): Extended;
var
  i: Byte;
  s: Extended;
begin
  s := 0;
  for i := 0 to 15 do
  { 15 ggf. durch eine andere Zahl zu ersetzen;
    je nachdem, wie gut / schlecht die Näherung
    sein muss / darf bzw. die Zeit. }
  begin
    s := s + fac(2 * i) / (Math.Power(2, 2 * i) * Math.Power(fac(i), 2)) *
         Math.Power(r, (2 * i + 1)) / (2 * i + 1);
  end;
  Result := DegToRad(pi / 2 - s);
end;