Write a MATLAB program that computes the solution to the (dimensionless) unsteady diffusion equation

subject to an initial condition c(x, 0) = 0 and the boundary conditions c(0, t) = 1 and ∂c/∂x = 0 at x = 1. In other words, this is the solution for the concentration (or temperature profile) in an initially empty (cold) slab if we instantaneously put it in contact with a reservoir at the left that is at the maximum concentration (temperature). Use centered finite differences for the spatial derivative and RK4 for the time integration. You should have 51 nodes and a time step of 0.0001. Make a plot that has the concentration profile at the 11 dimensionless times t = 0, 0.025, 0.05, . . . , 0.25.

at || a2c 2