オープンキャンパス用 コード

Clear[F, m, a, v, x, x0, g, k, h, t1];
F[x_, v_, t_] := -k*x - h*v;
v[t_] := x'[t]; a[t_] := v'[t];
x[t_] = DSolve[{m a[t] == F[x[t], v[t], t], v[0] == 0, x[0] == x0},x[t], t][[1, 1, 2]];
x[t] // FullSimplify // TraditionalForm // Print;
g = 9.8; x0 = 0.1; k = 500; h = 0.32; t1 = 8;
mass[m0_] = (x[t] /. m -> m0);
Plot[{mass[1], mass[4]}, {t, 0, t1}, PlotRange -> All, 
 PlotLegends -> "Expressions"]