map3.m

This software gave volume MoM solution by the CG-FFT method

p(1)= -4.680597730198440e-009;
p(2)= -0.00396491390056;
p(3)= 9.505870907267184e+005;
p1(1)=  2.338081569120618e-009;
p1(2)= 0.00199999792780;
p1(3)= 2.748207907009300e+005;
p2(1)=  2.575016187701492e-010;
p2(2)= -0.00262750109238;
p2(3)= 8.346181327244387e+002;
break
for j=1:110,
   disp(j);
   s=['load h36',num2str(3000+(j-1)*10),'.mat'];
   eval(s);
   aharon36=aharon;
   s=['load h37',num2str(3000+(j-1)*10),'.mat'];
   eval(s);
   aharon37=aharon;
   s=['load h37',num2str(3000+(j-2)*10),'.mat'];
   eval(s);
   aharon371=aharon;
   in=abs(round(polyval(p2,(3000+(j-1)*10)*1000)/100));
   aharon37=[aharon371(101-in:101,:) ; aharon37(2:101-in-1,:)];
   clear aharon;
   clear aharon371;
   
   sizex=round((polyval(p,(3000+(j-1)*10)*1000)-160000)/100)+1;
   sizex=sizex+round((-polyval(p1,(3000+(j-1)*10)*1000-in*100)+860000)/100)+1;
%   aharon1=-32000*ones(101,sizex);
   for j1=1:100,
      sizex1=round((polyval(p,(3000+(j-1)*10)*1000)-160000)/100)+1;
      sizex=sizex1+round((-polyval(p1,(3000+(j-1)*10)*1000-in*100)+860000)/100)+1;
      aharon1(j1,1:sizex)=[aharon36(j1,1:sizex1-1)  aharon37(j1,7001-(sizex-sizex1):7001)];
   end;
   s=['save h367',num2str(3000+(j-1)*10),'.mat aharon1'];
   eval(s);
end;