create ode solver for 2 dof in matlab

main.m
% Define the equations of motion for your system
function dYdt = equationsOfMotion(t, Y)
    m1 = 1; % mass of first DOF
    m2 = 1; % mass of second DOF
    k1 = 10; % spring constant of first DOF
    k2 = 10; % spring constant of second DOF
    c1 = 1; % damping coefficient of first DOF
    c2 = 1; % damping coefficient of second DOF
    
    % Unpack the state vector
    y1 = Y(1);
    y2 = Y(2);
    dy1dt = Y(3);
    dy2dt = Y(4);
    
    % Define the equations of motion
    ddy1dt = (-k1*y1 - c1*dy1dt + k2*(y2-y1) + c2*(dy2dt-dy1dt))/m1;
    ddy2dt = (-k2*(y2-y1) - c2*(dy2dt-dy1dt))/m2;
    
    % Pack the derivatives into a vector
    dYdt = [dy1dt; dy2dt; ddy1dt; ddy2dt];
end

% Set initial conditions and time span
Y0 = [0; 1; 0; 0]; % initial conditions
tspan = [0 10]; % time span

% Solve the ODE using MATLAB's built-in solver
[t,Y] = ode45(@equationsOfMotion, tspan, Y0);

% Plot the results
plot(t,Y(:,1),t,Y(:,2));
xlabel('Time (s)');
ylabel('Displacement (m)');
legend('DOF 1', 'DOF 2');
1014 chars
36 lines