pdf("destroyer.pdf")

d0 <- read.fwf("C:\\data\\destroyer.dat", width=c(6,9,7,7,4,4),
col.names=c("stbdrpm", "portrpm", "hpmain", "hpcllct", "extfoul", "modfoul"))

rpm = (d0$stbdrpm + d0$portrpm)/100

foul = d0$extfoul + d0$modfoul

hp = d0$hpcllct/100


d1 <- data.frame(rpm=rpm, foul=foul, hp=hp)

attach(d1)

d1

d1.nl1 <- nls(hp ~ b0*foul + b1*(rpm^b2), start=c(b0=1, b1=1, b2=1))

summary(d1.nl1)

dev.off()