#include <math.h>

double dist(int dim, double* p1, double* p2){

        int i;
        double d = 0.0;

        for(i=0; i < dim; i++){
                d += (p1[i] - p2[i])*(p1[i] - p2[i]);
        }

        return sqrt(d);
}