data = read.table("cooling.data")
png("plot.png")

hnu = 4.138e-10;
A = 2.20e-6;
alpha = 7.10e5;
hnu_T = 1.4E4;
EinK =  29168;

plot(data,
    log="xy",
    xlim=c(1e3,1e5),
    ylim=c(1e-18,1e-17),
    type="l",
    xlab="Temperature [K]",
    ylab=" [Cooling Rate ergs cm⁻³ s⁻¹]"
    )

abline(v=hnu_T)

dev.off()
hnu_T