Would you be so kind to change algorithm for solving this system
d2_x/dt_2=-GMx/(x^2+y^2)^(3/2)
d2_y/dy_2=-GMy/(x^2+y^2)^(3/2)
algorithm .It is called Runge-Cute Method
#include <stdio.h>
#include <stdafx.h>
#include <math.h>
#include <iostream>
using namespace std;
double f(int i,double x,double y[4]){
switch (i){
case 1:return y[4];break;
case 2...
2008-03-17 18:42:47 UTC in Solar Model
