Athanasios,
sorry - obviously I misunderstood your question!
As far as I could understand on the quick the "improved Euler method" I simply would do it like so:
\[Alpha] = 0.9;
\[Beta] = 0.2;
\[Gamma] = 1.2;
f[{x_, y_, z_}] := {z + (y - \[Alpha])*x, 1 - \[Beta]*y - x^2, -x - \[Gamma]*z}
h = 0.1;
impEuler[v_List] := v + h/2 (f[v] + f[v + h f[v]])
data = NestList[impEuler, {1, 2, 3}, 1000];
Graphics3D[{Arrowheads[{.03, #} & /@ Range[0, 1, .05]], Arrow[data]},
BoxRatios -> Automatic, Axes -> True, AxesLabel -> {"X", "Y", "Z"}]

I do hope that this now is more helpful! At least my "solution" above can be used as reference.
Best regards -- Henrik