[u, v, w, y]

c = 1 - y
wrate = 21000 * (1 - c) * exp (- 26.7 / (1 + 3 * c))
dcdx = dx (y)