function [A]=UpdateFemMatrix(Agrad,Kb,M,S,sigma);%UpdateFemMatrix Assembles the system matrix for EIT% Function [A]=UpdateFemMatrix(Agrad,Kb,M,S,sigma);% updates the system matrix A given a new conductivity vector sigma.%% INPUT% Agrad = the gradient part of the system matrix (see FemMatrix.m)% Kb,M and S = other blocks of the system matrix (see FemMatrix.m)% sigma = conductivity (or admittivity) vector%% OUTPUT% A = the updated system matrix rAgrad=sqrt(size(Agrad,1));A=[reshape(Agrad*sigma,rAgrad,rAgrad)+Kb,M;M.',S];