#include <math.h>
#include <algorithm>
#include <iostream>
using namespace std;

int num=0;
double sigma, eksp, delta, start, stop,curx,nextx,areal;
double part1,part2,x,power;

double e = 2.7182818284590452353602874713526624977572;
double pi = 3.1415926535897932384626433832795028841971693993751058209749445923078;



int main(){
    while(true){
        areal = 0;
        cout << "Forventenings verdi: ";
        cin >> eksp;
        cout << "Sigma: ";
        cin >> sigma;
        cout << "Start: ";
        cin >> start;
        cout << "Slutt: ";
        cin >> stop;

        delta = (stop-start)/1000000;
        x = start;
        part1 = 1/(sigma*sqrt(2*pi));


        while(x+delta < stop){

            power = -1*(pow((x-eksp),2)/(2*pow(sigma,2)));
            part2 = pow(e,power);
            curx = part1*part2;

            power = -pow((x+delta-eksp),2)/(2*pow(sigma,2));
            part2 = pow(e,power);
            nextx = part1*part2;

            if(curx > nextx){
                areal += delta*nextx; //square
                areal += (delta*(curx-nextx))/2; // Triangle


            }
            else{

                areal += delta*curx; //Square
                areal += (delta*(nextx-curx))/2; //Triangle

            }

            x += delta;
            if(areal > 1.0){
                break;
            }
        }
        cout << "X: " << x << "   ";
        cout << "Area: " << areal << endl << endl;

    }
    return(0);
}