結果 : implicit euler method code