double GAULOR(double f, double a, double b, double f0, double type)

{

double GAUSS1(),LORTZN();

if(type>0.5){ return LORTZN(f,a,b,f0); }

else { return GAUSS1(f,a,b,f0); }

}