Загрузка данных


A=[[11,12,13],[21,22,23],[31,32,33]]
A

import numpy as np
A=np.array(A)
A[:,0]

# решение краевой задачи методом конечной разницы на пайтон
n=5
x1=0
xn=1
h=(xn-x1)/(n-1)
X=[i*h for i in range(n)]
P=[23/1+4*x*(1-x) for x in X]
F=[-1.88*(2-P[i]*X[i]+(1-X[i])) for i in range(n)]
P=[round(p,3) for p in P]
import numpy as np
A=np.zeros((n,n))
for i in range(n):
    for j in range(n):
        if i==j:
            A[i][j]=-2/h**2+P[i]
        if i==j-1 or i==j+1:
            A[i][j]=1/h**2
#учитываем краевые условия
F[0]=0
F[n-1]=0
A[0],A[n-1]=0,0
A[:,0],A[:,n-1]=0,0
A[0][0],A