% script to generate the data for assignment two
% data for part 1
diary sol_two_data
format short e;
L = 3.75
t1 = 1.2457
[xt1,vt1] = cork(t1,0,0);
err1 = xt1 - L
% data for part 2
t2 = 1.2494
[xt2,vt2] = cork1(t2,0,0);
err2 = xt2 - 3.76
dt = t2-t1
% data for part 3
[xt3,vt3] = cork(t1,0.0,1.0e-4);
err3 = xt3 - L
[xt4,vt4] = cork(t1,1.0e-4,0.0);
err4 = xt4 - L
diary off