Perhaps you can do this way:
tspan = [0 0.1];
x0 = [2 1.5 1 0.5];
[t, x] = ode45(@odefcn, tspan, x0);
plot(t, x), grid on,
xlabel({'$t$'}, 'interpreter', 'latex')
ylabel({'$\mathbf{x}(t)$'}, 'interpreter', 'latex')
legend({'$x_{11}(t)$', '$x_{12}(t)$', '$x_{21}(t)$', '$x_{22}(t)$'}, 'interpreter', 'latex', 'fontsize', 14, 'location', 'northwest')
function xdot = odefcn(t, x)
xdot = zeros(4, 1);
A1 = [0 2; -2 0];
A2 = diag([0 2]*[1; 1]*[15 20]);
A = [A1 zeros(2); zeros(2) A2];
xdot = A*x;
end
