[求助]如何用标准C编写系统阶跃响应输出图形？

各位大侠们，请多多帮忙：一二阶线性系统16/(s^2+2.4s+16)，用标准C语言完成：给它个阶跃输入，编程完成输出，并且画出图。怎么设计啊，请各位指点。
注：全用标准C实现。
先谢谢各位了！！

为什么不用matlab？？？
用c我想也行，先把系统函数变换到頻域，和阶跃函数的傅里叶变换相乘，就是它的頻域输出，就可以画出输出频谱图，也可以变换到时域画出时域图。过程可能会有些麻烦，因为标准c函数库中没有相关的函数，要自己写。可以试试，会是一道很好的练习题