SOLVED:Write a MATLAB program that uses the method of lines, RK4, and finite differences to solve the diffusion equation (∂c)/(∂t)=(∂^2 c)/(∂x^2) subject to an initial condition c(x, 0)=1 and boundary conditions c(0, t)=1 and c(1, t)=0. Your program should automatically plot the solution at the times t=0.0005, 0.005,0.02,0.05,0.1, and 0.5 .

Runge-Kutta method (Order 4) for solving ODE using MATLAB