Hydrated soft tissues, such as articular cartilage, are often modeled as biphasic systems with individually incompressible solid and fluid phases, and biphasic models are employed to fit experimental data in order to determine the mechanical and hydraulic properties of the tissues. Two of the most common experimental setups are confined and unconfined compression. Analytical solutions exist for the unconfined case with the linear, isotropic, homogeneous model of articular cartilage, and for the confined case with the non-linear, isotropic, homogeneous model. The aim of this contribution is to provide an easily implementable numerical tool to determine a solution to the governing differential equations of (homogeneous and isotropic) unconfined and (inhomogeneous and isotropic) confined compression under large deformations. The large-deformation governing equations are reduced to equivalent diffusive equations, which are then solved by means of finite difference (FD) methods. The solution strategy proposed here could be used to generate benchmark tests for validating complex user-defined material models within finite element (FE) implementations, and for determining the tissue's mechanical and hydraulic properties from experimental data.