crank nicolson implementation

I am in need of help with an implementation of Crank Nicolson using MATLAB for a project. Time is of the essence, so the work must be completed within a week. I have specific reference materials and sources that I will provide the freelancer, to ensure a successful outcome. The language should be MATLAB. If you have experience and a solid understanding of this programming language and the theory of Crank Nicolson then please reach out and let me know. I want the equations attached solved for theta, theta_c and u.