a = 0.2; b = 0.2; c = 9.25; del = 0.01;
xe = (c - Sqrt[c^2 - 4*a*b])/2;
ye = -(c - Sqrt[c^2 - 4*a*b])/(2*a);
ze = (c - Sqrt[c^2 - 4*a*b])/(2*a);
jac = {{0, -1, -1}, {1, a, 0}, {ze, 0, xe - c}};
vec = Eigenvectors[jac];
gramschmidt[w_?MatrixQ] := Module[{v = ConstantArray[0, Length[w]]},
Table[
v[[n]] = w[[n]] -
Sum[(v[[i]].w[[n]]/v[[i]].v[[i]])*v[[i]], {i, 1, n - 1}],
{n, 1, Length[w]}];
v
]
gramschmidt[vec]