please solve with matlab code

c) Now use FDD, BDD, or CDD (whichever is appropriate) to find the 1st and 2nd derivative of Temperature (T) with respect to depth (z) with O(h2). Plot these results separately, and locate the thermocline from d2T/dz2 plot. (Formulations for FDD, BDD, and CDD are given in Figures 21.3, 21.4, and 21.5, respectively). Remember Thermocline is at a depth where (d2T/dz2)=0. Jse cubic splines to determine the thermocline depth for a Lake for which the Temperature vs Depth data is given as follows