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

/* [wxMaxima: title   start ]
Nullstellen einer Polynomfunktion
   [wxMaxima: title   end   ] */

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

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

Die Nullstellen einer Polynomfunktion sind zu bestimmen.
   [wxMaxima: comment end   ] */

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

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

/* [wxMaxima: input   start ] */
f:read("Gib die Funktion y=f(x) ein");
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
g:rhs(f)=0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
l:realroots(g);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
n:length(l);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
X:makelist(ev(x,l[i]),i,1,n);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
X:makelist(floor(X[i]*1000+0.5)/1000.0,i,1,n);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
P:makelist([X[i],0],i,1,n);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Nullstellen:makelist(N[i]=P[i],i,1,n);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Nullstellen:transpose(Nullstellen);
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
print("Die Nullstellen sind",Nullstellen)$
/* [wxMaxima: input   end   ] */

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