Question: how to paint four 3D-lines in one plot with points in the 4 matrixs?

I have 4 matrixes such as  P, P1, P2,P3

what can i do paint four 3D-lines  in one plot with points in the 4 matrixs?

 P := Matrix(51, 3);

R := proc (s) options operator, arrow; 0.6e-1*(arctan(10*s/1.3-3.0)/Pi+1/2) end proc;

y := proc (s) options operator, arrow; evalf((-1)*0.6e-1*(arctan(10*s/1.3-3.0)/Pi+1/2)*cos(6*Pi*s/1.3)) end proc;

z := proc (s) options operator, arrow; evalf(0.6e-1*(arctan(10*s/1.3-3.0)/Pi+1/2)*sin(6*Pi*s/1.3)) end proc

for i to 51 do

P[i, 1] := evalf(Int(sqrt(1-((6*Pi/(2.6)*0.6e-1)*R(s))^2-(D(R))(s)^2), s = 0 .. 0.26e-1*i));

P[i, 2] := y(0.26e-1*i);

P[i, 3] := z(0.26e-1*i)

end do

for i to 51 do

P[i, 1] := evalf(Int(sqrt(1-((6*Pi/(2.6)*0.6e-1)*R(s))^2-(D(R))(s)^2),s = 0 .. 0.26e-1*i));

P[i, 2] := y(0.26e-1*i);

P[i, 3] := z(0.26e-1*i)

end do

P1 := Matrix(51, 3):

for i to 51 do

P1[i, 1] := P[i, 1];

P1[i, 2] := P[i, 2]+0.2e-1*P[i, 2]/(P[i, 2]^2+P[i, 3]^2)^(1/2);

P1[i, 3] := P[i, 3]+0.2e-1*P[i, 3]/(P[i, 2]^2+P[i, 3]^2)^(1/2)

end do

P2 := Matrix(51, 3);

for i to 51 do

P2[i, 1] := P[i, 1];

P2[i, 2] := evalf(P[i, 2]+0.2e-1*P[i, 2]*cos(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2)-0.2e-1*P[i, 3]*sin(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2));

P2[i, 3] := evalf(P[i, 3]+0.2e-1*P[i, 2]*sin(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2)+0.2e-1*P[i, 3]*cos(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2))

end do

 P3 := Matrix(51, 3);
 for i to 51 do

P2[i, 1] := P[i, 1];

P2[i, 2] := evalf(P[i, 2]+0.2e-1*P[i, 2]*cos(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2)+0.2e-1*P[i, 3]*sin(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2));

P2[i, 3] := evalf(P[i, 3]-0.2e-1*P[i, 2]*sin(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2)+0.2e-1*P[i, 3]*cos(2*Pi*1/3)/(P[i, 2]^2+P[i, 3]^2)^(1/2))

 end do;
 

 Thank you!

liu