// integro2.c
// Integroidaan f(x) valilta [0,pi]
// Käytetään funktio-osoitinta
// vl -92
#include <stdlib.h>
#include <stdio.h>
#include <math.h>


typedef double (*tFunktio)(double x);

double integroi(tFunktio f,double x1, double x2, int tiheys)
{
  double x,dx,summa=0;

  dx = (x2- x1)/tiheys;

  for (x=x1+dx/2 ; x<x2; x+=dx)
    summa += f(x)*dx;

  return summa;
}

int main(void)
{
  printf("Integraali sin(x) valilta [0,pi] on noin %7.5lf\n",
          integroi(sin,0,M_PI,100));
  return 0;
}