Ok so I see the pattern...
v_i ^2 on the diagonal and v_ij in the off diagonal entries. Now the problem remains  how do you solve this using matrix techniques (ie in Matlab)? It is now nonlinear so I guess there is no chance of doing this (ie. Linear algebra)?
Dave
