Je bloque sur un programme Python depuis quelques jours.
Tout d'abord, dans ce qui va suivre, $ h=\dfrac{1}{n+1} $ et $ x_i = ih $.
Etant donnée une fonction f(x,y) donnée par l'utilisateur, faire un programme calculant la matrice $ (b_{i,j})_{1 \leq i,j \leq n} $ des valeurs $ b_{i,j} = f(x_i,x_j) $, puis le vecteur colonne $ b $ correspondant, dans $ \mathbb{R}^N $. On pourra faire une fonction suivant le format bvector(f,n) renvoyant le vecteur colonne $ b $ désiré.
J'ai commencé ainsi :
Code : Tout sélectionner
def bvector(f,n):
h=1/(n+1)
x=np.zeros(n)
for i in range(n):
x[i]=i*h
b=np.zeros((n,n))
for i in range(n):
for j in range(n):
b[i,j]=f(x[i],x[j])
return b.reshape
def f(x,y):
f=x+y
return f
print(bvector(f(3,4),2))