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))