# define main function
main := proc()

    print("Maple output");
    print("Chapter 3");
    prob05();
    prob05B(); # general solution
end;

prob05  := proc()
local a, t1,t2, v1,v2, eq1, eq2, eq3, eq4,eq5, eq6, eq7;

    print("Problem 5");
    
    eq1 := 0.0 = 4.0 -16.0 * t1 * t1;
    eq2 := v1 = -32.0 * t1;
    eq3 := 3.0 = 0.0 + v2 * t2 - 16.0 * t2 * t2;
    eq4 := 0.0 = v2 - 32.0 * t2;
    eq5 := a = (v2 - v1) / 0.010;
    eq6 := t1 > 0.0;
    eq7 := t2 > 0.0;

    print(solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7] , [t1,t2,v1,v2,a]));
end;
prob05B  := proc()
local a, t1,t2, v1,v2, y0, y2, g, deltat, eq1, eq2, eq3, eq4,eq5, eq6, eq7;

    print("Problem 5B");
    
    eq1 := 0.0 = y0 - 16.0 * t1 * t1;
    eq2 := v1 = - 32.0 * t1;
    eq3 := y2 = 0.0 + v2 * t2 - 16.0 * t2 * t2;
    eq4 := 0.0 = v2 - 32.0 * t2;
    eq5 := a = (v2 - v1) / 0.010;
    eq6 := t1 > 0.0;
    eq7 := t2 > 0.0;

    print(solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7] , [t1,t2,v1,v2,a]));
end;

# call main
main();