编写程序来计算不定积分可以通过以下步骤进行：
导入必要的数学库
对于数值计算，可以使用 `numpy` 库。
对于符号计算，可以使用 `sympy` 库。
定义被积函数
可以通过用户输入或者直接在代码中定义函数形式。
选择计算方法
数值方法：如梯形法则、辛普森法则、龙贝格法则等。
符号计算方法：如使用 `sympy` 库的 `integrate（）` 函数。
实现数值方法
例如，实现梯形法则的函数：
```python
def trapezoid_integral（f, a, b, n）:
h = （b - a） / n
integral_sum = （f（a） + f（b）） / 2
for i in range（1, n）:
x = a + i * h
integral_sum += f（x）
return integral_sum * h
```
实现符号计算方法
例如，使用 `sympy` 计算不定积分：
```python
import sympy as sp
def symbolic_integral（f, x）:
return sp.integrate（f, x）
```
调用函数并输出结果
根据用户输入或定义的参数，调用相应的函数并输出不定积分的近似值或符号结果。
示例
```python
import sympy as sp
定义变量
x = sp.symbols（'x'）
定义被积函数
f = 3*x2 + 2*x + 1
计算不定积分
indefinite_integral = symbolic_integral（f, x）
print（f"不定积分： {indefinite_integral} + C"）
```
建议
选择合适的工具：根据具体需求和精度要求选择数值方法或符号计算方法。
注意精度和效率：数值方法可能需要更多的计算资源，而符号方法可以提供精确结果但计算时间较长。
测试和验证：编写测试用例来验证程序的正确性和鲁棒性。
通过以上步骤和示例代码，你可以编写一个程序来计算不定积分。