function [t,y] = euler_impulse(A,B,C,D,dt,T)
    t = 0:dt:T;
    n = length(t);
    x = zeros(1,n+1);
    y = zeros(1,n);
    u = 1/dt;
    for i=1:n
        x(1,i+1)=(A*x(1,i)+B*u)*dt+x(1,i);
        y(1,i)=C*x(1,i)+D*u;
    end
[t,y] = euler_impulse(-1/7, 4/7, 1, 0, 1,  60);
plot(y,t);