Please see the below script I have created to solve the problem. It doesn't work. Can some one help?
clear; close all; clc;
syms x1(t) x2(t) x3(t) x4(t);
r=[-0.05 0 0.03 0.05;
0.05 -0.05 0 0;
0 0.02 -0.02 0;
0 0.03 0.02 -0.05];
D=[0;1;4;0];
X=[x1;x2;x3;x4];
odes=diff(X)==r*X+D...
