高斯型数值积分公式的Python程序

最新推荐文章于 2024-08-01 20:23:41 发布

Fo*(Bi)

最新推荐文章于 2024-08-01 20:23:41 发布

阅读量2.1k

点赞数

分类专栏：算法文章标签： python

本文链接：https://blog.csdn.net/weixin_48615832/article/details/111469075

版权

算法专栏收录该内容

69 篇文章 52 订阅

订阅专栏

在这里插入图片描述

from sympy import *

def f(t):
    f = 2000*log(140000/(140000-2100*t))-9.8*t
    return f

a = 8
b = 30
x = symbols('x')
truth = integrate(f(x),(x,a,b)).evalf()
print(truth) #真值

n = 2 #步长,就是将(a,b)区间分为多少个块
h = (b-a)/n
result = h* f(h*(-1/sqrt(3))+(b+a)/n) + h* f(h*(1/sqrt(3))+(b+a)/n)
print(result.evalf())