Consider the following code:
\[Psi][{n_, l_, m_}, {r_, \[Theta]_, \[Phi]_}] :=
With[{\[Rho] = 2 r/(n a)},
Sqrt[(2/(n a))^3 (n - l - 1)!/(2 n (n + l)!)] Exp[-\[Rho]/
2] \[Rho]^
l LaguerreL[n - l - 1, 2 l + 1, \[Rho]] SphericalHarmonicY[l,
m, \[Theta], \[Phi]]];
factor[n_] := 15 n;
n = 3; l = 1; m = 0; a = 1;
For[n = 1, n < 3, n++,
For[l = 0, l < n, l++,
For[m = 0, m <= l, m++,
g = DensityPlot[
4 \[Pi] (x^2 + y^2 +
z^2) (Abs@\[Psi][{n, l, m}, {Sqrt[x^2 + y^2 + z^2],
ArcTan[z, Sqrt[x^2 + y^2]],
ArcTan[x, y]}])^2 /. {y -> 0}, {x, -factor[n] a,
factor[n] a},
{z, -factor[n] a, factor[n] a}, Mesh -> False, Frame -> False,
PlotPoints -> 200, ColorFunctionScaling -> True,
ColorFunction -> "SunsetColors",
PlotRange -> All, PlotLabel ->
Style[
"(" <> ToString[n] <> "," <> ToString[l] <> "," <>
ToString[m] <> ")",
FontSize -> 24]];
Print[g];
]
]
]
If I change PlotRange->All to PlotRange->Automatic, I'll get my favorite size, but the quality of output is not as well as that with "All". But the size of plot is too small with the first one. Is there any way to take an output with the size of "Automatic" but with the quality of "All"?

PlotRange -> {## & @@ ({#, #} &@{-#, #} &@factor[n] a/2), All}? $\endgroup$Automatic,All,Fulletc work. $\endgroup$