Displaying 1 result from an estimated 1 matches for "delz".
Did you mean:
dell
2013 Apr 10
0
Problem with ode
...r){
if (hr <=16){
450*sin(hr*pi/16)
} else {
0
}
}
i = which(z == DEPTH)
I = I0(hr=hr)*0.43*exp(-z*Kw - Kc*(cumsum(delz*P)[i]
- delz*P[i]/2))
FI = alpha*I/(um^2 + (alpha*I)^2)^.5
return(FI)
}
J1 = function(x){
n = 10^3
FI2 <- function(y)F.I(z = x, y)
J1 = sum(sapply(runif(n, 0, 16), FI2))*16/24/n
return(J1)
}
J = sapply(DEPTH,...