|
 |
Inspired by Draco's n-body system anims, I did a quick anim last night of such a
system using my user defined force code. This time, all of the particles are
attracted to one another, with the force proportional to (m1*m2)/(distance^2). I
only used quick and dirty figures to get a result close to what I wanted. The sun
weighs about 20000 times more that the planets, so it pretty much stays put in
the centre. The planets orbiting an the x-z plane were made to have an almost
circular orbit, whilst the planet in the y-z plane has some ecentricity.
Comments, complaints, etc. always welcome.
MJL
--
prism{0,.1,30#local I=1;#while(I<30)#local B=asc(substr(// Mark James Lewin
"#K?U_u`V[RG>3<9DGPL.0EObkcPF'",I,1))-33;<div(B,10)-4mod(B,10)+5*div(I,21)-
6>#local I=I+1;#end,-4pigment{rgb 9}rotate-x*90translate 15*z}//POV-Ray 3.5
Post a reply to this message
Attachments:
Download 'solar.mpg' (280 KB)
|
|
|
|
> Inspired by Draco's n-body system anims, I did a quick anim last night of such
a
> system using my user defined force code. This time, all of the particles are
> attracted to one another, with the force proportional to (m1*m2)/(distance^2).
I
> only used quick and dirty figures to get a result close to what I wanted. The
sun
> weighs about 20000 times more that the planets, so it pretty much stays put in
> the centre. The planets orbiting an the x-z plane were made to have an almost
> circular orbit, whilst the planet in the y-z plane has some ecentricity.
> Comments, complaints, etc. always welcome.
Just let it run for a couple of minutes, I'm sure there will be not much left of
your solar system :-/
Nice animation btw!
cu!
--
camera{location-z*3}#macro G(b,e)b+(e-b)*(C/50)#end#macro L(b,e,k,l)#local C=0
;#while(C<50)sphere{G(b,e),.1pigment{rgb G(k,l)}finish{ambient 1}}#local C=C+1
;#end#end L(y-x,y,x,x+y)L(y,-x-y,x+y,y)L(-x-y,-y,y,y+z)L(-y,y,y+z,x+y)L(0,x+y,
<.5,1,.5>,x)L(0,x-y,<.5,1,.5>,x) // ZK http://www.povplace.be.tf
Post a reply to this message
|
|
|
|
Zeger Knaepen wrote:
> Nice animation btw!
Thanks!
MJL
--
text{ttf"timrom.ttf"concat(#local O=1;#while(O<7)chr(val(substr(concat(#local Q=
1;#while(Q<7)str(asc(substr("???<?>",Q,1))-56,0,0),#local Q=Q+1;#end""),O,2))),#
local O=O+2;#end"").1,0pigment{rgb 9}translate-<1,.3,-2>} // MJL
Post a reply to this message
|
|
