function t,y = euler_stepp(A, B, C, D, dt, T)
    t = 0:T:dt
    n = len(t)
    x = zeros(1, n + 1)
    y = zeros(1, n)

    for i=1:n
        x(0, i) = (A * x(0, i) + B) * dt + x(0, i)
        y(0, i) = C * x(0, i) + D
    end