No Title  |POLVER05_0 |1
d(X)/d(W) = -ra / Fa0
X(0) = 0
ra = -k * (Ca * Cb - (Cc ^ 3 / Kc))
Ca = Ca0 * (1 - X) / (1 + eps * X) * y
Cb = Ca0 * (thetaB - X) / (1 + eps * X) * y
Cc = 3 * Ca0 * X / (1 + eps * X) * y
d(y)/d(W) = -alpha * (1 + eps * X) / (2 * y)
y(0) = 1
k = 1000
alpha = 0.01
thetaB = 2
eps = 0.33
Ca0 = 0.01
Fa0 = 15
W(0) = 0
W(f) = 90
Kc = 0.05