Uso il file descritto qui.
x <- c(55,70,155,160,155,115,105,165,110,115,85,165,110,155,105) y <- c(85,115,205,230,185,185,145,240,140,155,125,290,170,195,135) range(x,y) # asc 55 165 ord 85 290 N 15 plot(c(0,170),c(0,300),type="n") abline(h=0,v=0,lty=2,col="blue") abline(h=axTicks(2),v=axTicks(1),lty=3,col="blue") points(x,y) yy <- function(U,x) U*x; Fun <- function(x) yy(U0,x) dif <- 1e300 m1 <- 0; M1 <- 100; calc1(x,y,1e4); dif; U0 # 6757.679 1.425915 # Proseguo restringendo l'intervalli per U; aumento le prove m1 <- 0; M1 <- 10; calc1(x,y,1e4); dif; U0 # 6750.542 1.420354 m1 <- 1.41; M1 <- 1.43; calc1(x,y,1e4); dif; U0 # 6750.539 1.420462 curve(Fun, add=TRUE, col="red")