Das deutsche Python-Forum

Hallo,

ich belege dieses Semester Numerik und tue mich leider mit den Programmierübungen sehr schwer. Daher benötige ich etwas Hilfe.

Ich soll eine Funktion Function(x,f) erstellen, die die Matrix D aus R^(n+1,n+1) mit

D(i,k)= {[xi,...,xi+k]f ; k=0,...n, i=0,...,n-k.............0 ; sonst

folgendes habe ich bisher, bin mir aber nicht sicher:

import numpy as np
import matplotlib.pyplot as plt

def divDiff(x,f):
n = len(x)-1
Delta = np.zeros((n+1,n+1))

for i in range(0,n+1):
Delta[i,0] = f

for k in range(1,n+1):
for i in range (0, n-k):
Delta[i,k] = (Delta[i+1, k-1]-Delta[i,k-1])/(x[i+k]-x)

for i in range(0,n-k):
Delta[i,0] = f

return Delta

Servus,

Ich verstehe leider nicht ganz was genau die Funktion tun soll?
Außerdem wäre es schön, wenn du deinen Code in entsprechende Tags packst, dann ist das auch schön formatiert .

Warum ziehst Du von n 1 ab, wenn Du überall, wo Du n benutzt wieder +1 draufaddierst?
Aus Deiner Beschreibung werde ich nicht schlau, was Du denn nun für Matrizen bilden sollst. Hast Du Beispiele?