Verification of Steefan's Law

cat("Steefans Law");
voltage <- c(seq(25, 47, 0.4));voltage;
length(voltage)
current <- c(62,62,62,62,62,63,63,63,63,63,63,63,64,64,64,64,64,64,64,65,65,65,65,65,65,65,66,66,66,66,66,66,67,67,67,67,67,67,68,68,68,68,68,69,69,69,69,69,69,69,69,69,70,70,70,70);current;
length(current)
resistance <- c(voltage/current);resistance
power <- (current * voltage); power
LogP <- c(log10(power));LogP
LogR <- c(log10(resistance));LogR
dataSet <- data.frame(voltage,current,resistance,power,LogP,LogR);dataSet
plot(LogR, LogP, "l",col = "green", main = "Verification of Steefan's Law");
fitLM <- lm(LogP ~ LogR);fitLM
slope <- coef(fitLM)[2];slope;
abline(fitLM)