追赶法\Thomas算法 在常微分方程两点边值问题的差分解法，紧差分解法中的应用（附python代码实现及作图）

偏微分数值解法课程作业系列

Chapter 1 常微分方程两点边值问题的差分解法

题目引入：

也可为
注：差分法紧差分法最终实现形式都为追赶法求解 代码通用

分析：

6题分析

第7题紧差分法分析

下附上第七题紧差分法代码实现：

import matplotlib.pyplot as plt
import numpy as np

def thomas(N,a,b,c,d):
    l = np.zeros((N+1,1))
    u = np.zeros((N+1,1))
    x = np.zeros((N+1,1))
    y = np.zeros((N+1,1))
    u[0] = b[0]
    for i in range(N):
        l[i+1] = a[i+1]/u[i