計算むずかしい

円柱系の表面拡散による濃度プロファイルの変化を計算していました。
3次の項までで考慮したんですけど、拡散時間が長くなると中心付近での計算結果に異常が(´・ω・｀)

何次まで考慮すればいいんだよ(´；ω；｀)