rpt2func(x)

function fx = rpt2func(x)

% rz objective function

% required by repeat2.m

% Orbital Mechanics with Matlab

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

global tiwrk yi yfinal

% number of differential equations

neq = 6;

% step size guess

h = 60;

% integrate from ttmp to requested time = x

tetol = 1.0e-8;

yfinal = rkf78('j2eqm', neq, tiwrk, x, h, tetol, yi);

% compute current value of z component of position vector

fx = yfinal(3);