got results of the first global newton iteration

updates: 
1) global residual vector and jacobian matrix calculation was confirmed and compared to matlab results based on their norm values (for the first newton's iteration).
2) direct solver is used in vprime function.