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

/* [wxMaxima: title   start ]
Parabel: maximale Krümmung
   [wxMaxima: title   end   ] */

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

/* [wxMaxima: comment start ]
******************************************
Dokumentnummer: D1036
Quelle: Johann Weilharter, D1035
******************************************

Es ist zu vermuten, dass eine Parabel ihre maximale
Krümmung im Scheitelpunkt hat.
   [wxMaxima: comment end   ] */

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

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

/* [wxMaxima: input   start ] */
f:read("Gib eine quadratische Funktion y=a*x**2+b*x+c ein");
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
f:''rhs(f);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ab:diff(f,x);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ab2:diff(f,x,2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
k:ab2/sqrt((1+ab**2)**3);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
abk:diff(k,x);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
m:solve(abk=0,x);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
xm:ev(x,m[1]);
/* [wxMaxima: input   end   ] */

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

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

/* [wxMaxima: input   start ] */
ks:ev(k,x=xs);
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
if xs=xm then print("Die maximale Krümmung einer Parabel ist im Scheitelpunkt")$
/* [wxMaxima: input   end   ] */

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