抛体可以是任意角度、任意初速度的！

/*オイラー陽解法を用いて一時反応の濃度を計算するプログラム*/

#include <stdio.h>
#define T 10.0 /*計算領域　T＝10.0の定義*/
#define ni 20 /*格子点公分割　niの定義*/

main()
{
    int i; /*格子点番号　整数型*/
    float t[ni + 1];/*時刻tの一次元配列*/
    float u[ni + 1];/*濃度uの一元元配列*/

    float dt = T / ni; /*時間間隔Δtの定義*/

    FILE *fp; /*ファイルポインタ*/

        printf("%f\n", dt);

        /* 時刻t[i]の計算 */
        for (i = 0; i <= ni; i++ ) {
            t[i] = i*dt;
        }

        /*濃度uの初期値設定*/
        u[0] = 1.0;

        /*濃度u[i]の繰り返し計算*/
        for (i = 1; i <= ni; i++) {
            u[i] = (1.0 - dt)*u[i - 1];
        }

        /*計算結果のファイル出力*/
        fopen_s(&fp, "./c41dt05.data", "w");

        for (i = 0; i <= ni; i++) {
            fprintf(fp, " %-6.3f %e\n", t[i], u[i]);
        }

        fclose(fp);
    return(0);
 这是我们老师讲的课题，除了日语，你能看懂么？

T版，您晚上不睡觉？还是有时差

for (i = 1; i <= ni; i++) {
            u[i] = (1.0 - dt)*u[i - 1];
        }

这一句我还没懂，这算的是什么？你贴错代码了吧，这好像是另一题的解法呀

[此贴子已经被作者于2016-1-7 09:17编辑过]

就算是日语，你也应该知道这代码和题目没啥关系吧

オイラー陽解法を用いて一時反応の濃度を計算するプログラム
油杯陽解法用一时反应浓度计算的程序

百度翻译

那是给你示范微元法，但不是匀变速运动。解法是一样的。

excel能做成轨迹和轨迹视频？还真没用过这种功能！

可以的，就用描點、折缐那個東西。當然直接在C寫也可以。

我的代码：

#include <stdio.h>

#define g 9.8
#define dt 0.1
#define T 5.0

int main(void)
{
    double x = 0, y = 0;
    double vx = 20, vy = 20;
    double t;

    FILE *fp = fopen("13t7032kc46.data", "w");

    printf("time\tx\ty\n");
    fputs("time\tx\ty\n", fp);

    for (t = 0; t <= T; t += dt)
    {
        printf("%.2f\t%.2f\t%.2f\n", t, x, y);
        fprintf(fp, "%.2f\t%.2f\t%.2f\n", t, x, y);

        x += vx * dt;

        y += vy * dt;
        vy -= g * dt;
        if (y < 0) // 如果落到地面了，重置vy的初始速度（忽略阻力）
        {
            y = 0;
            vy = 20;
        }
    }

    fclose(fp);

    return 0;
}