diff --git a/按方法整理/常微分方程-经典龙格-库塔格式.py b/按方法整理/常微分方程-经典龙格-库塔格式.py index fa3deef..f04207b 100644 --- a/按方法整理/常微分方程-经典龙格-库塔格式.py +++ b/按方法整理/常微分方程-经典龙格-库塔格式.py @@ -19,7 +19,7 @@ if __name__=="__main__": ########################################################################## x0 = 0 # x的初始值(左边界)换成题干里面的################# y0 = -1 # y的初始值换成题干里面的################# - fxy = lambda x,y: x + y #f(x,y)换成题干里面的################# + fxy = lambda x,y: x+y #f(x,y)换成题干里面的################# real_fx = lambda x: -x-1 #真实函数换成题干里面的################# h = 0.1 #步长换成题干里面的################# xk = 2 #x的右边界换成题干里面的################# diff --git a/按方法整理/常微分方程-阿达姆斯预测-校正方法.py b/按方法整理/常微分方程-阿达姆斯预测-校正方法.py index d3dc496..85040a4 100644 --- a/按方法整理/常微分方程-阿达姆斯预测-校正方法.py +++ b/按方法整理/常微分方程-阿达姆斯预测-校正方法.py @@ -49,7 +49,7 @@ def AdamusFix(k,first_ys,x0,y0,h,xk,fxy,fx_real): if __name__ == "__main__": # 定义初始值和参数##################################################################################### - k = 3 # 阿达姆斯k+1步显式方法,精度为k+1阶,最常用k=3,其他阶数我没试过 P253 + k = 3 # 阿达姆斯k+1步校正方法,精度为k+1阶,最常用k=3,其他阶数我没试过 P253 first_ys = [1, 1.0954, 1.1832, 1.2649] # 前几个y的值,可用龙格-库塔计算或者知道精确解自己算出来 x0 = 0.0 # 初始x值 y0 = 1.0 # 初始y值