Plot3D[(1/gamma(|alpha|)) y^(alpha-1) e^(-y), {y, 0, 5}, {α, -0, 5}]