/*ffun1 derive flux at a specific energy, for one continuum */
double ffun1(double xnu);