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

题目：碳原子在液态金属下的扩散模型 `C] t2^  
qBKIl= ne  
  t[iE >  
,}`II|.oB  
目的：通过研究碳原子在液态金属下的扩散模型,可以根据需要对金属进行渗碳，从而改善金属的硬度，以适合不同的需要。 #/5jW
H7U  
hofZpM  
  x*bM C&Ea  
In^$+l%O[  
原理：用偏微分方程来模拟碳原子的扩散过程，本次研究的是二维扩散，所用方程和边界条件如下： yy/wSk  
oDM}h +  
3x'BMAA+  
Ojie.+'SB  
).-B@&Eu%  
其中u表示碳的质量分数。对此扩散方程采用交替方向隐格式（ADI）方法进行求解。对于多维差分格式，采用交替方向隐格式可以减少计算量和复杂度。而且该方法的误差阶和稳定性都很好。 JBi<TDm/  
H.w
p{m{  
  GaM#a[p  
'iUfr@  
方法：在将偏微分方程化成差分方程的过程中采用P-R格式，P-R格式是在第K层到第K+1层中引入一个过渡层，即第K+1/2层。从第K层到第K+1/2层时，对 在过渡层，对 在第K层上用二阶中心差商逼近，然后从第K+1/2层到第K+1层，对 仍在过渡层，对 在第K+1层上用二阶中心差商逼近。然后用追赶法解具有三对角阵的方程组。 4wBMBCJ;P  
K~v"%sG{`  
编程框图：有限差分求解的程序框图如下图所示。 X^fMt]  
>8>`-  
S37Bl5W  
g3Kc? wTC  
  ~~:i+-[  
y=5s~7]  
  7(Y!w8q&^  
&``;1/J*W  
  u -CY-  
s{`r$:!  
实现步骤：1、偏微分方程化为差分方程 duT2:~H2  
nmS3  
          2、对差分方程用追赶法求解 0(:"q!h  
3EF|1B/5  
          3[/font .. V0ulIKck  
Mq6.!j  

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

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

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

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