Function getLineInfo#[](point1#[], point2#[])     'Rückgabewert: [0] = Steigung, [1] = Y-Achsenabschnitt
                                                  '=> pointN[1] = result[0]*pointN[0]+result[1]
  Local result#[] = [0.0, 0.0]
  result[0] = (point1[1]-point2[1])/(point1[0]-point2[0])
  result[1] = point1[1]-point1[0]*result[0]
  Return result
End Function
Function pointInPoly%(poly#[][], point#[])    'NUR konvexe, geordnete Polygone, da lineare Ungleichungen
  SetColor 0, 255, 0
  Local line#[], cmp1%, cmp2%, third%
  For i=1 To poly.length
    cmp1  =  i    Mod poly.length  'cmp1->cmp2 wird mit point verglichen
    cmp2  = (i+1) Mod poly.length
    third = (i+2) Mod poly.length   'Der 3. Punkt gibt die Seite an, auf der point liegen soll
    If poly[cmp1][0] = poly[cmp2][0] Then    'X-Werte gleich => kein getLineInfo möglich
      If Sgn(point[0]-poly[cmp1][0]) <> Sgn(poly[third][0]-poly[cmp1][0]) Then
            '3. Punkt und Vergleichspunkt auf unterschiedlichen Seiten
        Return False
      EndIf
    Else
      line = getLineInfo(poly[cmp1], poly[cmp2])
      '=> line[0]*triangle[cmp1/2][0]+line[1] = triangle[cmp1/2][1]
      If Sgn(point[1]-(line[0]*point[0]+line[1])) <> Sgn(poly[third][1]-(line[0]*poly[third][0]+line[1])) Then
          '3. Punkt und vergleichspunkt auf unterschiedlichen Seiten
        Return False
      EndIf
    EndIf
  Next
  Return True
End Function



'''TEST'''
Local points#[][] = [[5.0, 0.0], [100.0, 100.0], [100.0, 200.0], [5.0, 300.0]]

Local point#[] = [1.0, 2.0]
Graphics 800, 600
Local poly#[] = [1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0]
For i=0 To 3
  poly[i*2]   = points[i][0]
  poly[i*2+1] = points[i][1]
Next
Local timer:TTimer = CreateTimer(30)
Repeat
  WaitTimer timer
  point[0] = MouseX()
  point[1] = MouseY()
  SetColor 255, 255, 255
  DrawPoly poly
  If pointInPoly(points, point) Then
    SetColor 0, 0, 255
  Else
    SetColor 255, 0, 0
  EndIf
  DrawOval point[0]-2, point[1]-2, 4, 4
  Flip
  Cls
Until AppTerminate() Or KeyHit(KEY_ESCAPE)
End