source("http://macosa.dima.unige.it/r.R") tt <- NULL; for(t in (0:24)*5) tt <- c(tt,t) pp <- c(1.02,1.10,1.35,1.53,1.75,2.02,2.12,2.59,2.77,3.17,3.84,4.12,4.89, 5.57,6.16,6.98,8.00,9.32,11.19,12.84,14.30,15.75,18.73,21.22,23.02) Plane(0,120, 0,25) P <- function(t) 1.0268^t graph1(P, 0,120, "blue") POINT(tt,pp, "brown") polyli(c(105,105,0),c(0,16,16),"blue") type(105,1.3,"105") type(4,16,"16") # # ovvero: tt <- NULL; for(t in (0:24)*5) tt <- c(tt,t) pp <- c(1.02,1.10,1.35,1.53,1.75,2.02,2.12,2.59,2.77,3.17,3.84,4.12,4.89, 5.57,6.16,6.98,8.00,9.32,11.19,12.84,14.30,15.75,18.73,21.22,23.02) plot(c(0,120),c(0,25),type="n",xlab="", ylab="") abline(v=axTicks(1), h=axTicks(2), col="blue",lty=3) points(tt,pp) P <- function(t) 1.0268^t plot(P,0,120, add=TRUE, col="red") lines(c(105,105,0),c(0,16,16),lty=2,col="brown")