#define  A 10.0  // Integral starting position
#define  B 22.0  // Integral ending position
#define  W   .1  // Integral step size

#include <stdio.h>
#include <mpi.h>
#define  round(v) (float)(int)(v / W) * W

double f(double x) {
    return x * x;
}

int main(int argc, char* argv[]) {
    int i;
    int id, p;
    MPI_Init(&argc, &argv);
    MPI_Comm_size(MPI_COMM_WORLD, &p);
    MPI_Comm_rank(MPI_COMM_WORLD, &id);

    // Calculate job boundaries and print them in order
    double a, b, da, db, wpp;
    wpp   = (B - A) / p;
    da    = wpp * id;
    db    = wpp * (id + 1);
    a     = A + round(da);
    b     = A + round(db);
    for (i = 0; i < p; i++) {
        if (i == id)
            printf("Process #%d: f(%f) to f(%f)\n", id, a, b);
        MPI_Barrier(MPI_COMM_WORLD);
    }

    MPI_Finalize();
}