I having problems with MATLAB ode45 solver. For the following code, Matlab is giving me index errors. All I am trying to do is integrate the differential equations in a ode45 function. But I don't see what the problem here is. Can you help me fix this code.   % Initial Conditions w = [2;-1;1]; T = 30; I = [150; 400; 400]; EP = [0;0;0;1]; C = eye(3,3); tspan = 0:20; % Using ode45 to integrate w dot options = odeset('RelTol',1e-10,'AbsTol',1e-10); result = ode45(@K_DDE, tspan, [w';C], options); % Extracting information from ode solver t = result.x; w_ode = result.y(1,:); C_ode = result.y(2:4,:) function dwCdt = K_DDE(~,w_C)     w = w_C(1,:);     C = w_C(2:4,:);     % Kinematic Differentia Equations for DCM     dCdt = eye(3,3);     dCdt(1,1) = C(1,2)*w(3) - C(1,3)*w(2);     dCdt(1,2) = C(1,3)*w(1) - C(1,1)*w(3);     dCdt(1,3) = C(1,1)*w(2) - C(1,2)*w(1);     dCdt(2,1) = C(2,2)*w(3) - C(2,3)*w(2);     dCdt(2,2) = C(2,3)*w(1) - C(2,1)*w(2);     dCdt(2,3) = C(2,1)*w(2) - C(2,2)*w(1);     dCdt(3,1) = C(3,2)*w(3) - C(3,3)*w(2);     dCdt(3,2) = C(3,3)*w(1) - C(3,1)*w(3);     dCdt(3,3) = C(3,1)*w(2) - C(3,2)*w(1);     dwCdt = [w'; dCdt];     end

I having problems with MATLAB ode45 solver. For the following code, Matlab is giving me index errors. All I am trying to do is integrate the differential equations in a ode45 function. But I don't see what the problem here is. Can you help me fix this code.

% Initial Conditions
w = [2;-1;1];
T = 30;
I = [150; 400; 400];
EP = [0;0;0;1];
C = eye(3,3);
tspan = 0:20;

% Using ode45 to integrate w dot
options = odeset('RelTol',1e-10,'AbsTol',1e-10);
result = ode45(@K_DDE, tspan, [w';C], options);

% Extracting information from ode solver
t = result.x;
w_ode = result.y(1,:);
C_ode = result.y(2:4,:)

function dwCdt = K_DDE(~,w_C)

    w = w_C(1,:);
    C = w_C(2:4,:);

    % Kinematic Differentia Equations for DCM
    dCdt = eye(3,3);
    dCdt(1,1) = C(1,2)*w(3) - C(1,3)*w(2);
    dCdt(1,2) = C(1,3)*w(1) - C(1,1)*w(3);
    dCdt(1,3) = C(1,1)*w(2) - C(1,2)*w(1);
    dCdt(2,1) = C(2,2)*w(3) - C(2,3)*w(2);
    dCdt(2,2) = C(2,3)*w(1) - C(2,1)*w(2);
    dCdt(2,3) = C(2,1)*w(2) - C(2,2)*w(1);
    dCdt(3,1) = C(3,2)*w(3) - C(3,3)*w(2);
    dCdt(3,2) = C(3,3)*w(1) - C(3,1)*w(3);
    dCdt(3,3) = C(3,1)*w(2) - C(3,2)*w(1);

    dwCdt = [w'; dCdt];
   

end