Kindly help me out here!!

I have a system can be described by the second order differential equation in 1, where h is the convective heat transfer coefficient, p is the perimeter (circumference) of the pin, k is the conductive heat transfer coefficient and Ac is the cross-sectional area.

I want to develop a numerical model by applying Euler method to the differential equation in equation 1. I want to do this on matlab.