Перевод программы из Матлаба в Математику

s1=tf([0.08],[0.36 0.6 1]);          % непрерывная передаточная функция объекта
T_end=60;                                 % интервал измерений
dt=0.2;                                     %шаг дискретизации
t=0:dt:T_end;                            %массив дискретного времени
N=length(t);                               %размер выборки
u=ones(N,1);                              %моделирование единичного входного воздействия
y=lsim(s1,u,t);                            %моделирование выходного воздействия
n=2;                                           %порядок объекта
R=[y(n:N-1) y(n-1:N-2) u(n:N-1)]; %формирование расширенной матрицы
                                                         %данных
Y=y(n+1:N);                                  % формирование вектора выходных данных
betta=inv(R'*R)*R'*Y;
                                                     % расчет параметров непрерывной модели
T1=dt^2/(1-betta(1)-betta(2));
T2=(betta(2)*T1+T1+dt^2)/dt;
K=betta(3)*T1/dt^2;
 
s2=tf([K],[T1 T2 1]);
y2=lsim(s2,u,t);
plot(t,y,t,y2,':');
grid;