Question: how can i draw odeplot this numeric dsolve

i want to draw plot Y[1] . but i cant. thank for you

i want to take this equation

restart;
with(LinearAlgebra); Digits := 70; with(plots);

K := Matrix(5, 5, [[18278000.0000, 0., 0., 0., 0.], [0., 8478500.0000, 0., 0., 0.], [0., 0., 3874800.0000, 0., 0.], [0., 0., 0., 494710000.0000, 0.], [0., 0., 0., 0., 7707500.0000]]); M := Matrix(5, 5, [[4119700.0000, 0., 0., 0., 0.], [0., 175900.0000, 0., 0., 0.], [0., 0., 52796., 0., 0.], [0., 0., 0., 2002900.0000, 0.], [0., 0., 0., 0., 21711.]]); L := Vector[column](5, [400689.48, 36882.21, 24223.175, 7197.4865, 5007.6466]); V := Vector[column](5, [4.4368, 48.201, 73.39160, 246.999, 355.012]);
NS := ImportMatrix("D:/tabas.txt", datatype = float[8]);
t_NS := NS[() .. (), 1];
acc_NS := NS[() .. (), 2];
plot(NS, t = 0 .. 4.5, size = [800, 400]);
acc := unapply(CurveFitting:-Spline(t_NS, acc_NS/(9.81), t, degree = 1), t);
eq := seq(diff(x[i](t), t$2)+(2*0.5e-1)*sqrt(V[i])*(diff(x[i](t), t))+V[i]*x[i](t) = L[i]/M[i, i], i = 1 .. 5, 1);
ic := [seq({x[i](0) = 0, (D(x[i]))(0) = 0}, i = 1 .. 6, 1)];
res := evalf(seq(dsolve({eq[i], ic[i][]}, numeric, method = rkf45), i = 1 .. 5, 1));
Z := `<,>`(res[1], res[2], res[3], res[4], res[5]);
Y := Phi . Z;
#odeplot(rhs(Y[1]), t = 0 .. 4)

import file text

tabas.txt

Please Wait...