Funktion erstellen

[GP14]

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

[NPC]

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 .

[Sirius3]

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?