
function [dy] = fun(y,tt)

y_1 = y(1);
y_2 = y(2);

x=tt;

dy_1 = -16*y_1+12*y_2+16*cos(x)-13*sin(x);
dy_2 = 12*y_1-9*y_2-11*cos(x)+9*sin(x);

dy = [dy_1,dy_2];