**TI92P*mainmaRnewton_i°¥ZV() Prgm ClrIO Disp "Newton'sches" Disp "Interpolationspolynom" Disp "","","" Disp "Autor: Marc Niggemann" Disp "Homepage: www.m-niggemann.com" Disp "Version: 3.01 Jan 2003" Pause ClrIO Local p,a,i,j Input "Anzahl der Wertepaare?",p ClrIO newMat(p,2)a For i,1,p,1 Input "x"&string(i)&"=?",j ja[i,1] Input "y"&string(i)&"=?",j ja[i,2] ClrIO EndFor Disp "xi yi" Disp a Pause ClrIO ©--------Differenzen Disp "Differenzen" Local l,lr,ll For l,1,p-1,1 newMat(p-l,1)lr For i,1,p-l,1 a[i,1]lr[i,1] EndFor newMat(p-l,1)ll For i,p,l+1,­1 ©l+1 wegen den grenzen a[i,1]ll[i-l,1] EndFor Disp string(l)&". linke Spalte" ll-lr#("l"&string(l)) ©speicherung Disp #("l"&string(l)) Pause EndFor ©----------- dividierten Differenzen ClrIO Disp "dividierten Differenzen" Local r,rz,px,pxx,hil,indi 0pxx 0indi For l,1,p-1,1 newMat(p-l,1)r For i,1,p-l,1 If l=1 Then (a[i+1,2]-a[i,2])/(#("l"&string(l))[i,1])r[i,1] Else (rz[i+1,1]-rz[i,1])/(#("l"&string(l))[i,1])r[i,1] EndIf EndFor Disp string(l)&". rechte Spalte" Disp r rrz Pause 1hil For j,1,l,1 hil*(x-a[j,1])hil EndFor rz[1,1]*hilpx Disp "Teil der Summe von Px",px Pause pxx+pxpxx EndFor Disp "y0 + die Summe ",a[1,2] Pause pxx+a[1,2]pxx ClrIO Disp "p(x)=" Disp pxx ©--------loeschung For l,1,p-1,1 DelVar #("l"&string(l)) EndFor EndPrgmLäåÜ—