v.ice2 = c(0.09, 0.13, 0.17, 0.21, 0.26, 0.33, 0.45, 0.66, 1.19, NaN)
v.ice10 = c(0.10, 0.25, 0.42, 0.60, 0.78, 1.00, 1.32, 1.78, 2.60, 4.39)
v.dic = c(0.414, 0.771, 1.115, 1.454, 1.798, 2.18, 2.61, 3.12, 3.75, 4.54)

plot(v.dic)
points(v.ice10, col='red')
points(v.ice2, col='blue')


v.mis <- function(v0, N, n, R, RV) {
  rp = n*R*RV / (n*R + RV)
  return( v0 * rp /(rp + (N-n)*R) )
}

plot( v.mis (4.54, 10, seq(1,10), 0.100, 3))
plot( v.mis (4.54, 10, seq(1,10), 0.100, 10))
points( v.mis (4.54, 10, seq(1,10), 0.100, 3), col='red')
points( v.mis (4.54, 10, seq(1,10), 0.100, 1), col='blue')
points( v.mis (4.54, 10, seq(1,10), 0.100, 0.3), col='green')
points( v.mis (4.54, 10, seq(1,10), 0.100, 0.04), col='yellow')