def PicardPlot(ax, s, c, d):
    """
    s : array of nonincreasing singular values
    c : array of coefficients |<u_n, b>|
    d : array of coefficients for the solution |<u_n, b>| / sigma_n
    """
    ax.plot(s, '.', clip_on=False, markersize=9, label='$\\sigma_n$', color
        =ax.pointcolor)
    ax.plot(c, 'dm', clip_on=False, alpha=ax.alpha, label=
        '$\\vert\\langle \\phi_n, u^i_{\\mathbf{z},\\tau} \\rangle\\vert$')
    ax.plot(d, 'Pc', clip_on=False, markersize=7, alpha=ax.alpha, label=
        '$\\vert\\langle \\phi_n, u^i_{\\mathbf{z},\\tau} \\rangle\\vert/\\sigma_n$'
        )
    ax.legend(loc='upper center', bbox_to_anchor=(0.5, 1.175), ncol=3,
        fancybox=True, shadow=True, fontsize='large')
    ax.set_xlabel('n', color=ax.labelcolor)
    ax.set_xlim([0, len(s)])
    ax.set_ylim(bottom=0)
    ax.locator_params(axis='y', nticks=6)
    ax.ticklabel_format(style='sci', axis='y', scilimits=(0, 0))
    return ax
