#include <cstdlib>
#include <cmath>
#include <iostream>
#include <stdio.h>
#include "ext_value.h"
using namespace std;

// для упрощения переносимости явно включил все нужные файлы
// более "профессиональным" было бы написать сборку проекта из нескольких файлов
// (в большом проекте сборка из нескольких файлов действительно удобнее)
// в представленном виде я не ожидаю никаких проблем ни у кого ни в каких компиляторах 
#include "fcn-ext.cpp"
#include "d87-ext.cpp"

int main()
{  
int i,n=6;
// t0 - момент начала движения
// x0,y0,z0 - координаты
// vx0,vy0,vz0 - скорости в момент выстрела
double t=0,x0,vx0,y0,vy0,z0,vz0,tend=1800,eps=1.e-12,hmax=600,h=100,mu=3.9860044e5,Vkr;

x0=6578.250; y0=0; z0=0; vx0=0; Vkr=sqrt(mu/x0); vy0=1.05*Vkr*cos(0.9); vz0=1.05*Vkr*cos(0.9);

ext_value2 y[n];
for(i=0;i<n;i++)  y[i].set_dim(n); // установка размерности
y[0].set(x0); y[1].set(y0); y[2].set(z0); y[3].set(vx0); y[4].set(vy0); y[5].set(vz0);
for(i=0;i<n;i++)  y[i].set(1.0,i); 

// в отличие от типа double
// значение всех производных получается одним вычислением с использованием соответствующего типа данных:
	dopri8(n,t,y,tend,eps,hmax,h);

// здесь вывод результатов
    printf("Получен ответ:\n");
	printf(" x0=%lf\n",y[0].get() );
	printf(" y0=%lf\n",y[1].get() );
	printf(" z0=%lf\n",y[2].get() );
	printf("vx0=%lf\n",y[3].get() );
	printf("vy0=%lf\n",y[4].get() );
	printf("vz0=%lf\n",y[5].get() );
    printf("Производная (ext_value):\n");
    for(i=0;i<n;i++) printf("d2(%d)=%lg\n",i,y[i].get(1,3) );

return 0;
}