import math


def cylinder_plane(F, L, D1, E1, E2):
    """
    Calculate line contact pressure a cylinder and a plane.

    Arguments;
        F: Force, e.g. in [N].
        L: Contact length, e.g. in [mm].
        D1: Diameter of cylinder in e.g. [mm].
        E1: Young's modulus of first cylinder in e.g. [MPa].
        E2: Young's modulus of the plane in e.g. [MPa].

    Returns:
        Maximum contact pressure, e.g. in [MPa]
    """
    A = F / (2.86 * L)
    B = 1 / D1
    C = 2 * E1 * E2 / (E1 + E2)
    Pmax = math.sqrt(A * B * C)
    return Pmax