# [analisi discussa in "Curve Approssimanti" negli Oggetti Matematici]
# Posizionati nella cartella in cui c'e' il file; se e' il caso
# cambia la riga seguente
dove <- "http://macosa.dima.unige.it/om/voci/curvap"
# ovvero usa  dove <- "C:/.." con l'indicazione della cartella in cui
# hai salvato il file "semilog.txt"
# Carica ed esamina il file
nome <- paste(dove,"semilog.txt",sep="/")
dati <- read.table(nome,sep=",")
dim(dati)
[1] 16  2
min(dati$V1); max(dati$V1); min(dati$V2); max(dati$V2)
[1] 2.5
[1] 77.5
[1] 1e-04
[1] 0.0735
x1 <- -1; x2 <- 80; y1 <- 0; y2 <- 0.08
plot(c(x1,x2), c(y1,y2), type="n",xlab="",ylab="")
abline(v=0,h=0,col="blue")
abline(v=seq(10,80,10),h=seq(0.01,0.08,0.01),lty=3)
points(dati, pch=20)

dati2 <- dati
dati2$V2 <- log(dati2$V2)
min(dati2$V1); max(dati2$V1); min(dati2$V2); max(dati2$V2)
[1] 2.5
[1] 77.5
[1] -9.21034
[1] -2.61047
x1 <- -1; x2 <- 80; y1 <- -10; y2 <- 0
plot(c(x1,x2), c(y1,y2), type="n",xlab="",ylab="")
abline(v=0,h=0,col="blue")
abline(v=seq(10,80,10),h=seq(-10,-1,1),lty=3)
points(dati2); lines(dati2)
points(dati2[1,c(1,2)], pch=20)
points(dati2[7,c(1,2)], pch=20)
lines(c(dati2[1,1],dati2[1,1],dati2[7,1]),c(dati2[1,2],dati2[7,2],dati2[7,2]),col="red")

dati2[1,c(1,2)]
   V1       V2
1 2.5 -2.61047
dati[7,c(1,2)]
    V1     V2
7 32.5 0.0037
(-2.6+5.6)/(-30); 0.25-2.6
[1] -0.1
[1] -2.35
f <- function(x) -0.1*x-2.35
plot(f,0,80,add=TRUE,lty=3,col="red")
exp(-2.35)
[1] 0.09536916