Please come out with a Matlab code solution

9.4. In the quenching operation used to harden a metal component, the compone or bar is heated to a high temperature 7 and immersed in an oil (or water) bath at temperature Tu as shown in Fig. 9.19. The consideration of ther equilibrium of the bar and the oil lead to the equations [9.15] and ha di myc dTu hA di + Tb = Tw = (EZ) = with initial conditions T(0) Tio and T(0) T20, where Tb(x) = temperature of the bar (oil), my (m) = mass of the bar (oil), cs (c) = specif heat of the bar (oil), h = convection heat transfer coefficient and A area of the bar. surface (E1) + Tw=Tb Determine the variations of T, and Tu as functions of time, f, for 0 for the following data: h = 500 BTU/hr-ft-F, A=0.1 ft, m = 0.02 lb, c = 0.22 BTU/lb - F. T10 = 1200F, T20 = 65F, m=2.5 lb, and cw = 1.00 BTU/lb - F. 10 sec