2- and 3-Body Demo

Two and three body problems.

function f2body (x,y,m1,m2)

  y=f(x,y) for the differential equations

function test2body (y1,m1=1,m2=1,n=100,tmax=360,fixed=0,h=1)

  2-body demo
 
  Intial y is a row vector with two positions and two speed vectors
  in the form [x1,y1,x2,y2,vx1,vy1,vx2,vy2]. m1 and m2 are the masses
  of the two bodies. 
 
  For slower animation increase h. For more accuracy increase n. The
  simulation runs until one body leaves the definition area (unless
  fixed is true), or tmax is reached.
 
  >fixed fixes the display at the first mass, which is positioned
  in (0,0).

Functions

function f3body (x,y,m1,m2,m3)

  y=f(x,y) for the differential equations

function test3body (y1,m1=1,m2=1,m3=1,n=100,fixed=0,tmax=360,h=1)

  3-body demo
 
  Works like test2body, but takes three positions and speed vectors,
  and three masses.

Set Planets

function setplanets

  Set the planets and initial vectors with the mouse.
 
  3-body which sets the positions with the mouse.

Examples