The solutions of the space–time fractional diffusion equations and that of the space–
time fractional Fokker–Planck equation are probabilities evolving in time and stable in the sense
of Lévy. The fundamental solution, Green function, of the space–time fractional diffusion equation,
is early obtained by using the scale invariant method. In this paper, I use this reduced Green
functions and the scale invariant method to obtain the fundamental solution, Green function, of
the fractional diffusion equation and henceforth I obtain the solution of the space–time fractional
Fokker–Planck equation, by applying the Biller ´s transformation between the independent spatial coordinates
of these fractional differential equations. Henceforth, I simulate these solutions in the 3D
for all the possible values of the space and time fractional orders and also for different values of the
skewness.