d_Z3_free = + k_Z2_Z3_free * Z2 - k_Z3_free_sink * Z3_free - k_Z3_free_bound * Z3_free + k_Z3_bound_free * Z3_bound
d_Z3_bound = + k_Z3_free_bound * Z3_free - k_Z3_bound_free * Z3_bound

d_Z0_free = - 0 - k_Z0_free_bound * Z0_free + k_Z0_bound_free * Z0_bound - k_Z0_free_Z1 * Z0_free
d_Z0_bound = + k_Z0_free_bound * Z0_free - k_Z0_bound_free * Z0_bound
d_Z1 = + k_Z0_free_Z1 * Z0_free - k_Z1_sink * Z1

d_Z0_free = - 0 - k_Z0_free_bound * Z0_free + k_Z0_bound_free * Z0_bound - k_Z0_free_Z1 * Z0_free
d_Z0_bound = + k_Z0_free_bound * Z0_free - k_Z0_bound_free * Z0_bound
d_Z1 = + k_Z0_free_Z1 * Z0_free - k_Z1_sink * Z1 - k_Z1_Z2 * Z1
d_Z2 = + k_Z1_Z2 * Z1 - k_Z2_sink * Z2

d_Z0_free = - 0 - k_Z0_free_bound * Z0_free + k_Z0_bound_free * Z0_bound - k_Z0_free_Z1 * Z0_free
d_Z0_bound = + k_Z0_free_bound * Z0_free - k_Z0_bound_free * Z0_bound
d_Z1 = + k_Z0_free_Z1 * Z0_free - 0 - k_Z1_Z2 * Z1
d_Z2 = + k_Z1_Z2 * Z1 - k_Z2_sink * Z2 - k_Z2_Z3_free * Z2
d_Z3_free = + k_Z2_Z3_free * Z2 - k_Z3_free_sink * Z3_free - k_Z3_free_bound * Z3_free + k_Z3_bound_free * Z3_bound
d_Z3_bound = + k_Z3_free_bound * Z3_free - k_Z3_bound_free * Z3_bound

Relative error in the sum of squares is at most `ftol'.

d_Z0_free = - 0 - k_Z0_free_bound * Z0_free + k_Z0_bound_free * Z0_bound - k_Z0_free_Z1 * Z0_free
d_Z0_bound = + k_Z0_free_bound * Z0_free - k_Z0_bound_free * Z0_bound
d_Z1 = + k_Z0_free_Z1 * Z0_free - k_Z1_sink * Z1

d_Z0_free = - 0 - k_Z0_free_bound * Z0_free + k_Z0_bound_free * Z0_bound - k_Z0_free_Z1 * Z0_free
d_Z0_bound = + k_Z0_free_bound * Z0_free - k_Z0_bound_free * Z0_bound
d_Z1 = + k_Z0_free_Z1 * Z0_free - k_Z1_sink * Z1 - k_Z1_Z2 * Z1
d_Z2 = + k_Z1_Z2 * Z1 - k_Z2_sink * Z2

d_Z0_free = - 0 - k_Z0_free_bound * Z0_free + k_Z0_bound_free * Z0_bound - k_Z0_free_Z1 * Z0_free
d_Z0_bound = + k_Z0_free_bound * Z0_free - k_Z0_bound_free * Z0_bound
d_Z1 = + k_Z0_free_Z1 * Z0_free - 0 - k_Z1_Z2 * Z1
d_Z2 = + k_Z1_Z2 * Z1 - k_Z2_sink * Z2 - k_Z2_Z3_free * Z2
d_Z3_free = + k_Z2_Z3_free * Z2 - k_Z3_free_sink * Z3_free - k_Z3_free_bound * Z3_free + k_Z3_bound_free * Z3_bound
d_Z3_bound = + k_Z3_free_bound * Z3_free - k_Z3_bound_free * Z3_bound