The rotated parabolic equation [J. Acoust. Soc. Am. 87, 1035-1037 (1990)] is generalized to problems involving ocean-sediment interfaces of variable slope. The approach is based on approximating a variable slope in terms of a series of constant slope regions. The original rotated parabolic equation algorithm is used to march the field through each region. An interpolation-extrapolation approach is used to generate a starting field at the beginning of each region beyond the one containing the source. For the elastic case, a series of operators is applied to rotate the dependent variable vector along with the coordinate system. The variable rotated parabolic equation should provide accurate solutions to a large class of range-dependent seismo-acoustics problems. For the fluid case, the accuracy of the approach is confirmed through comparisons with reference solutions. For the elastic case, variable rotated parabolic equation solutions are compared with energy-conserving and mapping solutions.