很急，关于fortran的编程问题

题目：用四阶Runge-kutta方法，求解阻尼振动方程md^2x/dt^2+cdx/dt+kx=0,已知质量m=100,劲度系数k=100,阻尼系数c=20,初速v=0，初始位置x=100,数值计算物体的位置和速度。（用fortran语言)还有画图！
有急用，谢谢大家了！！！