北太天元求解常微分方程

以下是使用北太天元求解常微分方程的示例代码：

import baltamatica as bt

# 定义微分方程
def f(t, y):
    return y - t**2 + 1

# 定义初始条件
y0 = 0.5

# 定义求解区间
t_span = [0, 2]

# 使用solve_ivp函数求解微分方程
sol = bt.solve_ivp(f, t_span, y0)

# 输出结果
print(sol.t)
print(sol.y[0])

在上面的代码中，我们首先导入了北太天元库，并定义了一个微分方程f(t, y)，然后定义了初始条件y0和求解区间t_span。接着，我们使用solve_ivp函数求解微分方程，并将结果存储在sol变量中。最后，我们输出了求解结果。