История изменений

int[n_] := Evaluate[Assuming[n >= 0,
    Integrate[Exp[-x^2] x^n,
     {x, 0, Infinity}]]];
pint[p_] :=
  Total[MapIndexed[#*int[First[#2] - 1] &,
    CoefficientList[p, x]]];
f[i_] := N[HermiteH[i, x], 400];
F[n_] := Module[{i, v},
   v = Table[0, {i, n}];
   v[[1]] = f[1];
   For[i = 2, i <= n, i++,
    v[[i]] = Expand[f[i] -
       Sum[
        v[[j]] pint[f[i]*v[[j]]]/pint[v[[j]]*v[[j]]],
        {j, 1, i - 1}]]];
   v];

psi = Integrate[(F[20] /. x -> t), {t, 0, x}];
Table[
  Table[pint[(x^2 + 2*x + 2)*psi[[i]]*psi[[j]]], {i, 1, 
    Length[psi]}],
  {j, 1, Length[psi]}];
N[%]