
simp=dlmread('simp_test2d.txt');
vert=dlmread('vtx_test2d.txt');
arg_=[1, 100, 110, 75,75, 110,-1, 1]
matS_2=[];
[mat1_]=fantom_InNodes_v1(vert,simp,[],ones(size(vert)), arg_);
matS_2=[matS_2 ECT3D_FEMMatNodes2Simp(vert,simp,mat1_)];
z=20;
z1=100;
for i=1:10
    i
    z=100;
    perm=3;
    z1=z1-10;
    arg_=[1, z, z1, 75, 75, z1, -1, perm]
    %z1 wysokosc walca wewnatrz bazy
    %perm wartosc permittivity
    %z rozmiar walca wewnatrz --> 75+60 dla z=60 
    
    mat1=fantom_InNodes_v1(vert,simp,[],mat1_, arg_);
    matS_2=[matS_2 ECT3D_FEMMatNodes2Simp(vert,simp,mat1)];
%trisurf(simp,vert(:,1),vert(:,2),vert(:,3),matS_2);
   
end

%1; 50; 110; 75;  110;  -1; 1