Scribd Logo
Upload

Welcome to Scribd!

  • 21 views

    Exp 8 Simpsonrule

    Uploaded by

    dhruv
    The document describes a C program that uses Simpson's 1/3 rule to numerically integrate functions. The program takes in the lower and upper bounds of integration and the interval width, calculates the function values at sample points, and applies the Simpson's 1/3 rule formula to estimate the integral. The example function integrated is f(x) = 1/(1+x) from 0 to 1, and the program outputs the estimated integral value. Numerical integration allows calculating quantities such as areas, volumes, and fields that involve integrals.

    Copyright:

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

    Exp 8 Simpsonrule

    Uploaded by

    dhruv
    21 views3 pages
    The document describes a C program that uses Simpson's 1/3 rule to numerically integrate functions. The program takes in the lower and upper bounds of integration and the interval width, calculates the function values at sample points, and applies the Simpson's 1/3 rule formula to estimate the integral. The example function integrated is f(x) = 1/(1+x) from 0 to 1, and the program outputs the estimated integral value. Numerical integration allows calculating quantities such as areas, volumes, and fields that involve integrals.

    Original Title

    exp8simpsonrule

    Copyright

    © © All Rights Reserved

    Available Formats

    DOCX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document describes a C program that uses Simpson's 1/3 rule to numerically integrate functions. The program takes in the lower and upper bounds of integration and the interval width, calculates the function values at sample points, and applies the Simpson's 1/3 rule formula to estimate the integral. The example function integrated is f(x) = 1/(1+x) from 0 to 1, and the program outputs the estimated integral value. Numerical integration allows calculating quantities such as areas, volumes, and fields that involve integrals.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    21 views3 pages

    Exp 8 Simpsonrule

    Uploaded by

    dhruv
    The document describes a C program that uses Simpson's 1/3 rule to numerically integrate functions. The program takes in the lower and upper bounds of integration and the interval width, calculates the function values at sample points, and applies the Simpson's 1/3 rule formula to estimate the integral. The example function integrated is f(x) = 1/(1+x) from 0 to 1, and the program outputs the estimated integral value. Numerical integration allows calculating quantities such as areas, volumes, and fields that involve integrals.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    You are on page 1of 3

    Experiment 8

    Program for solving numerical integration by Simpson’s 1/3 rule

    Integration is an integral part in science and engineering to calculate things such as area,


    volume, total flux, electric field, magnetic field and many more. Here, we are going to take
    a look at numerical integration method  (Simpson’s 1/3 rule in particular using C
    language) to solve such complex integration problems.

    For this, let’s discuss the C program for Simpson 1/3 rule for easy and accurate
    calculation of numerical integration of any function which is defined in program.

    In the source code below, a function f(x) = 1/(1+x) has been defined. The calculation


    using Simpson 1/3 rule in C is based on the fact that the small portion between any two
    points is a parabola. The program follows the following steps for calculation of the
    integral.

     As the program gets executed, first of all it asks for the value of lower boundary
    value of x i.e. x0, upper boundary value of x i.e. xn and width of the strip, h.
     Then the program finds the value of number of strip as n=( x n – x0 )/h and checks
    whether it is even or odd. If the value of ‘n’ is odd, the program refines the value
    of ‘h’ so that the value of ‘n’ comes to be even.
     After that, this C program calculates value of f(x) i.e ‘y’ at different intermediate
    values of ‘x’ and displays values of all intermediate values of ‘y’.
     After the calculation of values of ‘c’, the program uses the following formula to
    calculate the value of integral in loop.
    Integral =  *((y0 + yn ) +4(y1 + y3 + ……….+ yn-1 ) + 2(y2 + y4 +……….+ yn-2 ))

     Finally, it prints the values of integral which is stored as ‘ans’ in the program.
    If f(x) represents the length, the value of integral will be area, and if f(x) is area, the output
    of Simpson 1/3 rule C program will be volume. Hence, numerical integration can be
    carried out using the program below; it is very easy to use, simple to understand, and
    gives reliable and accurate results.

    f(x) = 1/(1+x)

    #include<stdio.h>

    #include<conio.h>

    float f(float x)

    {
    return(1/(1+x));

    void main()

    int i,n;

    float x0,xn,h,y[20],so,se,ans,x[20];

    printf("\n Enter values of x0,xn,h: ");

    scanf("%f%f%f",&x0,&xn,&h);

    n=(xn-x0)/h;

    if(n%2==1)

    n=n+1;

    h=(xn-x0)/n;

    printf("\n Refined value of n and h are:%d %f\n",n,h);

    printf("\n Y values: \n");

    for(i=0; i<=n; i++)

    x[i]=x0+i*h;

    y[i]=f(x[i]);

    printf("\n %f\n",y[i]);

    }
    so=0;

    se=0;

    for(i=1; i<n; i++)

    if(i%2==1)

    so=so+y[i];

    else

    se=se+y[i];

    ans=h/3*(y[0]+y[n]+4*so+2*se);

    printf("\n Final integration is %f",ans);

    getch();

    You might also like