%[V1]=subV1(Y1,Im1,Ip1,E); function [V1]=subV1(Y1,Im1,Ip1,E) I1=[-Im1;-Ip1;0;0;0;0]; V1=Y1\I1; V1=V1+E;