No Title  |POLVER05_0 |1
d(X) / d(W) = -rA/FA0 
X(0) = 0
d(p) / d(W) = -alpha/2/p*(1+epsilon*X) 
p(0) = 1
CA=CA0*(1-X)*p/(1+epsilon*X)
epsilon=-0.5
CA0=0.2
k=5 
alpha=0.0162 
FA0=2.5 
rA=-k*CA^2
W(0) = 0
W(f) = 100