run a simulation in matlab

main.m
% For example, to simulate a simple harmonic oscillator:
m = 1;      % mass
k = 5;      % spring constant
x0 = 1;     % initial position
v0 = 0;     % initial velocity

% Equation of motion: m * x'' + k * x = 0
f = @(t, y) [y(2); -k/m * y(1)];

tspan = [0, 10];    % time interval to simulate
y0 = [x0; v0];      % initial conditions
334 chars
12 lines

main.m
% Solve the differential equation using ode45
[t, y] = ode45(f, tspan, y0);

% Extract position and velocity
x = y(:, 1);
v = y(:, 2);
135 chars
7 lines

main.m
% Plot position and velocity
figure;
subplot(2, 1, 1);
plot(t, x);
ylabel('Position');

subplot(2, 1, 2);
plot(t, v)
xlabel('Time');
ylabel('Velocity');
153 chars
11 lines