public class ExponentialCDF2 {


    public static void main (String[] argv)
    {
        Function F = makeExponentialCDF ();
	double a = 0, b = 3;          // Limit this to the interval [0,3]
	int M = 20;                   // Number of intervals.
	double delta = (b-a) / M;     // Interval size.

        double exValue = 0;

        // INSERT YOUR CODE HERE.

        System.out.println ("Exponential ex: " + ex);
    }

    static Function makeExponentialCDF ()
    {
	double a = 0, b = 3;
	int M = 50;                   // Number of intervals.
	double delta = (b-a) / M;     // Interval size.

	double[] intervalCounts = new double [M];
	double numTrials = 1000000;

	for (int t=0; t<numTrials; t++) {
            // Random sample:
	    double y = RandTool.exponential (2);
            // Find the right interval:
            int k = (int) Math.floor ((y-a) / delta);
            // Increment the count for every interval above and including k.
            for (int i=k; i<M; i++) {
                intervalCounts[i] ++;
            }
	}

	// Now compute probabilities for each interval.
	double[] cdf = new double [M];
	for (int k=0; k<M; k++) {
	    cdf[k] = intervalCounts[k] / numTrials;
	}

        // Build the CDF. Use mid-point of each interval.
        Function F = new Function ("Uniform cdf");
	for (int k=0; k<M; k++) {
	    double midPoint = a + k * delta + delta/2;
            F.add (midPoint, cdf[k]);
	}

        return F;
    }

}