% testHW6

fprintf('\nTest Begin\n')

V0 = 50; h0 = 3;
fprintf('\nInput: V0 = %g, h0 = %g\n',V0,h0);
if exist('mytent1');
    [S,x,y,z,v] = mytent1(V0,h0);
    fprintf('\nMy solution: S = %g\n',S);
    fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);
end
[S,x,y,z,v] = tent1(V0,h0);
fprintf('\nYZ solution: S = %g\n',S);
fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);

S0 = S;
fprintf('\nInput: S0 = %g, h0 = %g\n',S0,h0);
if exist('mytent2');
    [V,x,y,z,v] = mytent2(S0,h0);
    fprintf('\nMy solution: V = %g\n',V);
    fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);
end
[V,x,y,z,v] = tent2(S0,h0);
fprintf('\nYZ solution: V = %g\n',V);
fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);


S0 = 80; h0 = 2;
fprintf('\nInput: S0 = %g, h0 = %g\n',S0,h0);
if exist('mytent2');
    [V,x,y,z,v] = mytent2(S0,h0);
    fprintf('\nMy solution: V = %g\n',V);
    fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);
end
[V,x,y,z,v] = tent2(S0,h0);
fprintf('\nYZ solution: V = %g\n',V);
fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);

V0 = V;
fprintf('\nInput: V0 = %g, h0 = %g\n',V0,h0);
if exist('mytent1');
    [S,x,y,z,v] = mytent1(V0,h0);
    fprintf('\nMy solution: S = %g\n',S);
    fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);
end
[S,x,y,z,v] = tent1(V0,h0);
fprintf('\nYZ solution: S = %g\n',S);
fprintf('[x,y,z,v] = [%9.4e, %9.4e, %9.4e, %9.4e, %9.4e]\n\n',x,y,z,v);