抛体运动，用有限差分法，求求各位啦
对不起，还是我这个小白，我重新发下我的作业（我是学建筑的，希望有人能帮我一把）以下是要求（我自己翻译的日语，有的词语可能翻译的不准确）：
二次元直角坐标系(x, y) 里，由牛顿的运动第二法则，质量m[kg]的质点的运动方程式
m*d^2*x/(dt^2)=fx
m*d^2*y/(dt^2)=fy
在这里(fx, fy)是对质点的作用力[N]的x成分、y 成分。
上式是2阶线性微分方程，可以转化成1阶联系线性微分方程
m*du/dt=fx
dx/dt= u
m*dv/dt= fy
dy/dt= v
这里，(u, v)是质点的速度[m/s]的x成分、y成分。初期条件是时刻t=0的时候，设坐标为(x(0), y(0))，速度(u(0), v(0))。
举例子说明，以水平方向为x轴，垂直方向为y轴，无视空气阻力下，质点只在重力作用下的运动方程式是
m*du/dt= 0
dx/dt= u
m*dv/dt= −mg
dy/dt= v
这里，重力加速度g = 9.8[m/s2]。
用欧拉显式方法（这个不会翻译，英语好像是euler法，中文好像是差分法吧）做分离。时间的分离里添加字幕n
m*(u^(n+1) − u^n)/Δt= 0
(x^(n+1) − x^n)/Δt= u^n
m*(v^(n+1) − v^n)/Δt= −mg
(y^(n+1) − y^n)/Δt= v^n
所以差分式为
u^(n+1) = u^n
x^(n+1) = x^n + Δt*u^n
v^(n+1) = v^n − Δt*g
y^(n+1) = y^n + Δt*v^n
以上是老师的讲解。
课题如下：
時刻t=0 初期条件如下：
質点の初期座標(x0, y0)=(0,0)
質点の初期速度(u0, v0)=(20,20)
重力加速度はg=9.8
时间间隔dt=0.1，T=5.0
得出的数据文件已 13t7032kc46.data 为命名
得出的数据，用质点的投射运动的轨迹（横軸x、縦軸y），制作EXCEL的图像，然后做成视频
哪位能帮我解答下C语言的部分
最近连续通宵画建筑设计图纸，又赶上期末，谁帮帮小白我一把
尽可能在今天中午之前，麻烦大家了
