function dreikoerper_dopri
%Berechnung des Arenstorf-Orbits mit Dopri 5(4)-Verfahren


%Zeitraum festlegen
T = 17.06521656015796255889
% Anfangsbedingung setzen
%z_init = [0.994 0 0 -2.0015851063790825];
%z_init = [0.994 0 0 -2.01];
%z_init = [0.994 0 0 -2.02];
%z_init = [0.994 0 0 -2.03];
z_init = [0.994 0 0 -2.031732630];

z = zeros(1,4);
z(1,:) = z_init;
z_hat = zeros(1,4);
z_hat(1,:) = z_init;
z_old = z_init;
err = 0;
TOL = 10^(-6);
q = 4;
p=1/(q+1);
h = 0.1;
t_k = 0;
erg = [];
k=0;

%Runge-Kutta Koeffizienten festlegen
b = [35/384;0;500/1113;125/192;-2187/6784;11/84;0];
b_hat = [5179/57600; 0; 7571/16695;393/640;-92097/339200;187/2100;1/40];
c = [ 0; 1/5;3/10;4/5; 8/9; 1; 1];
A = [0 0 0 0 0 0;...
     1/5 0 0 0 0 0;...
     3/40 9/40 0 0 0 0;...
     44/45 -56/15 32/9 0 0 0;...
     19372/6561 -25369/2187 64448/6561 -212/729 0 0;...
     9017/3168 -355/33 46732/5247 49/176 -5103/18656 0;...
     35/384 0 500/1113 125/192 -2187/6784 11/84];

end