/Users/nathalie/Desktop/Capture d’écran 2013-05-03 à 17.24.11.png

V(h):=%pi*(R^2-h^2)*(R+h)/3;

ratsimp(%);

diff(V(h),h,1);

ratsimp(%);

expand(%o3);

declare(h,mainvar);

expand(%o3);

delta:%pi^2*((-2*R/3)^2-4*(-1)*R^2/3);

h1:(2*%pi*R/3-4*%pi*R/3)/(2*(-%pi));

h2:(2*%pi*R/3+4*%pi*R/3)/(2*(-%pi));

solve([%o3], [h]);

V(R/3);

R:1;

plot2d([V(h)], [h,0,1],
 [gnuplot_preamble, "set size ratio 1; set zeroaxis; set grid;"])$

plot2d([V(h)], [h,0,1],
 [plot_format, gnuplot],
 [gnuplot_preamble, "set size ratio 1; set zeroaxis; set grid;"], [gnuplot_term, ps],
 [gnuplot_out_file, "/Users/nathalie/Desktop/PS/plot2d-Ex1.eps"])$

kill(R);