Scribd Logo
Upload

Welcome to Scribd!

  • 8 views

    Simpson Rules: 1 Composite Simpson's 1/3 Rule

    Uploaded by

    Omar Ardila
    The document discusses Simpson's 1/3 rule for numerical integration and provides an example calculation using Python code to estimate the integral of a function from 0 to 0.8 using 4 intervals.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Simpson Rules: 1 Composite Simpson's 1/3 Rule

    Uploaded by

    Omar Ardila
    8 views1 page
    The document discusses Simpson's 1/3 rule for numerical integration and provides an example calculation using Python code to estimate the integral of a function from 0 to 0.8 using 4 intervals.

    Original Description:

    simpson1

    Original Title

    Simpson Rules

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document discusses Simpson's 1/3 rule for numerical integration and provides an example calculation using Python code to estimate the integral of a function from 0 to 0.8 using 4 intervals.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    8 views1 page

    Simpson Rules: 1 Composite Simpson's 1/3 Rule

    Uploaded by

    Omar Ardila
    The document discusses Simpson's 1/3 rule for numerical integration and provides an example calculation using Python code to estimate the integral of a function from 0 to 0.8 using 4 intervals.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 1

    Simpson rules

    June 10, 2022

    1 Composite Simpson’s 1/3 rule


    a=0 ; b = 0.8 ; n=4

    f (x) = 0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5

    ∫ 0.8
    (0.2 + 25x − 200x2 + 675x3 − 900x4 + 400x5 )dx = 1.640533
    0

    ∑
    n−1 ∑
    n−2
    f (x0 ) + 4 f (xi ) + 2 f (xj ) + f (xn )
    i=1,3,5 j=2,4,6
    I = (b − a) ∗
    3n

    [1]: import numpy as np

    a=0; b=0.8; n=4;


    x=np.linspace(a,b,n+1)
    y=0.2 + 25*x - 200*x**2 + 675*x**3 - 900*x**4 + 400*x**5

    s1=[]; sum1=0
    for i in range(1,n,2):
    s1.append(y[i])
    sum1=np.sum(s1)

    s2=[]; sum2=0
    for j in range(2,n-1,2):
    s2.append(y[j])
    sum2=np.sum(s2)

    I=(b-a)*(y[0]+4*sum1+2*sum2+y[-1])/(3*n)

    print("La integral estimada es: ", I)

    La integral estimada es: 1.6234666666666717

    You might also like