double b(double a, double fi) { double arg = a - fi; if (fabs(arg) < 1e-15) exit(0); return cos(arg) / sin(arg); }