{
gROOT->Reset();
const Int_t nmax=400;
Float_t t1[nmax], x1[nmax], y1[nmax], z1[nmax], v1[nmax];
Int_t ntimes, m = 0;
FILE * pFile;
pFile = fopen("../../work/ode45.dat","r");
for(Int_t i=0; i<nmax;i++){
    fscanf(pFile,"%f %f %f %f %f",&t1[i],&x1[i],&y1[i],&z1[i],&v1[i]);
    ntimes++;
}

TCanvas *c1 = new TCanvas("c1");
c1->Divide(2,1);
TGraph2D *gr1 = new TGraph2D(nmax,x1,y1,z1);
TGraph *gr = new TGraph(nmax,t1,v1);
gr->GetXaxis()->SetTitle("t");
gr->GetYaxis()->SetTitle("v");

c1->cd(1);
gr->Draw("AP");
c1->cd(2);
gr1->Draw("CP");
c1->Update();
c1->Modified();
fclose(pFile);
}