simpson

//Cubatura

#include <stdio.h>
#include <conio.h>
#include <math.h>
float I_Simpson(float,float,float,float,int,int);
float f(float,float);
int main(void)
{
    float a,b,c,d;
    int m,n;

    printf("\nIntroduceti a,b,c,d : ");
    scanf("%f,%f,%f,%f",&a,&b,&c,&d);
    printf("\nIntroduceti m,n : ");
    scanf("%d,%d",&m,&n);
    printf("\nRezultatul este = %f",I_Simpson(a,b,c,d,m,n));
    getch();

    return 1;   // din cauza functiei main()
}
float I_Simpson(float a,float b,float c,float d,int m,int n)
{
   float h,k;
   float sum;
   int i,j;
   h=(b-a)/n;
   k=(d-c)/m;
   sum=0;
   i=1;
   do {
      j=1;
      do {
    sum=sum+h*k*
        (
            f(a+(i-1)*h,c+(j-1)*k)+
            f(a+(i+1)*h,c+(j-1)*k)+
            f(a+(i-1)*h,c+(j+1)*k)+
            f(a+(i+1)*h,c+(j+1)*k)+
            4*f(a+i*h,c+(j+1)*k)+
            4*f(a+i*h,c+(j-1)*k)+
            4*f(a+(i-1)*h,c+j*k)+
            4*f(a+(i+1)*h,c+j*k)+
            16*f(a+i*h,c+j*k)
        )/9;
    j=j+2;
    } while(j<=m-1);
    i=i+2;
    } while(i<=n-1);
    return sum;
}
float f(float x,float y)
{
        return x+y ;
    return 1/(x*y);
}