function nma_solveProb_43
%solves HW problem 4.3

r0_vec = input('Enter r0 >');
v0_vec = input('Enter v0 >');
nu0    = input('Enter initial angle nu (in degrees) >');
nu     = input('Enter final angle nu (in degrees) >');
mu     = input('Enter mu (gravitational parameter, use 1 for earth) >');

nu0 = nu0 * (pi/180);
nu  = nu  * (pi/180);

 
tof=nma_getDeltaTimeFromDeltaNu(r0_vec,v0_vec,nu0,nu,mu)
[r_vec,v_vec] = nma_moveProbe(r0_vec,v0_vec,tof,mu);

fprintf('r0 = %6.6f I + %6.6f J + %6.6f K\n',r0_vec(1),r0_vec(2),r0_vec(3));
fprintf('v0 = %6.6f I + %6.6f J + %6.6f K\n',v0_vec(1),v0_vec(2),v0_vec(3));

fprintf('\nProbe is moving .....\n\n');

fprintf('r = %6.6f I + %6.6f J + %6.6f K\n',r_vec(1),r_vec(2),r_vec(3));
fprintf('v = %6.6f I + %6.6f J + %6.6f K\n',v_vec(1),v_vec(2),v_vec(3));