我尝试在MATLAB中求解带有节点热源的四面体有限元的热扩散问题,这个节点取决于解矢量。非线性方程系统如下: 乙U” + A U = Q(T) 与B是热capactiy矩阵,A是导电性基体,q是源项和U作为温度。我使用了一个Adams-Bashforth /梯形规则预估校正器方案,并进行了Picard迭代,接着是时间步长控制。源项的温度在最后一个时间步的温度和预测器温度之间精确评估。这是预测校正器
我想弄清楚如何在Fortran中调用求解器来求解一些n维非线性方程。我遇到了具有一些非线性求解例程的SLATEC库,我尝试使用的库是SNSQE。 (旁白:也许你知道一些更好/更容易使用的解算器,请告诉!) 我的示例代码试图解决一个简单的二维问题。你可以手工解决它,但重点是找出解算器的工作原理。代码在这里: program test
implicit none
! declar