################################################
#
# Make a smooth function for a given gradient values
#
################################################
cleanerfun <- function(x, y0, Corz, Grad){
  val <- y0 + Grad[1] * x
  for( i in 2:(length(Corz)-1)){
    delta <- Corz[i+1] - Corz[i]
    buff  <- ( Grad[i] - Grad[i-1] ) / 2 * ( delta * log( 2 * cosh(( x - Corz[i] )/delta)) + x - Corz[i])
    val   <- val + buff
  }
  return( val )
}

# x : データの値がほしい x 座標
# y0 : x[1] におけるデータの値
# Corz: データの値の傾きを与える点の座標
# Grad: 上記座標でのデータの値の傾き

#
# Examples : give gradient varlues
#
# Delta: 傾きを与える点の間隔
# Corz: 傾きを与える点の座標
# Grad: 上記座標での傾き
#
Delta  <- 0.1
Corz   <- seq(0,10,Delta)+Delta/2
Grad   <- -sin(Corz)

#
# Examples : calculate the value
#
x      <- seq(0,10,0.01)
result <- cleanerfun(x, 1, Corz, Grad)
plot(x, result)