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

/* [wxMaxima: title   start ]
Das Newtonverfahren
   [wxMaxima: title   end   ] */

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

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

Für eine Gleichung soll eine Folge von 10 Näherungslösungen
nach dem Newtonschen Verfahren für einen gegebenen
Startwert ermittelt werden.
   [wxMaxima: comment end   ] */

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

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

/* [wxMaxima: input   start ] */
g:read("Gib eine Gleichung ein");
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
x[0]:read("Gib den Startwert ein");
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
f:lhs(g)-rhs(g);
/* [wxMaxima: input   end   ] */

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

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

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

/* [wxMaxima: input   start ] */
a(x):=''ab;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
g(n):=x[n]:x[n-1]-f(x[n-1])/a(x[n-1]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
X:makelist(g(n),n,1,10),numer;
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
X:makelist(a[i]=X[i],i,1,10);
/* [wxMaxima: input   end   ] */

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

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

/* [wxMaxima: input   start ] */
print("Folge der Näherungslösungen =",X)$
/* [wxMaxima: input   end   ] */

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