# Finding the functions that Mathematica plots while using BodePlot

I've the following Mathematica-code:

BodePlot[{1/(((5.601443110900379*^-8)*(5.601443110900379*^-8)*(1.\
6672176431702952*^-7)*(56000)*(56000 +
56000)*((56000*56000)/(56000 + 56000)))*(2*Pi*f*
I)^3 + ((5.601443110900379*^-8)*(5.601443110900379*^-8)*(\
56000)*(56000 +
56000) + \
(5.601443110900379*^-8)*(1.6672176431702952*^-7)*(56000)*(56000 +
56000))*(2*Pi*f*
I)^2 + ((5.601443110900379*^-8)*(56000) + \
(5.601443110900379*^-8)*(56000 + 56000 + 56000))*(2*Pi*f*I) +
1)}, {f, 0, 100}, PlotStyle -> {Black},
ScalingFunctions -> {{"Linear", "dB"}, {"Linear", "Degree"}},
FrameLabel -> {{" ", " "}, {HoldForm[Text[Frequentie[f]]],
HoldForm[Text[Fase[\[Degree]]]]}}, ImageSize -> Large]
`

And it give two plots: one of the amplitude and one the phase.

Now, I want to know what functions Mathematica is using to plot that phase. How can I found that out? So for the phase it is using the $$\arctan$$ functions but I do not know with what boundaries it is using it and how the function mathematically can be expressed, so can I ask Mathematica for that?