/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/

/* [wxMaxima: title   start ]
Zylinder mit minimaler Oberfläche
   [wxMaxima: title   end   ] */

/* [wxMaxima: section start ]
Problembeschreibung
   [wxMaxima: section end   ] */

/* [wxMaxima: comment start ]
******************************************
Dokumentnummer: D1034
Quelle: Johann Weilharter, Analysis
******************************************

Gegeben ist das Volumen eines Zylinders. Wie müssen
die Abmeldungen für r und h gewählt werden, damit die
Oberfläche des Zylinders minimal wird?
   [wxMaxima: comment end   ] */

/* [wxMaxima: section start ]
Problemlösung
   [wxMaxima: section end   ] */

/* [wxMaxima: comment start ]
EINGABE
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
V:read("Gib das Volumen des Zylinders ein");
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
VERARBEITUNG
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
NB:V=r**2*%pi*h;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
l:solve(NB,h);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
O:2*r**2*%pi+2*r*%pi*h;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
O:O,l[1];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ab:diff(O,r);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
extr:solve(ab=0,r);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
r_ext:r,extr[3];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
h:ev(h,l[1]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
h_ext:ev(h,extr[3]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
r_ext:r_ext,numer;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
h_ext:h_ext,numer;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
r_ext:floor(r_ext*100+0.5)/100.0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
h_ext:floor(h_ext*100+0.5)/100.0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ab2:diff(ab,r);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ab2:ev(ab2,extr[3]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
AUSGABE
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
if ab2>0 then print("Das Minimum")$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
if 2*r_ext=h_ext then print("Gleichseitiger Zylinder!")$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Radius =",r_ext)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Höhe =",h_ext)$
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$