用仿真程序求解的闭环输出响应y（t）.

用matlab语言编制单变量系统三阶龙格-库塔法求解程序，程序入口要求能接收状态方程各系数阵（A，B，C，D），和输入阶跃函数r（t）=R*1（t）；程序出口应给出输出量y（t）的动态响应数值解序列

单位反馈系统的开环传递函数已知如下

用matlab语句 、函数求取系统闭环零极点，并求取系统闭环状态方程的可控标准型实现。

已知二阶系统状态方程为

写出取计算步长为h时，该系统状态变量X=[12，xx]的四阶龙格-库塔法递推关系式。