Please solve with matlab code Consider an IVP for the diffusion equation in one dimension: Ju(x,t) at J²u(x,t) = D дх² on the interval x = [0, Z] with initial condition u(x, 0) = f(x), VxE [0, L] and Dirichlet boundary conditions u(0,t) = u(L,t) = 0 Vt > 0. Write Formulation and Equations than solve with MATLAB codes and Show model and analyze