programming a PDE in Matlab

Who can help me with programming a partial differential equation?

The equation is: d(Ci)/dt = b * 1/(n^2) * d/d(n) * ( n^2 * d (Ci)/dn )

- With Ci being a concentration internally. - b some sort of constant - n the dimensionless variable for a radial distance, varying from 0 to 1.

Answer 1

[Ivan van der Kroon]

You already have written it in a ode-form, i.e. first temporal derivative on LHS. The RHS only needs the correct implementation of the spatial derivatives. You can use functions like gradient, diff or spectral methods to compute these.