please solve with matlab code

d) If the heat flux from the surface to the bottom layer can be computed with Fourier's LAW J=kdzdT Compute the flux at this thermocline interface (k=0.01cal/(s.cm.C)) NOTE: Please check "help spline" and "help diff" to get more information on using built in functions Jse cubic splines to determine the thermocline depth for a Lake for which the Temperature vs Depth data is given as follows