Crank Nicolson MATLAB

Question

Есть ли в MATLAB встроенный решатель, который решает данный ODE (порядок>1) с использованием метода Кранка-Николсона? (Обратите внимание, что у меня есть система ODE, поэтому я мог бы решить каждый из них, используя Crank-Nicolson. Я пробовал ode45, но не получаю правильных результатов)