In this talk I will discuss the construction of approximate solutions for the Non-linear Fokker-Planck equation. We utilize the $L^2$-Wasserstein gradient flow structure of this PDEs to perform a semi discretization in time by means of the variational BDF2 method. Our approach can be considered as the natural second order analogue of the Minimizing Movement or JKO scheme. In comparison to our own recent work on constructing solutions to $\lambda$-contractive gradient flows in abstract metric spaces, the technique presented here exploits the differential structure of the underlying $L^2$-Wasserstein space. We directly prove that the obtained limit curve is a weak solution of the non-linear Fokker-Planck equation without using the abstract th...