Уравнение Шредингера для нестационарного гамильтониана имеет вид:

Я пытаюсь реализовать решатель для уравнения Шредингера для нестационарного гамильтониана вode45, Однако, поскольку гамильтониан $ H (t) $ зависит от времени. Я не знаю, как сделать интерполяцию вode45, Можете ли вы дать мне несколько советов?

psi0 = [0 1];
H = [1 0;0 1]*cos(t); %this is wrong, I do not know how to implement this and pass it to ode45
hbar = 1;
t    = [0:1:100];
[T, psi] = ode45(dpsi, t, psi);
function dpsi = f(t, psi, H, psi0)
dpsi = (1/i)*H*psi;

Я также пытаюсь найти решение матричной интерполяции вMATLAB: интерполяция с использованием матрицы.