def gamma_get_shape_scale(mean,stdev):
    """ Getting gamma distribution shape and scale
    
    Parameters
    ----------
    mean (float): mean of the gamma distribution
    stdev (float): stander deviation of the gamma distribution
    
    Returns
    ----------
    shape (float): shape of the gamma distribution
    scale (float): scale of the gamma distribution
    """
    shape = (mean**2)/(stdev**2)
    scale = (stdev**2)/mean
    
    return shape,scale