from math import sin, cos, pi


def coordinates_degree_to_stereographic(Lambda_degree, Phi_degree):
    """
    Converts geographic coordinates [°] to cartesian coordinates [m] in stereographic RADOLAN projection.
    
    :Parameters:
    ------------
    
        Lambda_degree : float
            Degree of latitude [°N / °S]
        Phi_degree : Float
            Degree of longitude [°E / °W]
    
    :Returns:
    ---------
    
        (x, y) : Tuple with two elements of type float
            Cartesian coordinates x and y in stereographic projection [m]
            
    """
    
    # Earth Radius in m
    R = 6370.04 * 1000
    
    # Convert decimal degrees to radian
    Phi = Phi_degree * pi/180 
    Lambda = Lambda_degree * pi/180
    
    # Phi0 = 60°N --> Plane of projection subtends terrestrial sphere at 60°N 
    # Lambda0 = 10°E --> Cartesian coordinate system is aligned at 10°E meridian
    phi0 = 60 * pi/180
    lambda0 = 10 * pi/180
    # M = Stereographic Scaling Factor
    M = (1 + sin(phi0))/(1 + sin(Phi)) 
    
    x = R * M * cos(Phi) * sin(Lambda - lambda0)
    y = -R * M * cos(Phi) * cos(Lambda - lambda0)
    return (x, y)