#include <math.h>

double get_solid_angle(double *pos, double *pmt, double *n, double r);