Use matlab

Consider steady-state one-dimensional heat conduction in a wall with specified temperatures on the two wall boundaries (see below). The governing equation is dx2d2T=kq, where T is the temperature (C),q is the heat source strength (W/m3), and k is the thermal conductivity (W/mC). BC (Boundary conditions) are: x=L,T=TL=30Cx=+L,T=T+L=10C The analytical solution to the governing equation with the boundary conditions applied is T(x)=2kq(L2x2)+2LT+L(L+x)+2LTL(Lx) The analytical solution is the exact solution to this case. Using a finite difference approximation of the equation, use MATLAB to derive the numerical and the analytical solution results for the following parameters: k=0.3W/roC,L=0.2m,q=1000W/m3 Compare the two solutions for several values of the numerical grid cell size. For guidance, try the following grid cell sizes: x=0.1m,0.05m,0.01m,0.001m. Show your result graphically. An illustration of the shape of the temperature profile (from the analytical solution) for various values of Q is shown below. Example profiles