python改进欧拉法

最新推荐文章于 2022-10-19 22:38:00 发布

厄巧克力

最新推荐文章于 2022-10-19 22:38:00 发布

阅读量413

点赞数

分类专栏：计算方法文章标签： python 开发语言

本文链接：https://blog.csdn.net/weixin_49530111/article/details/126284885

版权

计算方法专栏收录该内容

11 篇文章 2 订阅

订阅专栏

import sympy
x=sympy.symbols("x")
y=sympy.symbols("y")
fy=input("f(x,y)=")
a=float(input("a="))
b=float(input("b="))
n=int(input("n="))
α=float(input("α="))
fy=sympy.simplify(fy)
X=[]
Y=[]
h=(b-a)/n
X.append(a)
Y.append(α)
for i in range(n):
    x1=X[i]+h
    X.append(x1)
    yp=Y[i]+h*(fy.evalf(subs={x:X[i],y:Y[i]}))   #对x,y赋值
    yc=Y[i]+h*(fy.evalf(subs={x:X[i+1],y:yp}))
    y1=(1/2)*(yp+yc)
    Y.append(y1)
print('x:',X)
print('y:',Y)