• Support PF! Buy your school textbooks, materials and every day products Here!

C code, restricted three body

  • Thread starter Pixter
  • Start date
  • #1
30
0

Homework Statement




http://www2.warwick.ac.uk/fac/sci/physics/teach/module_home/px270/projects/c4s_projects07.pdf [Broken]

it's the first example

Homework Equations


see above


The Attempt at a Solution



right i shouldn't have any probs with the programming but i need to understand the calcs that are put on the site..

as i understand i'm suppose to the first eqn to find a,, then usign the second eqn to find rn+1.. this is done in components so for x and y..

but i tried it for a third mass that is at x=0 and y=positive integer(both the stationary massses are of equal magnitude).. what should happen is that it should be attracted and just go towards negative integers but I get acceleration in x direction.. ((note the printf function are just there for me to see some stuff))..

could someone maybe see where i've gone wrong and maybe give me som pointers:

here is the code:

#include <stdlib.h>
#include <stdio.h>
#include <math.h>


main(int argc, char* argv[])
{
float x[101], y[101], mneg, mnegm, mpos, rpnx, rpny, rppx, rppy, ax[101], ay[101], xzero, yzero, vxzero, vyzero, dt;
int n;
FILE *output;

output = fopen("data.dat", "w");

sscanf(argv[1], "%f", &mneg);
sscanf(argv[2], "%f", &mpos);
sscanf(argv[3], "%f", &xzero);
sscanf(argv[4], "%f", &yzero);
sscanf(argv[5], "%f", &vxzero);
sscanf(argv[6], "%f", &vyzero);



dt = 0.005;
x[0] = xzero;
y[0] = yzero;
x[1] = xzero + vxzero*dt;
y[1] =yzero + vyzero*dt;
rpnx= -1;
rpny= 0;
rppx= 1;
rppy= 0;
printf("%f %f\n", x[1], y[1]);

for(n=1; n < 100; ++n)
{
ax[n]=-mneg*((x[n]+rpnx)/(sqrt(pow((x[n]+rpnx),6))))-mpos*((x[n]-rppx)/(sqrt(pow((x[n]-rppx),6))));
ay[n]=-mneg*((y[n]+rpny)/(sqrt(pow((y[n]+rpny),6))))-mpos*((y[n]-rppy)/(sqrt(pow((y[n]-rppy),6))));

x[n+1]=((ax[n])*pow(dt,2))-(x[n-1])+2*(x[n]);
y[n+1]=((ay[n])*pow(dt,2))-(y[n-1])+2*(y[n]);
fprintf(output, "%f %f\n", x[n+1],y[n+1]);

}
printf("%f %f", ax[1], ay[1]);

fclose(output);
}


edit: note that the first mass is inputted as a "negative mass" for the equation to work..
also a wird problem if setting inital velocities to zero then the third mass doesn't move at ll which it should
 
Last edited by a moderator:

Answers and Replies

  • #2
Nothing appears when I open the web page...
 
  • #3
AlephZero
Science Advisor
Homework Helper
6,994
291
The attachment worked for me though it was a bit slow to open.

You seem to have misunderstood equation (1) in the attachment. The |r+r_p|^3 and |r-r_p|^3 terms are not vectors. |r+r_p| is the scalar length of the vector r+r_p.
 
Last edited:

Related Threads on C code, restricted three body

Replies
1
Views
16K
  • Last Post
Replies
3
Views
5K
  • Last Post
Replies
14
Views
880
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
9
Views
742
Replies
7
Views
1K
  • Last Post
Replies
3
Views
1K
Replies
9
Views
4K
Replies
4
Views
2K
Top