Synchrotron Bending Magnet Radiation

plot=LogLogPlot[1.33*(10^4)*(((1.7))^2)*200*(10^(-3))(x^2)((BesselK[(2/3),(x/2)])^2),{x,0.0001,15}, AxesLabel->{Style["E/(Subscript[E, c])",14],Style["\!\(\*FractionBox[\(\*SuperscriptBox[\(d\), \(3\)]\\ \*SubscriptBox[\(F\), \(B\)]\), \(dΘ\\ dϕ\\ dω/ω\)]\)",14],}, ImageSize->500]

Assuming[x>0,{G1[x_]=x*(∫_x^∞)BesselK[(5/3),y]\[DifferentialD]y}]

plot2=LogLogPlot[2.46*(10^13)*(((1.7))^2)*200*(10^(-3))*G1[x],{x,0.001,15},PlotRange->{(10^9),2*(10^13)}, AxesLabel->{Style["E/(Subscript[E, c])",14],Style["Flux",14]}, ImageSize->500]

Flux is in ((photons/s)/(0.1% BW))

Export["flux.png", 
  plot2, "PNG"]