def METHOD_NAME():
    """Validate zoeppritz and approximations at zero incident angle"""
    # Create composite parameters
    ai0, si0, _ = vp0 * rho0, vs0 * rho0, vp0 / vs0
    ai1, si1, _ = vp1 * rho1, vs1 * rho1, vp1 / vs1
    # Zoeppritz
    rpp_zoep = zoeppritz_pp(vp1, vs1, rho1, vp0, vs0, rho0, theta[0])
    rpp_zoep_approx = approx_zoeppritz_pp(vp1, vs1, rho1, vp0, vs0, rho0, theta[0])
    # Aki Richards
    rvp = np.log(vp0) - np.log(vp1)
    rvs = np.log(vs0) - np.log(vs1)
    rrho = np.log(rho0) - np.log(rho1)
    G1, G2, G3 = akirichards(theta[0], vs1 / vp1)
    rpp_aki = G1 * rvp + G2 * rvs + G3 * rrho
    # Fatti
    rai = np.log(ai0) - np.log(ai1)
    rsi = np.log(si0) - np.log(si1)
    G1, G2, G3 = fatti(theta[0], vs1 / vp1)
    rpp_fatti = G1 * rai + G2 * rsi + G3 * rrho
    assert_array_almost_equal(rpp_zoep, rpp_zoep_approx, decimal=3)
    assert_array_almost_equal(rpp_zoep, rpp_aki, decimal=3)
    assert_array_almost_equal(rpp_zoep, rpp_fatti, decimal=3)