17   PetscReal beta;  member
30   PetscReal          alpha, beta;  in RHSFunction()  local
40   beta     = user->beta;  in RHSFunction()
55     da = alpha * a * a / (1. + beta * h) + rho_a - mu_a * a;  in RHSFunction()
79   PetscReal          alpha, beta;  in RHSJacobian()  local
90   beta     = user->beta;  in RHSJacobian()
103     va[0]         = 2 * alpha * a / (1. + beta * h) - mu_a;  in RHSJacobian()
106     va[1]         = -beta * alpha * a * a / ((1. + beta * h) * (1. + beta * h));  in RHSJacobian()
215   user.beta     = 1.;  in main()
226   PetscCall(PetscOptionsReal("-beta", "Inhibition factor", "ex42.c", user.beta, &user.beta, NULL));  in main()
237   PetscCall(PetscPrintf(PETSC_COMM_WORLD, "beta:  %5.5g\n", (double)user.beta));  in main()