package piirto;

/**
 * Ohjelmalla integroidaan numeerisesti funktio-olio.
 * @author Vesa Lappalainen
 * @version 1.0, 25.03.2003
 */
public class Integroi2 {

  interface FunktioRR {
    public double f(double x);
  }

  public static class SinFun implements FunktioRR {
    public double f(double x) { return Math.sin(x); }
  }

  public static class ExpFun implements FunktioRR {
    public double f(double x) { return Math.exp(x); }
  }

  public static class OmaFun implements FunktioRR  {
    public double f(double x) { return 2*x- 5;      }
  }

  public static class Poly2 implements FunktioRR  {
      double a,b,c;
      Poly2(double a, double b, double c) {
          this.a = a;
          this.b = b;
          this.c = c;
      }
      public double f(double x) { return a*x*x + b*x + c;   }
    }

  public static double integroi(FunktioRR 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.f(x)*dx;

    return summa;
  }


  public static void main(String[] args)  {
    double ifx;
    ifx = integroi(new OmaFun(),0,5,100);
    System.out.printf("Integraali omafun(x) väliltä [0,5] on noin %7.5f%n",ifx);

    FunktioRR sinOlio = new SinFun();
    ifx = integroi(sinOlio,0,Math.PI,100);
    
    double y = sinOlio.f(1);  ///  double y = Math.sin(1);
    
    System.out.printf("Integraali sin(x) väliltä [0,pi] on noin %7.5f%n",ifx);
    ifx = integroi(new ExpFun(),0,1,500);
    System.out.printf("Integraali exp(x) väliltä [0,1] on noin %7.5f%n",ifx);
    
    FunktioRR p2 = new Poly2(1,2,2);
    FunktioRR p2b = new Poly2(4,-2,2);
    
    ifx = integroi(p2,0,1,500);
    System.out.printf("Integraali P2(x) väliltä [0,1] on noin %7.5f%n",ifx);
  }
}