Matlab - Problema flujo couette metodo Crank Nicolson

Hola! Tengo que resolver el problema de flujo de Couette con el método de Crank Nicolson de diferencias finitas y mi nivel de Matlab es muy básico, me vendría bien cualquier ayuda sobre como declarar las condiciones iniciales y, sobre todo, como como programar la matriz para resolver las ecuaciones.

Un saludo y gracias,