/*pow1 evaluate luminosity of continuum */
double pow1(double freq);