# Ejercicio 11.20

kill(all)$

## Primera parte

Designamos: x, radio de la base de la lata (forma de cilíndro recto circular)
                       y, altura de la lata
V = 330 cl

Superficie y volumen de la lata: 

S(x,y):=2*%pi*x*y + 2*(%pi*x^2);

V(x,y):=%pi*x^2*y;

Función de Lagrange:

L(x,y,L) := S(x,y) + L*(V(x,y) - V) ;

Derivadas parciales de la función de Lagrange:

Eq1 : diff(L(x,y,L),x)=0;  
Eq2 : diff(L(x,y,L),y)=0;  
Eq3 : V(x,y) - V = 0;
algsys([Eq1, Eq2, Eq3], [x,y,L]);

La única solución con sentido es la primera

x1 = x1:V^(1/3)/(2^(1/3)*%pi^(1/3));
y1 = y1:(2^(2/3)*V^(1/3))/%pi^(1/3);

Superficie mínima:

S(x1,y1) ;

## Segunda parte

V = V:330;

Eq1: V = %pi*x^2*y;
solve (Eq1, y);

define (f(x) , S(x, 330/(%pi*x^2) ));

Dimensiones de la lata de superficie mínima:

x1,  numer;
y1, numer;

Gráfica de la superficie en función del radio de la base de la lata

wxplot2d([f(x)], [x,1,10], [y, 0, 900], [ylabel, ""] , [legend, "f(x)"] , [style, [lines, 2]] ,
 [gnuplot_preamble, "set grid;"])$