represent the following differential equation as a function: d2y(t)/dt2 + 1000*dy(t)/dt + 1000000*y(t) = 1000000*x(t) in matlab assuming x(t) is a step function using the math toolkit in matlab in matlab

main.m
function dy = myODE(t, y)
    dy = zeros(2,1);
    dy(1) = y(2);
    dy(2) = 1000000*x(t) - 1000000*y(1) - 1000*y(2);
end

% Define x(t) as a step function
function x_val = x(t)
    x_val = 1 * (t >= 0); % Step function for t >= 0, x(t) = 1
end

% Define the time vector
tspan = [0 5]; % From t=0 to t=5

% Solve the differential equation
[t, y] = ode45(@myODE, tspan, [0 0]);
377 chars
17 lines