function [mat]=ECT2D_FEMCreate_Inhomogenity(vtx,simp,mat_ref);


rad=[(abs(max(vtx(:,1)))+abs(min(vtx(:,1))))/2,(abs(max(vtx(:,2)))+abs(min(vtx(:,2))))/2];

IleF = input (['Ile walcow:']);
mat(1:size(simp,1),1) = mat_ref;

for ii=1:IleF
    %The Radius of the base of cylinder
    RInh(ii,1) = input (['Outer Radius of the base ([dx,dy]= [' num2str(rad(1)) ',' num2str(rad(2)) '] ): ']);
    RInh(ii,2) = input (['Inner Radius of the base ([dx,dy]= [' num2str(rad(1)) ',' num2str(rad(2)) '] ): ']);
    % Co-ordinates of the centre of the inhomogeneity
    CntInh(ii,1) = input(['X co-ordinate (max=' num2str(max(vtx(:,1))) ', min=' num2str(min(vtx(:,1))) '): ']);
    CntInh(ii,2) = input(['Y co-ordinate (max=' num2str(max(vtx(:,2))) ', min=' num2str(min(vtx(:,2))) '): ']);
    % Conductivity of the inhomogenity
    CndInh(ii,1) = input(['Conductivity within current inhomogeneity : ']);
    % Initiate tetrahedra's counter
    
    

    for i=1:size(simp,1)
        x_c = (vtx(simp(i,1),1)+vtx(simp(i,2),1)+vtx(simp(i,3),1))/3;
        y_c = (vtx(simp(i,1),2)+vtx(simp(i,2),2)+vtx(simp(i,3),2))/3;
        x = [x_c y_c];
        d = sqrt(sum((x(1:2)-CntInh(ii,1:2)).^2));
         if (d <=RInh(ii,1) && d >=RInh(ii,2))
             mat(i,1)=CndInh(ii,1);
         end
    end
end

