(fitteia's extended models library)
double GLmix(double f, double a, double b, double f0, double m)
{
double GAUSS1(),LORTZN();
return m*GAUSS1(f,a,b,f0)+(1-m)*LORTZN(f,a,b,f0);
}