# przykład 2.3

plot(mmax,log10(perf))
cpus.lm.m<-lm(log10(perf)~mmax)
abline(cpus.lm.m)

xp<-seq(64.0,64000,length=25)
yp<-seq(0.0,256,length=25)
cpusT<-expand.grid(mmax=xp,cach=yp)

cpus.lm<-lm(log10(perf)~mmax+cach,cpus)
cpusT$perf<-predict(cpus.lm, newdata=cpusT,type="response")
pp<-matrix(cpusT$perf,25,25)
persp(xp,yp,pp,theta = -45, phi = 30, expand = 0.8, d=2, 
ltheta = 30, shade = 0.75,ticktype="detailed",nticks=6,xlab = "mmax",ylab = "cach",zlab = "log10(perf)")->rr
points(trans3d(cpus$mmax, cpus$cach, 10, pm = rr), col = 2, pch =16)