我正在研究我的毕业项目,特别是关于流体动力学的项目,我有一个非线性方程组要解决,我选择牛顿方法,所以我必须通过matix的雅可比行列式(实际上是12x12矩阵).该矩阵中的每个元素都是在某个时刻评估的函数的导数,很难手动编写所有这些元素并计算每个导数; 系统看起来像:
f 1(x 1,x 2,x 3,...,x 12)= 0
其中x 1,x 2,x 3是变量(温度,压力......等)
我可以自动执行此操作吗？如果在Fortran中不可能,我可以使用其他脚本语言作为Python(sympy模块)吗？