# -*- coding: utf-8 -*- """ Created on Sun Dec 24 15:29:08 2017@author: www """ #本程序是用四階龍格庫塔法求解課本（數(shù)值計(jì)算方法 馬東升）P242頁的例7-3 #fun為指定的導(dǎo)數(shù)的函數(shù) #rf4為四階龍格庫塔法def fun(x,y):f = y - (2 * x / y)return f#input # x0,y0:初始給出的x0值，y0值 # h ：步長 # N ：迭代次數(shù) #print # x1,y1:每次迭代輸出的結(jié)果 def rf4(x0,y0,h,N):n = 1while(n != N):x1 = x0 + hk1 = fun(x0, y0)k2 = fun(x0+h/2, y0+h*k1/2)k3 = fun(x0+h/2, y0+h*k2/2)k4 = fun(x1, y0+h*k3)y1 = y0 + h * (k1 + 2* k2 + 2 * k3 + k4) / 6print("%.2f, %.6f" %(x1, y1))n = n + 1x0 = x1y0 = y1def main():rf4(0,1,0.2,5)main()

?

總結(jié)

以上是生活随笔為你收集整理的python实现四阶龙格库塔法的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。