帮忙做一下这个问题  多谢了

题目：碳原子在液态金属下的扩散模型 uCUBs(iD  
[;X YT  
  }1$8)zH  
fA+M/}=  
目的：通过研究碳原子在液态金属下的扩散模型,可以根据需要对金属进行渗碳，从而改善金属的硬度，以适合不同的需要。 s&fU|Jk8  
c85B-/  
  r&TxRsg{  
c%B=TAs5c  
原理：用偏微分方程来模拟碳原子的扩散过程，本次研究的是二维扩散，所用方程和边界条件如下： O050Q5zy  
B8IfE`  
AJmS1 B  
(@&|  
wvq<5gy}  
其中u表示碳的质量分数。对此扩散方程采用交替方向隐格式（ADI）方法进行求解。对于多维差分格式，采用交替方向隐格式可以减少计算量和复杂度。而且该方法的误差阶和稳定性都很好。 M)b`~|Wt  
se)I
2T{J  
  [KjQW/sb'  
QZwRg&d<o  
方法：在将偏微分方程化成差分方程的过程中采用P-R格式，P-R格式是在第K层到第K+1层中引入一个过渡层，即第K+1/2层。从第K层到第K+1/2层时，对 在过渡层，对 在第K层上用二阶中心差商逼近，然后从第K+1/2层到第K+1层，对 仍在过渡层，对 在第K+1层上用二阶中心差商逼近。然后用追赶法解具有三对角阵的方程组。 _IY)<'d  
}f
({03$  
编程框图：有限差分求解的程序框图如下图所示。 ]=_BK!O  
.b]oB_  
$~`(!pa:  
-7+Fb^"L  
  !p9B
H6$`  
sH&8"5BT%  
  DPe`C%Oc1  
bUqO.FZ[  
  a%-Yl%#  
ZW\h,8%  
实现步骤：1、偏微分方程化为差分方程 bxyU[`  
Pb~S{):  
          2、对差分方程用追赶法求解 x3WY26e  
?2G^6>O`  
          3[/font .. 1hMk\ -3S  
'.^JN
@  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

这个……不会是要人帮你做程序吧？如果是免费的，你还是用商业软件算吧；要程序的话，嘿嘿，不好意思，那可是要收费的！ iI_Fbw8  
     看了下的的附件，那个方程是抛物型的，求解上应该没有太大困难（要是想偷懒，Visual Fortran、Intel Fortran的IMSL库都有专门的函数直接求解，根本不需要太复杂的编程：做个主程序，定义好数组、边界条件，调用相应函数就可以了的。

非常感谢啊！我也是刚开始学，所以不太懂。

你查查IMSL库手册吧，里边有解这种方程的函数