Program Polinom Newton dengan Python 3.2

from numpy import*

#Inisiasi Matriks (x,y)
x=array([[80.0],[78.0],[76.0],[74.0],[72.0],[70.0],[68.0]])
y=array([[30.65],[29.46],[28.23],[26.98],[25.03],[23.87],[22.23]])
n=len(x)
ST=zeros((n,n))

#Mencari Selisih Terbagi
for k in range (0,n):
    ST[k][0]=y[k]
for k in range (1,n):
    for i in range (0,n-k):
        ST[i][k]=(ST[i+1][k-1]-ST[i][k-1])/(x[i+k]-x[i])
print(‘Tabel Selisih Terbagi’)
print(ST)

#Proses
ulang=”yes”
while ulang==”yes”:
    if ulang==”yes”:        
        print(‘—————————————————————‘)
        a=float(input(‘Masukkan x yang ingin dicari nilai y nya: ‘))
        jumlah=ST[0,0]
        for i in range (1,n):
            suku=ST[0][i]
            for k in range (0,i):
                suku=suku*(a-x[k])
            jumlah=jumlah+suku
        print(‘Dengan polinomial Newton: f(‘,a,’) = ‘,”%.2f”%jumlah)

        ulang=input(“Coba lagi? (yes/no)”)