/**********************************************************************
 *
 * File: riemann.c (10-Mar-2006)
 *                 (12-Mar-2007)
 *                 (02-Mar-2008)
 *
 * Riemannsumme als Approxmation des Integrals der Funktion
 *
 *   f(x) = 4/(1+x*x)
 *
 * a,b...Intervallgrenzen
 * n.....Anzahl der Teilintervalle
 * s.....Wert des Integrals
 *
 **********************************************************************/

#include <stdio.h>

/**********************************************************************/

float f(float x) {

/**********************************************************************/

  float y;

  y=4.0/(1.0+x*x);

  return y;

}

/**********************************************************************/

int main() {

/**********************************************************************/

  int i,n;
  float a,b,dx,s,x;

  printf(" a,b,n=");
  scanf("%f,%f,%d",&a,&b,&n);

  dx=(b-a)/n;

  s=0.0;
  for(i=0;i<=n-1;i++) {
    x=a+i*dx;
    s=s+f(x);
  }
  s=dx*s;

  printf(" integral=%f\n",s);

  return 0;

}

/**********************************************************************/