/* integro2.c */
#include <stdlib.h>
#include <stdio.h>
#include <math.h>

typedef double (*funktio_tyyppi)(double x);

double oma_funktio(double x)
{
  return 2*x-5;
}

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

  dx = (x2-x1)/tiheys;

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

  return summa;
}

int main(void)
{
  printf("Integraali sin(x) v„lilt„ [0,pi] on noin %7.5lf\n",
          integroi(sin,0,M_PI,100));
  printf("Integraali 2x-5 v„lilt„ [0,0.5] on noin %7.5lf\n",
          integroi(oma_funktio,0,0.5,100));
  printf("Integraali e^x v„lilt„ [0,1] on noin %7.5lf\n",
          integroi(exp,0,1,500));
  return 0;
}