Added function "check_concave"

This is the function to check if the upper and lower bounds are not actually bounding the density, i.e. the log-concavity is violated. The input arguments are u_k and l_k.