Crank-Nicolson method (Matlab)

[intel2]

Hello,

can someone help me write a matlab code for the Crank-Nicolson method using the attached integration equation:


thank you