// intergroi2.cpp
#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=x1+dx/2 ; x<x2; x+=dx)
    summa += f(x)*dx;

  return summa;
}

int main(void)
{
  double ifx;
  ifx = integroi(oma_funktio,0,5,100);
  printf("Integraali väliltä [0,5] on noin %7.5lf\n",ifx);

  ifx = integroi(sin,0,M_PI,1000);
  printf("Integraali sin(x) väliltä [0,pi] on noin %7.5lf\n",ifx);

  ifx = integroi(exp,0,1,500);
  printf("Integraali exp(x) väliltä [0,1] on noin %7.5lf\n",ifx);

  return 0;
}