*********************************
As you say: grad(f) = <-2xy, -x^2, 3z^2> .
This means grad(f) (5, -1, 1)=<-2(5)(-1), -(5^2), 3(1^2)> = <10, -25, 3> unlike what you say. It looks like you made the mistake of substituting -1 for x in x^2, but really you needed to substitute x=5 into -x^2.
Then since v=<-3,3,4>, v^ = <-3,3,4>/sqrt(34), as you say, so
D_v^ f(5, -1, 1)=<10, -25,3>.<-3,3,4>/sqrt(34) =(-30-75+12)/sqrt(34)=-93/sqrt(34)