def METHOD_NAME(theta:float, unit='deg'):
    """
    PTM of rotation of theta degrees along the X axis
    """
    if unit=='deg':
        theta = np.deg2rad(theta)
    Y = np.array([[1, 0, 0, 0],
                  [0, np.cos(theta), 0, np.sin(theta)],
                  [0, 0, 1, 0],
                  [0, -np.sin(theta), 0, np.cos(theta)]], dtype=float)
    return Y