double f1(){
    double denominator = g/sqr(v0)-4.*H/sqr(L);
    if(denominator <= 0)
        return -1;
    double x = sqrt(2.*p)/denominator;
    if(0 <= x && x < L/2.)
        return x;
    return -1;
}