Orbital Physics

-=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- (c) WidthPadding Industries 1987 0\|602\|0 -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=- -=+=-
SoCoder -> Snippet Home -> Collisions and Physics

shroom_monk	Created : 27 April 2007 Edited : 27 April 2007 System : Windows Language : Blitz Orbital Physics Continuous Freefall
Basically, it's a little red box using continuous freefall to orbit around a green circle. Like a planet and a sun. AppTitle "Orbit" Graphics 400,300,24,2 SetBuffer BackBuffer() Local sx# = 100,sy# = 150 Local sxv# = 0,syv# = Pi timer = CreateTimer(50) While Not KeyHit(1) WaitTimer(timer) sxv# = sxv# + (Cos#(ATan2#(sy# - 150,sx# - 200) + 180) * 0.1) syv# = syv# + (Sin#(ATan2#(sy# - 150,sx# - 200) + 180) * 0.1) sx# = sx# + sxv# sy# = sy# + syv# ;Cls Color 0,128,0 Oval 198,148,5,5,True Color 255,0,0 Oval sx# - 1,sy# - 1,3,3,True Flip Color 0,0,255 Oval sx# - 1,sy# - 1,3,3,True Wend FreeTimer timer End --v Uncomment the Cls to remove the trails. If anyone can find the perfect start point, gravity, or start speed, to make a perfect/near-perfect cirlce, then please let me know. Things you can change: Start position: Local sx# = 100,sy# = 150 --v Start velocity: Local sxv# = 0,syv# = Pi --v Gravity: sxv# = sxv# + (Cos#(ATan2#(sy# - 150,sx# - 200) + 180) * 0.01) syv# = syv# + (Sin#(ATan2#(sy# - 150,sx# - 200) + 180) * 0.01) --v ^Change the 0.01 to the gravity wanted. Must stay the same.^ \|edit\| Found an almost perfect orbit. Code updated. \|edit\| \|edit\| Found a better one. Interesting how Pi makes for a practically perfect cirlce... \|edit\|

Comments
Friday, 27 April 2007, 15:27
Jayenkai	Aww, man, I am gonna have SO much fun playing with this!! Thanks AppTitle "Orbit" Graphics 400,300,24,2 SetBuffer BackBuffer() Local sx# = 100,sy# = 150 Local sxv# = 0,syv# = Pi Dim Planet#(3,5) planet(1,1)=100 planet(2,1)=50 Planet(3,1)=0.1 planet(1,2)=300 planet(2,2)=50 Planet(3,2)=0.1 planet(1,3)=200 planet(2,3)=150 Planet(3,3)=0.1 planet(1,4)=300 planet(2,4)=250 Planet(3,4)=0.1 planet(1,5)=100 planet(2,5)=250 Planet(3,5)=0.1 timer = CreateTimer(50) While Not KeyHit(1) WaitTimer(timer) For n=1 To 5 sxv# = sxv# + (Cos#(ATan2#(sy# - planet(2,n),sx# - planet(1,n)) + 180) * planet(3,n)) syv# = syv# + (Sin#(ATan2#(sy# - planet(2,n),sx# - planet(1,n)) + 180) * planet(3,n)) Color 0,128,0 Oval planet(1,n)-2,planet(2,n)-2,5,5,True Next sx# = sx# + sxv# sy# = sy# + syv# ;Cls Color 255,0,0 Oval sx# - 1,sy# - 1,3,3,True Flip Color 0,0,255 Oval sx# - 1,sy# - 1,3,3,True Wend FreeTimer timer End --v WEEEEeeee!!!
Friday, 27 April 2007, 16:23
Blitz3Dman	It's an atom!! AppTitle "Orbit" Graphics 400,300,24,2 SetBuffer BackBuffer() Local sx# = 100,sy# = 150 Local sxv# = 0,syv# = Pi timer = CreateTimer(50) While Not KeyHit(1) WaitTimer(timer) sxv# = sxv# + (Cos#(ATan2#(sy# - 150,sx# - 200) + 180) * 0.3) syv# = syv# + (Sin#(ATan2#(sy# - 150,sx# - 200) + 180) * 0.3) sx# = sx# + sxv# sy# = sy# + syv# ;Cls Color 0,128,0 Oval 198,148,5,5,True Color 255,0,0 Oval sx# - 1,sy# - 1,3,3,True Flip Color 0,0,255 Oval sx# - 1,sy# - 1,3,3,True Wend FreeTimer timer End --v
Friday, 27 April 2007, 17:13
magicman	Thats pretty cool.
Friday, 27 April 2007, 17:17
Blitz3Dman	I'm making a really awesome one here, so hold on to your britches cause it's almost finished
Friday, 27 April 2007, 19:19
power mousey	ahhhhhhhh!!! thanks shroom. I'm beginning to see the light. Serious and True! whoa!! and wow!! light particles and radiated wavelengths over time and space...individual particles.photons grouped in packets forming waves as part of function in both time and space interacting with other objects and relativities. And in the grand expansion of the Cosmos relative to the temporal locales. And as to the point of initial singularity. True, I might point out. The singularity point. whoa I see, I see. but there is also the z factor in our 3 dimensions...and perhaps beyond as well. Utilizing and using one dimension above. Its been theorized around at least 7 and up to 24 dimensions. But thats a different story and topic for a later time. thank you for sharing this program. I have to ponder this. And check this out. Think it thru. good night,
Saturday, 28 April 2007, 00:02
shroom_monk	Glad you all enjoyed it!
Saturday, 28 April 2007, 00:10
shroom_monk	Here's one for ya, with 3 planets! AppTitle "Orbit 2" Graphics 400,300,24,2 SetBuffer BackBuffer() Local s1x# = 0,s1y# = 100 Local s1xv# = 0,s1yv# = 1 Local s2x# = 100,s2y# = 150 Local s2xv# = 0,s2yv# = 2 Local s3x# = 400,s3y# = 200 Local s3xv# = 3,s3yv# = 0 timer = CreateTimer(50) While Not KeyHit(1) WaitTimer(timer) s1xv# = s1xv# + (Cos#(ATan2#(s1y# - 150,s1x# - 200) + 180) * 0.1) s1yv# = s1yv# + (Sin#(ATan2#(s1y# - 150,s1x# - 200) + 180) * 0.1) s1x# = s1x# + s1xv# s1y# = s1y# + s1yv# s2xv# = s2xv# + (Cos#(ATan2#(s2y# - 150,s2x# - 200) + 180) * 0.1) s2yv# = s2yv# + (Sin#(ATan2#(s2y# - 150,s2x# - 200) + 180) * 0.1) s2x# = s2x# + s2xv# s2y# = s2y# + s2yv# s3xv# = s3xv# + (Cos#(ATan2#(s3y# - 150,s3x# - 200) + 180) * 0.1) s3yv# = s3yv# + (Sin#(ATan2#(s3y# - 150,s3x# - 200) + 180) * 0.1) s3x# = s3x# + s3xv# s3y# = s3y# + s3yv# ;Cls Color 0,128,0 Oval 198,148,5,5,True Color 255,0,0 Oval s1x# - 1,s1y# - 1,3,3,True Oval s2x# - 1,s2y# - 1,3,3,True Oval s3x# - 1,s3y# - 1,3,3,True Flip Color 0,0,255 Oval s1x# - 1,s1y# - 1,3,3,True Color 0,255,0 Oval s2x# - 1,s2y# - 1,3,3,True Color 255,255,0 Oval s3x# - 1,s3y# - 1,3,3,True Wend FreeTimer timer End --v How long do you think it'll be before they crash? And Jay, try making one, or even two, of your planets have a big gravity. The effects are interesting. \|edit\| Like this: \|edit\| planet(1,1)=100 planet(2,1)=50 Planet(3,1)=5 planet(1,2)=300 planet(2,2)=50 Planet(3,2)=0.1 planet(1,3)=200 planet(2,3)=150 Planet(3,3)=0.1 planet(1,4)=300 planet(2,4)=250 Planet(3,4)=5 planet(1,5)=100 planet(2,5)=250 Planet(3,5)=0.1 --v
Saturday, 28 April 2007, 02:42
TheMadProff	I love the first one of those two... And what about particle effects instead of the lines it creates... Thanks, thats really nice! Also, I love the second one, its weired!
Saturday, 28 April 2007, 02:50
shroom_monk	Well, of course, this is just a simple demo of the code. Feel free to take it, use it, and modify it how you want. Then post it here. I think one thing I will add is a setting where the gravity has less affect as the object gets further away, like in real life.
Saturday, 28 April 2007, 03:03
shroom_monk	Now, gravity effect gets less, the further away the block gets. If the effect gets to 0, the you've lost it! AppTitle "Orbit 3" Graphics 400,300,24,2 SetBuffer BackBuffer() Local sx# = 100,sy# = 150 Local sxv# = 0,syv# = Pi Local grav# = 0.1,spdis# = 0,gmul# = 1 timer = CreateTimer(50) While Not KeyHit(1) WaitTimer(timer) spdis# = Sqr#(((sx# - 200) ^ 2) + ((sy# - 150) ^ 2)) gmul# = 1 - (spdis# / 300) If gmul# < 0 Then gmul# = 0 sxv# = sxv# + (Cos#(ATan2#(sy# - 150,sx# - 200) + 180) * grav# * gmul#) syv# = syv# + (Sin#(ATan2#(sy# - 150,sx# - 200) + 180) * grav# * gmul#) sx# = sx# + sxv# sy# = sy# + syv# ;Cls Color 0,128,0 Oval 198,148,5,5,True Color 255,0,0 Oval sx# - 1,sy# - 1,3,3,True Color 255,255,255 Text 0,0,"Gravity Effect: " + (gmul# * 100) + "%" Text 0,10,"X Velocity: " + sxv# Text 0,20,"Y Velocity: " + syv# Flip Color 0,0,255 Oval sx# - 1,sy# - 1,3,3,True Color 0,0,0 Rect 0,0,200,35,True Wend FreeTimer timer End --v
Saturday, 28 April 2007, 06:37
Lineraf	That's pretty cool! Here's a BlitzMax version: AppTitle$ = "Orbit 3" Graphics 400,300 Local sx# = 100,sy# = 150 Local sxv# = 0,syv# = Pi Local grav# = 0.1,spdis# = 0,gmul# = 1 While Not KeyHit(KEY_ESCAPE) spdis# = Sqr(((sx# - 200) ^ 2) + ((sy# - 150) ^ 2)) gmul# = 1 - (spdis# / 300) If gmul# < 0 Then gmul# = 0 sxv# = sxv# + (Cos(ATan2(sy# - 150,sx# - 200) + 180) * grav# * gmul#) syv# = syv# + (Sin(ATan2(sy# - 150,sx# - 200) + 180) * grav# * gmul#) sx# = sx# + sxv# sy# = sy# + syv# 'Cls SetColor 0,128,0 DrawOval 198,148,5,5 SetColor 255,0,0 DrawOval sx# - 1,sy# - 1,3,3 SetColor 255,255,255 DrawText "Gravity Effect: " + (gmul# * 100) + "%",0,0 DrawText "X Velocity: " + sxv#,0,10 DrawText "Y Velocity: " + syv#,0,20 Flip(1) SetColor 0,0,255 DrawOval sx# - 1,sy# - 1,3,3 SetColor 0,0,0 DrawRect 0,0,200,35 Wend End --v
Saturday, 28 April 2007, 12:48
Blitz3Dman	Mwa ha ha!!! A self-sustaning universe featuring planets, moons/asteroids/space debris, explosions, and gravity based on shroom's code! Create planets by left mouseclicking! Create moons/asteroids/debris by right mouseclicking! AppTitle "Orbit" Graphics 800,600,24,2 SetBuffer BackBuffer() SeedRnd MilliSecs() Type moon Field sx#,sy#,sxv#,syv# End Type Type planet Field x#,y#,grav#,xv#,yv# End Type Type plosion Field x#,y#,xv#,yv#,time End Type For newmoon=1 To Rand(10,50) m.moon=New moon m\sx# = Rand(10,790):m\sy# = Rand(10,590) m\sxv# = Rnd(-0.1,0.1):m\syv# = Rnd(-0.1,0.1) Next For newplan=1 To Rand(10,20) p.planet=New planet p\x=Rand(100,700) p\y=Rand(100,500) p\grav=Rnd(0.001,0.003) p\xv=Rnd(-0.2,0.2) p\yv=Rnd(-0.2,0.2) Next timer = CreateTimer(50) While Not KeyHit(1) Cls WaitTimer(timer) If MouseDown(1) Then FlushMouse() p.planet=New planet p\x=MouseX() p\y=MouseY() p\xv=Rnd(-0.2,0.2) p\yv=Rnd(-0.2,0.2) p\grav=0.001 EndIf If MouseHit(2) Then FlushMouse() m.moon=New moon m\sx=MouseX() m\sy=MouseY() m\sxv#=Rnd(-0.1,0.1) m\syv#=Rnd(-0.1,0.1) EndIf For ex.plosion=Each plosion Color ex\time1.5,ex\time1.5,0 Rect ex\x,ex\y,2,2 ex\time=ex\time-1 ex\x=ex\x+ex\xv ex\y=ex\y+ex\yv For m.moon=Each moon If distance(ex\x,ex\y,m\sx,m\sy)<2 Then Delete m Next If ex\time<0 Then Delete ex Next For p.planet=Each planet goway=False p\x=p\x+p\xv p\y=p\y+p\yv If p\x>800 Then p\x=0 Else If p\x<0 Then p\x=800 If p\y>600 Then p\y=0 Else If p\y<0 Then p\y=600 Color 0,128,0 Oval p\x-1-(p\grav500),p\y-1-(p\grav500),3+(p\grav1000),3+(p\grav1000),True For op.planet=Each planet If p<>Null Then If p\x<>op\x And p\y<>op\y Then p\xv# = p\xv# + (Cos#(ATan2#(p\y# - op\y,p\x# - op\x) + 180) * (op\grav/p\grav/(Pi200))) p\yv# = p\yv# + (Sin#(ATan2#(p\y# - op\y,p\x# - op\x) + 180) (op\grav/p\grav/(Pi200))) If distance(p\x-1-(p\grav500),p\y-1-(p\grav500),op\x,op\y)<3+(p\grav1000) Then n.planet=New planet n\x=p\x n\y=p\y n\xv=p\xv/1.4 n\yv=p\yv/1.4 n\grav=p\grav+op\grav Delete op goway=True For boom=1 To Rand(3,12) ex.plosion=New plosion ex\x=p\x ex\y=p\y ex\xv=Rnd(-1.2,1.2) ex\yv=Rnd(-1.2,1.2) ex\time=Rand(50,150) Next EndIf EndIf EndIf Next For m.moon=Each moon If m<>Null And p<>Null Then spdis# = Sqr#(((m\sx# - p\x) ^ 2) + ((m\sy# - p\y) ^ 2)) gmul# = 1 - (spdis# / 300) m\sxv# = m\sxv# + (Cos#(ATan2#(m\sy# - p\y,m\sx# - p\x) + 180) * p\gravgmul) m\syv# = m\syv# + (Sin#(ATan2#(m\sy# - p\y,m\sx# - p\x) + 180) p\gravgmul) If distance(p\x-1-(p\grav500),p\y-1-(p\grav500),m\sx,m\sy)<3+(p\grav1000) Then For newmoon=1 To 2 m.moon=New moon m\sx=p\x+Rnd(-15,15) m\sy=p\y+Rnd(-15,15) m\sxv=Rnd(-0.1,0.1) m\syv=Rnd(-0.1,0.1) Next For boom=1 To Rand(1,2) ex.plosion=New plosion ex\x=m\sx ex\y=m\sy ex\xv=-m\sxv+Rnd(-0.2,0.2) ex\yv=-m\syv+Rnd(-0.2,0.2) ex\time=Rand(50,150) Next Delete m p\grav=p\grav-0.0015 If p\grav=<0 Then goway=True EndIf EndIf Next If p\grav>0.5 Or p\xv>50 Or p\yv>50 Then For boom=1 To Rand(p\grav/4,p\grav) ex.plosion=New plosion ex\x=p\x ex\y=p\y ex\xv=Rnd(-p\xv/3,p\xv/3) ex\yv=Rnd(-p\yv/3,p\yv/3) ex\time=Rand(50,150) Next For newmoon=1 To p\grav/10 m.moon=New moon m\sx=p\x m\sy=p\y m\sxv=Rnd(-p\xv/3,p\xv/3) m\syv=Rnd(-p\yv/3,p\yv/3) Next goway=True EndIf If goway=True Then Delete p Next For m.moon=Each moon m\sx# = m\sx# + m\sxv# m\sy# = m\sy# + m\syv# If Abs(m\sxv)>10 Or Abs(m\syv)>10 Then For boom=1 To Rand(10,18) ex.plosion=New plosion ex\x=m\sx ex\y=m\sy ex\xv=Rnd(-m\sxv/3,m\sxv/3) ex\yv=Rnd(-m\syv/3,m\syv/3) ex\time=Rand(50,150) Next Delete m Exit EndIf Color 128,128,128 Plot m\sx,m\sy If m\sx>800 Then m\sx=0 Else If m\sx<0 Then m\sx=800 If m\sy>600 Then m\sy=0 Else If m\sy<0 Then m\sy=600 For o.moon=Each moon If m<>Null Then If distance(m\sx,m\sy,o\sx,o\sy)<2 And m\sx<>o\sx And m\sy<>o\sy Then p.planet=New planet p\x=m\sx p\y=m\sy p\grav=Rnd(0.001,0.003) p\xv=Rnd(-((o\sxv+m\sxv)/2),((o\sxv+m\sxv)/2)) p\yv=Rnd(-((o\syv+m\syv)/2),((o\syv+m\syv)/2)) For boom=1 To Rnd(5,10) ex.plosion=New plosion ex\x=m\sx ex\y=m\sy ex\xv=Rnd(-1,1) ex\yv=Rnd(-1,1) ex\time=Rand(50,150) Next Delete m Delete o Exit EndIf Next Next Flip Wend FreeTimer timer End Function Distance(x, y, x2, y2) ;courtesy to Phoenix for function! x = Abs( x2 - x ) y = Abs( y2 - y ) Return Sqr((x * x) + (y * y)) End Function --v
Saturday, 28 April 2007, 18:27
JL235	I made a few changes to it. But I think I might have broken it. I'm pretty certain the numbers I'm using are Shrooms original numbers, but the planet movement I get is different. The changes I've made; added a Planet type and Planet functions (so it's easier to add more Planets), renamed 'sxv', 'sxy', 'sx' and 'sy', removed the timer and added delta time. But I don't understand why it's still not moving the way it was originally on Shrooms. Maybe me numbers are wrong, but I have discovered (by fluke) the perfect orbit. ;initialize AppTitle "Orbit" Graphics 400,300,24,2 SetBuffer BackBuffer() ; delta time Const PERFECT_DELTA# = 1000.0/60.0 Global delta# = 1 Global oldTime = MilliSecs() Function UpdateDeltaTime() delta# = (MilliSecs()-oldTime) / PERFECT_DELTA# End Function ; planet Type Planet Field x#, y# Field i#, j# End Type Local planet.Planet = NewPlanet(100, 150, 0, Pi) ; main loop While Not KeyHit(1) UpdatePlanet(planet) ;Cls ; draw sun Color 0,128,0 Oval 198,148,5,5,True ; draw planet Color 255,0,0 DrawPlanet(planet) Flip ; draw planet trail Color 0,0,255 DrawPlanet(planet) UpdateDeltaTime() Wend End Function NewPlanet.Planet(x#, y#, i#, j#) Local p.Planet p = New Planet p\x# = x# p\y# = y# p\i# = i# p\j# = j# Return p End Function Function UpdatePlanet(p.Planet) p\i# = p\i# + (Cos#(ATan2#(p\y# - 150,p\x# - 200) + 180) * 0.1);delta# p\j# = p\j# + (Sin#(ATan2#(p\y# - 150,p\x# - 200) + 180) 0.1);*delta# p\x# = p\x# + p\i# p\y# = p\y# + p\j# End Function Function DrawPlanet(p.planet) Oval p\x# -1, p\y# -1, 3, 3, True End Function --v
Sunday, 29 April 2007, 00:30
shroom_monk	@Blitz3DMan: Now that is good. Although, I think some of the collision detection needs checking, and it has a habit of crashing after a while (although I think that might just be my computer). @DiablosDevil: Your program is effectively the same as mine - same orbit, etc... because you commented out the delta, so it doesn't have any effect. Still pretty good, though. Has anyone found the perfect orbit for when gravity decreases as you move away (Orbit 3)?
Sunday, 29 April 2007, 06:09
Blitz3Dman	Although, I think some of the collision detection needs checking, and it has a habit of crashing after a while (although I think that might just be my computer). Yea, I know and it's not just your computer. It happens a lot when you go crazy creating planets. And this code is a nice bit faster modifying only 2 lines of code \|edit\| I fixed the collision detection. turns out I was detecting using the diameter instead of the radius. \|edit\| \|edit\| I found out that I forgot to set gmul at minimum 0 so I was having some odd problems... \|edit\| \|edit\| yet another edit: now planets use gmul and the collision still isn't fixed completely >.< \|edit\| AppTitle "Orbit" Graphics 800,600,24,2 SetBuffer BackBuffer() SeedRnd MilliSecs() Type moon Field sx#,sy#,sxv#,syv# End Type Type planet Field x#,y#,grav#,xv#,yv# End Type Type plosion Field x#,y#,xv#,yv#,time End Type For newmoon=1 To Rand(10,50) m.moon=New moon m\sx# = Rand(10,790):m\sy# = Rand(10,590) m\sxv# = Rnd(-0.1,0.1):m\syv# = Rnd(-0.1,0.1) Next For newplan=1 To Rand(10,20) p.planet=New planet p\x=Rand(100,700) p\y=Rand(100,500) p\grav=Rnd(0.001,0.003) p\xv=Rnd(-0.2,0.2) p\yv=Rnd(-0.2,0.2) Next timer = CreateTimer(50) While Not KeyHit(1) Cls WaitTimer(timer) If MouseDown(1) Then FlushMouse() p.planet=New planet p\x=MouseX() p\y=MouseY() p\xv=Rnd(-0.2,0.2) p\yv=Rnd(-0.2,0.2) p\grav=0.001 EndIf If MouseHit(2) Then FlushMouse() m.moon=New moon m\sx=MouseX() m\sy=MouseY() m\sxv#=Rnd(-0.1,0.1) m\syv#=Rnd(-0.1,0.1) EndIf If MouseDown(3) Then FlushMouse() m.moon=New moon m\sx=MouseX() m\sy=MouseY() m\sxv#=Rnd(-0.1,0.1) m\syv#=Rnd(-0.1,0.1) EndIf For ex.plosion=Each plosion Color ex\time1.5,ex\time1.5,0 Rect ex\x,ex\y,2,2 ex\time=ex\time-1 ex\x=ex\x+ex\xv ex\y=ex\y+ex\yv For m.moon=Each moon If distance(ex\x,ex\y,m\sx,m\sy)<2 Then Delete m Next If ex\time<0 Then Delete ex Next Color 0,128,0 For p.planet=Each planet goway=False p\x=p\x+p\xv p\y=p\y+p\yv If p\x>800 Then p\x=0 Else If p\x<0 Then p\x=800 If p\y>600 Then p\y=0 Else If p\y<0 Then p\y=600 Oval p\x-1-(p\grav500),p\y-1-(p\grav500),3+(p\grav1000),3+(p\grav1000),True For op.planet=Each planet If p<>Null Then If p\x<>op\x And p\y<>op\y Then spdis# = Sqr#(((p\x - op\x) ^ 2) + ((p\y - op\y) ^ 2)) gmul# = 1 - (spdis# / 300) If gmul<0 Then gmul=0 p\xv# = p\xv# + (Cos#(ATan2#(p\y# - op\y,p\x# - op\x) + 180) * (op\grav/p\gravgmul/(Pi200))) p\yv# = p\yv# + (Sin#(ATan2#(p\y# - op\y,p\x# - op\x) + 180) * (op\grav/p\gravgmul/(Pi200))) If distance(p\x-1-(p\grav500),p\y-1-(p\grav500),op\x,op\y)<3+(p\grav500) Then n.planet=New planet n\x=p\x n\y=p\y n\xv=p\xv/1.4 n\yv=p\yv/1.4 n\grav=p\grav+op\grav Delete op goway=True For boom=1 To Rand(3,12) ex.plosion=New plosion ex\x=p\x ex\y=p\y ex\xv=Rnd(-1.2,1.2) ex\yv=Rnd(-1.2,1.2) ex\time=Rand(50,150) Next EndIf EndIf EndIf Next For m.moon=Each moon If m<>Null And p<>Null Then spdis# = Sqr#(((m\sx# - p\x) ^ 2) + ((m\sy# - p\y) ^ 2)) gmul# = 1 - (spdis# / 300) If gmul<0 Then gmul=0 m\sxv# = m\sxv# + (Cos#(ATan2#(m\sy# - p\y,m\sx# - p\x) + 180) p\gravgmul) m\syv# = m\syv# + (Sin#(ATan2#(m\sy# - p\y,m\sx# - p\x) + 180) p\gravgmul) If distance(p\x-1-(p\grav500),p\y-1-(p\grav500),m\sx,m\sy)<3+(p\grav500) Then For newmoon=1 To 2 m.moon=New moon m\sx=p\x+Rnd(-15,15) m\sy=p\y+Rnd(-15,15) m\sxv=Rnd(-0.1,0.1) m\syv=Rnd(-0.1,0.1) Next For boom=1 To Rand(1,2) ex.plosion=New plosion ex\x=m\sx ex\y=m\sy ex\xv=-m\sxv+Rnd(-0.2,0.2) ex\yv=-m\syv+Rnd(-0.2,0.2) ex\time=Rand(50,150) Next Delete m p\grav=p\grav-0.0015 If p\grav=<0 Then goway=True EndIf EndIf Next If p\grav>0.5 Or p\xv>50 Or p\yv>50 Then For boom=1 To Rand(p\grav/4,p\grav) ex.plosion=New plosion ex\x=p\x ex\y=p\y ex\xv=Rnd(-p\xv/3,p\xv/3) ex\yv=Rnd(-p\yv/3,p\yv/3) ex\time=Rand(50,150) Next For newmoon=1 To p\grav/10 m.moon=New moon m\sx=p\x m\sy=p\y m\sxv=Rnd(-p\xv/3,p\xv/3) m\syv=Rnd(-p\yv/3,p\yv/3) Next goway=True EndIf If goway=True Then Delete p Next Color 128,128,128 For m.moon=Each moon m\sx# = m\sx# + m\sxv# m\sy# = m\sy# + m\syv# If Abs(m\sxv)>10 Or Abs(m\syv)>10 Then For boom=1 To Rand(10,18) ex.plosion=New plosion ex\x=m\sx ex\y=m\sy ex\xv=Rnd(-m\sxv/3,m\sxv/3) ex\yv=Rnd(-m\syv/3,m\syv/3) ex\time=Rand(50,150) Next Delete m Exit EndIf Plot m\sx,m\sy If m\sx>800 Then m\sx=0 Else If m\sx<0 Then m\sx=800 If m\sy>600 Then m\sy=0 Else If m\sy<0 Then m\sy=600 For o.moon=Each moon If m<>Null Then If distance(m\sx,m\sy,o\sx,o\sy)<2 And m\sx<>o\sx And m\sy<>o\sy Then p.planet=New planet p\x=m\sx p\y=m\sy p\grav=Rnd(0.001,0.003) p\xv=Rnd(-((o\sxv+m\sxv)/2),((o\sxv+m\sxv)/2)) p\yv=Rnd(-((o\syv+m\syv)/2),((o\syv+m\syv)/2)) For boom=1 To Rnd(5,10) ex.plosion=New plosion ex\x=m\sx ex\y=m\sy ex\xv=Rnd(-1,1) ex\yv=Rnd(-1,1) ex\time=Rand(50,150) Next Delete m Delete o Exit EndIf Next Next Flip Wend FreeTimer timer End Function Distance(x, y, x2, y2) ;courtesy to Phoenix for function! x = Abs( x2 - x ) y = Abs( y2 - y ) Return Sqr((x * x) + (y * y)) End Function --v
Sunday, 29 April 2007, 10:32
JL235	Dam, I must have been comparing my code to someone elses. I swear last night I wasn't getting a perfect orbit with Shrooms code.
Sunday, 29 April 2007, 11:53
shroom_monk	Diablos: You get a perfect orbit with orbit 1 program only. My orbit's 2 and 3 don't have perfect orbit.

Orbital Physics

Comments