def METHOD_NAME(output1, output2, scalefactor):
    """Make the plots"""
    matpow1 = os.path.join(output1, "powerspectrum-%.4f.txt" % scalefactor)
    pk1 = get_power(matpow1)
    matpow2 = os.path.join(output2, "powerspectrum-%.4f.txt" % scalefactor)
    pk2 = get_power(matpow2)
    ax_abs, ax2_rel = make_figures()
    ax_abs.loglog(pk1[0], pk1[1], label=output1)
    ax_abs.loglog(pk2[0], pk2[1], label=output2)
    ax2_rel.semilogx(pk2[0], pk2[1]/pk1[1])
    ax_abs.legend()
    plt.savefig("powerspectrum-%.4f.pdf" % scalefactor)
    plt.clf()