高手帮帮忙

hhj

有一系统，x(t)是基础位移输入，y(t)是位移系统响应，
系统描述为
my''+cy'+ky'=cx'+kx'       (式中y''表示对y去两次导;y'表示y对求一次导; x'对x求一次导)
式中，m、c、k 分别是系统质量、阻尼和刚度。假设m=1 kg，c=50
N·s/m，k=3600 N/m。
(1) 写出系统的传递函数，
(2) 用Matlab 计算系统的频响函数，作出频响函数的幅值和相位图；
(3) 用Matlab求出系统的零点和极点，作出零点和极点的分布图，判断系统稳定性。

这个题的程序我编了好久老出错，请论坛的朋友们帮帮忙做做。
谢谢