里查得森(Richardson)法解偏微分方程的程序。解的是定步長的拋物型偏微分方程:
du/dx - a * d2u/dx2 = 0
在程序中可以更改 a 的值以實現(xiàn)不同系數(shù)的解。由于它是一種顯式格式,解的結(jié)果不穩(wěn)定
在Matlab7下面運行發(fā)現(xiàn)竟然有一些值達到十的幾十次方,這是算法本身決定的。不是程序有錯誤
標(biāo)簽:
Richardson
dx
偏微分方程
程序
上傳時間:
2014-12-03
上傳用戶:1159797854