[求助]用差分法计算抛物方程的时候,这个边界条件怎么带进去啊?

[fengwanyao]

作者: amethystcat

方程:
du/dt=a*du^2/dz^2 a是常数,标准的抛物方程
边值条件: du(0,t)=b*(du/dz)|(z=0) b是常数,(du/dz)|z=0表示z=0时du/dz的值
初值条件:u(z,0)=c*f(z) c是常数,f是z的函数

编程的时候,这个边值条件怎么处理啊?
求教
有对这个比较熟的大侠留个联系方式吧,想详细请教
多谢了!~

边值条件使用中心差分离散，并要求方程在边界上成立，du^2/dz^2也采用中心差分离散，联立消去边界外的点，即可