Graphics 800, 600

Local X1:Float = 200.0
Local Y1:Float = 200.0
Local X3:Float = 500.0
Local Y3:Float = 400.0
Local H:Float = 50.0

Local OldY:Int = Y1
For Local X:Int = X1 + 5 To X3 Step 5
	Local Y:Int = GetParabolaY(X, X1, Y1, X3, Y3, H)
	
	DrawLine X - 5, OldY, X, Y
	
	OldY = Y
Next

SetColor 255, 0, 0
DrawLine 0, Y1 - H, 800, Y1 - H 'Maximalhöhe einzeichnen

SetColor 255, 255, 0
DrawOval X1 - 5, Y1 - 5, 10, 10 'Startpunkt
DrawOval X3 - 5, Y3 - 5, 10, 10 'Endpunkt

Flip 0
WaitKey
End

Function GetParabolaY:Float(ParabolaX:Float, X1:Float, Y1:Float, X3:Float, Y3:Float, H:Float)
	Local Y2:Float = Y1 - H 'Die Y-Achse ist beim Zeichnen ja umgekehrt, darum -H
	
	'Das hier kann man natürlich beschleunigen, indem man X2 und A vorausberechnet und die ^2 durch Multiplikationen ersetzt, aber ich lasse das mal als "Exercise for the reader"
	Local X2:Float = (-Sqr((Y2 - Y1)*(Y2 - Y3)*(X1 - X3)^2) + X1*(-Y3) + Y2*X1 + X3*Y1 - Y2*X3)/(Y1 - Y3)
	Local A:Float = -(Y2 - Y3)*(Y3 - Y1)^2/(Sqr((Y2 - Y1)*(Y2 - Y3)*(X1 - X3)^2) + X1*Y3 - Y2*X1 - X3*Y3 + Y2*X3)^2
	
	Return A*(ParabolaX - X2)^2 + Y2
End Function