复化求积公式计算积分

import numpy as np


def myfun(x):
    if x == 0:
        return 1
    else:
        return np.sin(x)/x


def tixing(fun, a, b, h):
    n = int((b-a)/h)+1
    x = np.linspace(a, b, n)
    sum1 = 0
    for i in range(1, n-1):
        sum1 = sum1 + h*fun(x[i])
    sum1 = sum1 + h/2*(fun(x[0])+fun(x[-1]))
    return sum1


def simpson(fun, a, b, h):
    n = int((b-a)/h)+1
    x = np.linspace(a, b, n)
    sum1 = 0
    for i in range(1, n-1):
        sum1 = sum1 + (h/6)*fun(x[i])*2+4*(h/6)*fun((x[i-1]+x[i])/2)
    sum1 = sum1 + h/6*(fun(x[0])+fun(x[-1]))+(2/3)*h*fun((x[i-1]+x[i])/2)
    return sum1


s1 = tixing(myfun, 0, 1, 0.125)
print('复化梯形求积公式积分结果:')
print(s1)
print('复化辛普森求积公式积分结果:')
s2 = simpson(myfun, 0, 1, 0.125)
print(s2)

在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

leetteel

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值