Plot of the numerical solutions to (54) given for various values of the power-law parameter k. The other parameters are fixed as ϵ = 1 and γ = 1, while boundary conditions are taken as ρ(0) = 1, ρ′(0) = 0. As we increase k, the solutions uniformly decrease in value, more rapidly tending toward negative infinity as η becomes large.