Scribd
Upload
47 views

Unit5.ipynb - Numerical Integration

Python Program using Trapezoidal rule Simpsons 1/3 rd rule and simpsons 3/8 th rule to solve definite integration problems

Uploaded by

BhikulalJaiswal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
47 views

Unit5.ipynb - Numerical Integration

Python Program using Trapezoidal rule Simpsons 1/3 rd rule and simpsons 3/8 th rule to solve definite integration problems

Uploaded by

BhikulalJaiswal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Double-click (or enter) to edit

Trapezoidal Rule

# Trapezoidal Method

def f(x):               # Define function to integrate
    return 1/(1 + x**2)
x0,xn,n=5,10,900
# Implementing trapezoidal method
h = (xn - x0) / n  # calculating step size
# Finding sum 
sum = f(x0) + f(xn)   
for i in range(1,n):
        k = x0 + i*h
        sum = sum + 2 * f(k)    
# Finding final integration value
integration = sum * h/2    
print("Integration result by Trapezoidal method is: %0.6f" % (integration) )    

Integration result by Trapezoidal method is: 0.097727

Simpsons 1/3 rd Rule

#Simpson's 1/3 Rule
def f(x):               # Define function to integrate
    return 1/(1 + x**2)
x0,xn,n=5,10,999
# Implementing Simpson's 1/3 
h = (xn - x0) / n  # calculating step size
# Finding sum 
sum = f(x0) + f(xn)   
for i in range(1,n):
    k = x0 + i*h
        
    if i%2 == 0:
            sum = sum + 2 * f(k)
    else:
            sum = sum + 4 * f(k)
   
 # Finding final integration value       
integration = sum * h/3    
print("Integration result by Simpsons 1/3 rd Rule is: %0.6f" % (integration) )    

I t ti lt b T id l th d i 0 097710
Simpsons 3/8 th rule

#Simpson's 3/8th Rule
def f(x):               # Define function to integrate
    return 1/(1 + x**2)
x0,xn,n=5,10,900
# Implementing Simpson's 3/8
h = (xn - x0) / n  # calculating step size
# Finding sum 
sum = f(x0) + f(xn)   
for i in range(1,n):
    k = x0 + i*h
        
    if i%3 == 0:
            sum = sum + 2 * f(k)
    else:
            sum = sum + 3 * f(k)
   
 # Finding final integration value       
integration = sum * 3*h/8    
print("Integration result by Simpsons 3/8 th Rule is: %0.6f" % (integration) )    

import sympy as sp
a=0
b=5
x = sp.Symbol('x')
f=.5*x+3
f.integrate((x,a,b))

21.25

import numpy as np
f=lambda x :400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
a=0
b=0.6
x=[a,b]
y=[f(a),f(b)]
np.trapz(y,x)

0.8591999999999943

import sympy as sp
a=0
b=0.6
x = sp.Symbol('x')
f=400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
f.integrate((x,a,b))

0.963600000000001

import sympy as sp
x = sp.Symbol('x')
f=400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
f.integrate((x,a,b))

576613.583333333

Double-click (or enter) to edit

import matplotlib.pyplot as plt
a=0
b=0.6
f=lambda x :400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
x =np.arange(a,b,0.01)
y=f(x)
plt.plot(x,y)
plt.fill_between([a,b],[f(a),f(b)],alpha = 0.5)
plt.show()

import matplotlib.pyplot as plt
import numpy as np
a=0
b=0.6
f=lambda x :400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
x =np.arange(a,b,0.01)
y=f(x)
plt.plot(x,y)
def SingleTrap(f,a,b):
  plt.fill_between([a,b],[f(a),f(b)],alpha=0.5)
  I=(b-a)*(f(a)+f(b))/2
  return I
I= SingleTrap(f,a,b)
print(I)
plt.show()
0.8591999999999943

import matplotlib.pyplot as plt
import numpy as np
def SingleTrap(f,a,b):
  plt.fill_between([a,b],[f(a),f(b)],alpha=0.5)
  I=(b-a)*(f(a)+f(b))/2
  return I
a,b,n=0,0.6,10
a0=a
h=(b-a)/n
I_sum = 0
f=lambda x :400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
x =np.arange(a,b,0.01)
y=f(x)
plt.plot(x,y)
for i in range(n):
  b0=a0+h
  I= SingleTrap(f,a0,b0)
  I_sum = I_sum + I
  a0 =b0
print(I_sum)
plt.show()

0.9548577599999999
import matplotlib.pyplot as plt
import numpy as np
def SingleSimpson13(f,a,b):
  h=(b-a)/2
  plt.fill_between([a,b],[f(a),f(b)],alpha=0.5)
  x0=a
  x1=x0+h
  x2=x1+h
  I=(b-a)*(f(x0)+4*f(x1)+f(x2))/6
  return I
a,b,n=0,0.6,10
a0=0
I_sum = 0
f=lambda x :400*x**5 -900*x**4+675*x**3 -200*x**2+25*x-0.2
x =np.arange(a,b,0.01)
y=f(x)
plt.plot(x,y)
for i in range(n):
  b0=a0+h
  I= SingleSimpson13(f,a0,b0)
  I_sum = I_sum + I
  a0 =b0
print(I_sum)
plt.show()

0.9635805599999998

You might also like