def discretize(pde, geo, time_nsteps=None, space_nsteps=None):
    """
        Discretize PDE and Geometry

        Parameters:
            pde (PDE): The partial differential equations.
            geo (Geometry): The geometry to be discretized.


        Returns
        --------
        pde_disc: PDE
            Reserved parameter

        geo_disc: DiscreteGeometry
            Discrte Geometry
    """

    # Geometry
    geo_disc = geo.discretize(time_nsteps, space_nsteps)

    return pde, geo_disc