El uso de las bacterias Gram-Schmidt proceso, incluso si se escribe correctamente, no es una gran idea para este tipo de pregunta, como se verá en general estropear el ganado propiedad de ser una base de vectores propios. Aquí es un enfoque oportunista que trabaja para esta particular de la matriz $A$.
Cosa que es fácil de notar es que el $A+I$ tiene todas las entradas $1$, por lo que no sólo su determinante es $0$, pero su rango es sólo $1$:
$$
A+I=\begin{pmatrix}1\\1\\1\end{pmatrix}\cdot\begin{pmatrix}1&1&1\end{pmatrix}
$$
Así, el espacio propio para $\lambda=-1$ $2$- dimensiones del núcleo de la derecha de la matriz $(1~~1~~1)$, el plano dado por la ecuación de $x+y+z=0$. Desde la traza de $A$ es la suma de sus valores propios (tomada con multiplicidades), y es $0$, el resto de autovalor es $2$, y un autovector de es $(1~~1~~1)^T$, que es ortogonal (como debe ser, dado que $A$ es simétrica) a la otra autoespacio $x+y+z=0$.
Ahora para encontrar $P$ tenemos una base ortogonal del subespacio propio $x+y+z=0$, y, a continuación, la normalización de todos los vectores propios. No puedo ver ningún particularmente agradable ortogonal pares en ese espacio propio, así que me voy a tomar el primer vector que viene a la mente $(1~~-1~~0)^T$ y tomar la cruz de producto con el autovector $(1~~1~~1)^T$$\lambda=2$, para dar a $(-1~~-1~~2)^T$. Dividiendo por las normas de estos vectores puedo encontrar
$$
P=\frac1{\sqrt6}
\begin{pmatrix}\sqrt3&-1&\sqrt2\\-\sqrt3&-1&\sqrt2\\0&2&\sqrt2\end{pmatrix}
\qquad \text{y}\qquad
D=\begin{pmatrix}-1&0&0\\0&-1&0\\0&0&2\end{pmatrix}.
$$