simulate a mass in between two springs connected to rigid frames in matlab

main.m
% System parameters
m = 1; % mass (kg)
k1 = 10; % spring constant 1 (N/m)
k2 = 5; % spring constant 2 (N/m)
L1 = 1; % rest length of spring 1 (m)
L2 = 1; % rest length of spring 2 (m)
d = 0.1; % damping coefficient

% Initial conditions
x0 = 0.5; % initial position of the mass
v0 = 0; % initial velocity of the mass

% Simulation time
tspan = [0 10]; % simulation time from 0 to 10 seconds

% Differential equation describing the system dynamics
f = @(t, y) [
    y(2);
    -(k1 * (y(1) - L1) + k2 * (y(1) - L2) + d * y(2)) / m
];

% Solve the differential equations
[t, y] = ode45(f, tspan, [x0; v0]);

% Plot the results
figure;
plot(t, y(:, 1));
title('Mass-Spring System Simulation');
xlabel('Time (s)');
ylabel('Position (m)');
734 chars
31 lines