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")