XProfan

Título de la ventana "Trashware: Planetenbahnsimulation para 'rumspielen"
'(CL)CopyLeft 2014-01 by P.Pájaro carpintero, Wien
Ventana 0,0-%maxx,%maxy-40:randomize:set("decimals",5):set("numwidth",9)
Font 2:var xx&=width(%hwnd)\2:var yy&=height(%hwnd)\2
var grad!=0.0174532925199432958:var rad!=57.2957795130823209
'=========================================================================
'    Experimentierfeld
'-------------------------------------------------------------------------
var G!=10^-17 : var m1!=10^20 : var v2!= 1'Periheldrehungen des Merkur
'var G!=10^-17 : var m1!=10^20 :var v2!= 2.2 ' Kreisbahn
'var G!=10^-15 : var m1!=10^20 :var v2!= 2.2 ' Más Gravitation?
'var G!=5*10^-18 : var m1!=10^20 :var v2!= 2 ' Pluto
'var G!=10^-15 : var m1!=10^30 :var v2!= 2.2 '
'=========================================================================
'var x2!=100:var y2!=0 ' Einschuss-Bahnpunkt
'     var w2!=90*grad! ' Einschusswinkel
var x2!=100:var y2!=0:var w2!=90*grad!
'=========================================================================
declarar a2!,rr!,alpha!:declarar a2x!,a2y!,v2x!,v2y!:declarar n&
cls $cdcdcd:line 20,yy& - 2*xx&-20,yy&:line xx&,20 - xx&,2*yy&-20

repeat

    locate 1,1:'imprimir n&,::inc n&
    imprimir x2!,y2!,v2!,w2!*rad!;:usepen 1,22,rgb(255,255,0):line xx&,(yy&) - xx&+1,yy&
    usepen 0,12,rnd(8^8):line xx&+x2!,(yy&-y2!) - xx&+1+x2!,yy&-y2!
    alpha!=ArcTan4(x2!,y2!)
    rr!=sqr(x2!)+sqr(y2!):caso rr!<10000:rr!=10000
    a2!= -G!*m1!/rr!
    v2x!=v2!*cos(w2!)+a2!*cos(alpha!)/2
    v2y!=v2!*sin(w2!)+a2!*sin(alpha!)/2
    w2!=ArcTan4(v2x!,v2y!)
    x2!=x2!+v2x!
    y2!=y2!+v2y!
    v2!=sqrt(sqr(v2x!)+sqr(v2y!))'*1.001
    waitinput 1:until %key=27:end

    proc ArcTan4 :parámetros x!,y!

        var pi!=3.1415926535897932:var w!=0'(CL)2014-01 P.Pájaro carpintero§gmx.at

        if x!=0:if y!>0:w!=pi!*0.5:elseif y!<0:w!=pi!*1.5:más :w!=0:endif :volver w!:elseif x!>0

            if y!=0:w!=0:volver w!:elseif y!>0:if x!>y!:w!=arctan(y!/x!):más :w!=pi!/2-arctan(x!/y!):endif

                volver w!:más :if x!<-y!:w!=pi!*1.5+arctan(x!/-y!):más :w!=2*pi!-arctan(-y!/x!):endif :volver w!
                endif :más :if y!>0:if x!>-y!:w!=pi!/2+arctan(-x!/y!):volver w!:más :w!=pi!-arctan(y!/-x!)
                volver w!:endif :elseif y!<0:if x!<y!:w!=pi!+arctan(-y!/-x!):más :w!=pi!*1.5-arctan(-x!/-y!)
                endif :volver w!:más :w!=pi!:volver w!:endif :endif :Imprimir " ArcTan4 ERROR":waitinput:waitinput

            ENDPROC