# Imponiamo un seed per avere un esempio riproducibile
set.seed(12345)

# Il numero di punti desiderati
num.points <- 100

# Creaimo dei valori di energia distribuiti normalmente
energy.x <- 1:num.points
energy.y <- rnorm(n=num.points, mean=0, sd=10)
# Creiamo il grafico
plot(energy.x, energy.y, t="o", pch=20, xlab="", ylab="Energia", las=1)
abline(h=0, col="darkgray", lwd=2)

# Dobbiamo trovare le intersezioni della nostra spezzata con l'asse x.
# I punti di intersezione stanno tra punti negativi e positivi
# Moltiplichiamo ciascun valore per il successivo: 
# quando il segno cambia il prodotto sarà negativo
# Questa linea ci trova tutti i punti prima dell'intersezione
crossings <- energy.y[-length(energy.y)] * energy.y[-1]
crossings <- which(crossings < 0)

# Possiamo disegnarli per confema (scommentare la linea successiva)
#  points(energy.x[crossings], energy.y[crossings], col="red", pch="X")

# Ora ci passiamo tutti i punti prima dell'intersezione e troviamo 
# l'effettiva intersezione usando una proporzione. 
# Visto che i problemi di geometria delle medie servivano a qualcosa?
intersections <- NULL
for (cr in crossings)
  {
  new.int <- cr + abs(energy.y[cr])/(abs(energy.y[cr])+abs(energy.y[cr+1]))
  intersections <- c(intersections, new.int)
  }

# Confermiamo che le intersezioni siano corrette (scommentare la linea successiva)
#  points(intersections, rep(0, length(intersections)), pch=20, col="red", cex=0.7)

last.intersection <- 0
for (i in intersections)
  {
  ids <- which(energy.x<=i & energy.x>last.intersection)
  poly.x <- c(last.intersection, energy.x[ids], i)
  poly.y <- c(0, energy.y[ids], 0)
  if (max(poly.y) > 0)
    {
    col="green"
    }
  else
    {
    col="red"
    }
  polygon(x=poly.x, y=poly.y, col=col)

  last.intersection <- i
  }

o--|
 \ |
  \|
---x--------
   |\
   |_o

new.int <- cr + abs(energy.y[cr])/(abs(energy.y[cr])+abs(energy.y[cr+1]))

abs(10) / (abs(10) + abs(-5)) = 10 / 15 = 2/3

  new.int <- energy.x[cr] + abs(energy.y[cr])/(abs(energy.y[cr])+abs(energy.y[cr+1])) *
    (energy.x[cr+1]-energy.x[cr])

grafico2.pdf

grafico3.pdf

# vettori che contengono le coordinate x e y dei vari punti
# le ordinate sono generate casualmente grazie alla funzione sample
vx<-c()
vy<-c()

for(i in 1:100){
 vx[i]=i
 vy[i]=sample(-10:10,1)
}

# il grafico viene disegnato

plot(vx,vy,main="Grafico di idropatia",xlab="Amminoacidi",ylab="Energia libera",type="l",col.axis="blue",las=1)
abline(h=0)


# ricerca di tutti i punti che precedono un'intersezione con l'asse delle x

int_x <- c()
int_y <- c()
d = 1 
int_x[d]=0 # inserisco un primo punto in 0,0 (mi serve per colorare la prima area)
int_y[d]=0
for(i in 1:(100-1)){
   if((vy[i]*vy[i+1]) <= 0){ # relazione da te suggeritami
     d=d+1 
     int_x[d]=vx[i]
     int_y[d]=vy[i]
    }
} 

int_x[d+1]=vx[100]
int_y[d+1]=0
points(int_x,int_y,col="red",pch=20)

# ricerca di tutte le intersezioni sull'asse delle x 

num_int = length(int_x) #ci si può mettere anche int_y
int_pol_x <- c()
int_pol_y <- c()
for(i in 1:num_int){
   int_pol_x[i] = ( int_x[i] + ( abs(int_y[i])/( abs(int_y[i])+abs(int_y[i+1]) ) ) )
   int_pol_y[i] = 0
}
points(int_pol_x,int_pol_y,col="blue",pch=20)

vx<-c()
vy<-c()

for(i in 1:100){
 vx[i]=i
 vy[i]=sample(-10:10,1)
}

vx = 1:100
vy = sample(-10:10, 100, replace=TRUE)

vy = sample(seq(-10,10,0.01), 100, replace=TRUE)