JAVA环境下带边界条件的线性常微分方程组（BVP）求解包咨询

以下是适配你需求的成熟Java实现选项：

• MTJ (Matrix Toolkit Java) 配合有限差分实现：MTJ是工业级的Java稀疏矩阵运算库，你所需要求解的方程n''(x) = -Dn(x) + n'(x)/x离散后会得到线性方程组，可直接调用MTJ的内置线性求解器计算结果，不需要自行实现底层矩阵运算，稳定性和计算效率都有保障。
• 基于Apache Commons Math 实现打靶法：虽然原生Apache Commons Math的ODE模块只支持初值问题求解，但你可以基于其已有的初值求解器封装打靶法逻辑，针对线性BVP收敛性非常稳定，开发量很小，适合快速实现需求。
• JFiniteElement 库：轻量级Java有限元实现库，原生支持一维常微分方程边值问题的离散与求解，你只需按文档代入方程系数和边界条件即可得到结果，不需要自行实现有限元离散逻辑。

注意：你给出的方程在x=0处存在奇点，若求解域包含x=0，需要单独推导该点的离散格式，避免数值计算错误。

内容的提问来源于stack exchange，提问作者bluesky