Wednesday - July 17, 2024
#Fortran #Algorithm #Blog #C Program

Fortran ~ Simpson’s 1/3rd Rule

Fortran ~ Simpson’s 1 3rd Rule

  • By putting n=2 into a basic equation by integration and taking the curve through (X0, Y0), (X1, Y1) and (X2, Y2) as a polynomial of degree 2 so that difference of order Three and greater than three become zero than we get.
  • F(X)=3.12*X**3+1.54*X**2-6.23
  • This is known as Simpson’s 1/3rd Rule.
simpsons 13rd rule example, simpsons 13 rule ppt, simpsons 13 rule in c, simpsons 13 rule formula derivation, simpsons 13 rule calculator, simpsons 13rd rule and direct integration gives same results if, composite simpson's rule example, limitations of simpson's rule, simpson's rule matlab, simpson's rule formula for volume, solved examples of numerical integration, simpson's 3/8 rule matlab, simpson 1/3 rule algorithm,

Remarks Simpson’s 1/3rd Rule

  • Simpson’s 1/3rd rule requires the division of the interval [X0, Xn] into an even number of sub-intervals of width h.In this rule, the interpolating polynomial is of degree 2.
  • Therefore this rule is also known as a parabolic rule. 
PROGRAM SIMPSON
C
C         PROGRAM TO DEMONSTRATE APPLICATION OF SIMPSONS 1/3rd RULE
C
C        THIS PROGRAM NUMERICALLY INTEGRATES GIVEN EXPRESSION WITHIN
C        SPECIFIED LIMITS OF INTEGRATION
C
           F(X)=3.12*X**3+1.54*X**2-6.23
C
           WRITE (*,*)'LIMITS OF INTEGRATION (A,B)'
           READ(*,*)A,B
           WRITE(*,*)'NUMBER OF STRIPES(Number must be Even)'
           READ(*,*)NS
           IF(NS/2*2.NE.NS)THEN
           WRITE(*,*)'NUMBER OF STRIPES (must be Even)'
           WRITE(*,*)'PROGRAM TERMINATED...'
           STOP
            ENDIF
C         STRIP SIZE
            H=(B-A)/NS
            AREA=F(A)+F(B)
            DO 100 X=A+H,B-H,2*H
                  AREA=AREA+4.0*F(X)
100      CONTINUE
            DO 200 X=A+2*H,B-2*H,2*H
                  AREA=AREA+2.0*F(X)
200      CONTINUE
            AREA=AREA*H/3.0
            WRITE(*,*)'AREA UNDER THE CURVE=',AREA
            STOP
            END

#OUTPUT

LIMIT OF INTEGRATION (A, B)
1 2
NUMBER OF STRIPS (Number must be Even)
2
AREA UNDER THE CURVE = 54.546665
Fortran ~ Simpson’s 1/3rd Rule

Fortran ~ Trapezoidal rule