def METHOD_NAME(FC, fill=True, key='Default',
                         title='Fragility curve', 
                         xmin=0, xmax=1, npoints=100, 
                         xlabel='x', 
                         ylabel='Probability of exceeding a damage state',
                         ax=None):
    """
    Plot fragility curve.
    Parameters
    -----------
    FC : wntr.scenario.FragilityCurve object
        Fragility curve
    fill : bool, optional
        If true, fill area under the curve
    key : string, optional
        Fragility curve state distribution key
    title : string, optional
        Plot title
    xmin : float, optional
        X axis minimum 
    xmax : float, optional
        X axis maximum 
    npoints : int, optional
        Number of points 
    xlabel : string, optional
        X axis label
    ylabel : string, optional
        Y axis label
    ax : matplotlib axes object, optional
        Axes for plotting (None indicates that a new figure with a single 
        axes will be used)
    Returns
    ---------
    ax : matplotlib axes object
    """
    if ax is None: # create a new figure
        plt.figure(figsize=[8,4])
        ax = plt.gca()
    ax.set_title(title)
    x = np.linspace(xmin,xmax,npoints)
    for name, state in FC.states():
        try:
            dist=state.distribution[key]
            if fill:
                plt.fill_between(x, dist.cdf(x), label=name)
            else:
                plt.plot(x, dist.cdf(x), label=name)
        except:
            pass
    ax.set_xlim((xmin,xmax))
    ax.set_ylim((0,1))
    ax.set_xlabel(xlabel)
    ax.set_ylabel(ylabel)
    ax.legend()
    plt.show(block=False)
    return ax