编写一个计算二阶常系数齐次线性微分方程

要编写一个计算二阶常系数微分方程的振荡解的程序，可以按照以下步骤：

1. 引入所需的库：例如numpy、matplotlib等。

2. 定义微分方程：在程序中，需要定义微分方程，可以使用lambda函数来定义方程。

3. 设定初值条件：需要设定初值条件，例如初始位置和初始速度。

4. 定义求解函数：可以使用odeint函数对微分方程进行求解，并定义一个函数来调用odeint函数。

5. 绘制结果：将求解得到的结果绘制出来，可以使用matplotlib库的plot函数来绘制曲线。

以下是一个可能的示例代码：

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint

# 定义微分方程
f = lambda y, t: [y[1], -y[0]]

# 设定初值条件
y0 = [1, 0]
t = np.linspace(0, 10, 101)

# 定义求解函数
def solve_ode():
    sol = odeint(f, y0, t)
    return sol[:, 0]

# 绘制结果
plt.plot(t, solve_ode())
plt.xlabel('t')
plt.ylabel('y')
plt.title('Solution of the Harmonic Oscillator')
plt.show()

在这个例子中，定义的微分方程为二阶常系数微分方程y'' + y = 0，其对应的振荡解为y = sin(t)。程序求解得到的结果也应该是这个振荡解的曲线。