function Jh=RemoveJacob(J,L);%RemoveJacob Removes the derivatives corresponding to the current  carrying electrodes% Function Jh=RemoveJacob(J,L);% removes all the derivatives that correspond to the% current carrying electrodes.% Change the voltages accordingly, see RemoveVolt.m.%% INPUT% J = Jacobian% L = number of electrodes% OUTPUT% Jh = new Jacobian[rJ,cJ]=size(J);J=reshape(J,L,rJ/L*cJ);kk=1;p=1;for jj=1:cJ for ii=1:L  J(:,p)=[J(kk:L,p);J(1:kk-1,p)]; kk=kk+1;p=p+1;endkk=1;endJ([1 2 L],:)=[];Jh=reshape(J,(L-3)*L,cJ);