改进欧拉法

最新推荐文章于 2024-06-19 10:23:56 发布

小盼你最萌哒

最新推荐文章于 2024-06-19 10:23:56 发布

阅读量4.6k

点赞数 2

分类专栏：数值分析

本文链接：https://blog.csdn.net/qq_27508477/article/details/84865011

版权

数值分析专栏收录该内容

9 篇文章 2 订阅

订阅专栏

#P236 例2
import numpy as np

x=0
y=0
yn=1
h=0.1
accuracy=6
for i in range(0,10):
    
    yp=round(y+h*(x-y),accuracy)
    yq=round(y+h*(x-yp+h),accuracy)
    #yn=round(x*(3+x**2)/3/(1+x**2),accuracy)
    yn=round(x-1+np.e**(-x),accuracy)
    ans=round(yn-y,accuracy)
    print(i,x,y,yn,ans,sep='\t\t')
    y=round((yp+yq)/2,accuracy)
    x=round(x+h,accuracy)

#P264 习题九 3.（2）
import numpy as np

def Fx(x,y):
    return (2*x+1)**0.5

def fx(x,y):
    return y-2*x/y
x=0
y=1
h=0.1
accuracy=4
for i in range(0,11):
    yp=round(y+h*fx(x,y),accuracy)
    yq=round(y+h*fx(x+h,yp),accuracy)
    #yn=round(x*(3+x**2)/3/(1+x**2),accuracy)
    yn=round(Fx(x,y),accuracy)
    ans=round(yn-y,accuracy)
    print(i,x,y,yn,ans,sep='\t\t')
    y=round((yp+yq)/2,accuracy)
    x=round(x+h,accuracy)