Radiation of a Relativistic Electron

values = {e-> 1, c-> 1, a-> 1};

γ= 1/(Sqrt[1-(β^2)]);

P[θ_, β_]=((e^2)/(4 π (c^3)))((a^2)/(((1-β Cos[θ]))^3))(1 - (((Sin[θ])^2)/((γ^2)(((1-β Cos [θ]))^2))))/.values//FullSimplify

polars=Table[
PolarPlot[
P[θ,β]/.values,
{θ,0,2π},
PlotRange->0.1+β,
PlotLabel-> "β=".β,
Ticks-> {{0.1+β}, {0.1+β}},
ImageSize->{300,300}
],{β,0,1,0.3}
]
regulars=Table[
Plot[
P[θ,β]/.values,{θ,-π,π},
PlotLabel-> "β=".β,
ImageSize->{300,300},
Ticks-> {{-π, ((-π)/2),0, (π/2),π},Automatic}
],{β,0,1,0.3}
]

Export["polar1.png", polars[[1]], "PNG"]
Export["polar2.png", polars[[2]], "PNG"]
Export["polar3.png", polars[[3]], "PNG"]
Export["polar4.png", polars[[4]], "PNG"]

Export["regular1.png", regulars[[1]], "PNG"]
Export["regular2.png", regulars[[2]], "PNG"]
Export["regular3.png", regulars[[3]], "PNG"]
Export["regular4.png", regulars[[4]], "PNG"]

fullpolar = PolarPlot[
P[θ,0.98]/.values,
{θ,0,2π},
PlotRange->{10000},
PlotLabel-> "β=0.98",
Ticks-> {Automatic, {-300, 300}},
ImageSize->{700,100}
]

Export["fullpolar.png",fullpolar, "PNG"]