def critical_pressure(Po, uni_comp_str, k):
    """
    Calculates and returns the critical pressure given Uniaxial Compressive Strength, Pressure (Vertical/Overburden), and k value.
    """
    return (2*Po-uni_comp_str)/(1+k)