% COMSOL Multiphysics Model M-file% Generated by COMSOL 3.5 snapshot (COMSOL 3.5.0.468, $Date: 2008/08/15 16:04:11 $)flclear fem% COMSOL versionclear vrsnvrsn.name = 'COMSOL 3.5';vrsn.ext = ' snapshot';vrsn.major = 0;vrsn.build = 468;vrsn.rcs = '$Name:  $';vrsn.date = '$Date: 2008/08/15 16:04:11 $';fem.version = vrsn;% Geometryg1=rect2('2.4','0.2','base','corner','pos',{'-1.2','1'},'rot','0');g2=rect2('2.4','0.2','base','corner','pos',{'-1.2','-1.2'},'rot','0');g3=rect2('0.2','2.4','base','corner','pos',{'-1.2','-1.2'},'rot','0');g4=rect2('0.2','2.4','base','corner','pos',{'1','-1.2'},'rot','0');g5=ellip2('0.1','0.3','base','center','pos',{'0','0'},'rot','0');g6=ellip2('0.3','0.1','base','center','pos',{'0','0'},'rot','0');g7=geomcomp({g5,g6},'ns',{'g5','g6'},'sf','g5+g6','edge','none');g8=geomdel(g7);g9=rect2('2','2','base','corner','pos',{'-1','-1'},'rot','0');% Analyzed geometryclear ss.objs={g1,g2,g3,g4,g8,g9};s.name={'R1','R2','R3','R4','CO2','R5'};s.tags={'g1','g2','g3','g4','g8','g9'};fem.draw=struct('s',s);fem.geom=geomcsg(fem);% Initialize meshfem.mesh=meshinit(fem, ...                  'hmax',[0.05]);% (Default values are not included)% Application mode 1clear applappl.mode.class = 'InPlaneWaves';appl.name = 'rfweh';appl.module = 'RF';appl.sshape = 2;appl.assignsuffix = '_rfweh';clear propprop.field='TES';prop.inputvar='lambda';appl.prop = prop;clear bndbnd.type = {'SC','cont'};bnd.ind = [1,1,1,2,1,2,1,2,1,2,2,2,2,1,2,1,2,2,2,2,1,1,1,1,2,2,2,2,2,2, ...  2,2];appl.bnd = bnd;clear equequ.coordOn = {{1;1},{1;0},{0;1},{0;0},{0;0}};equ.epsilonr = {1,1,1,1,3};equ.Stype = {'coord','coord','coord','none','none'};equ.ind = [1,2,1,3,4,3,1,2,1,5];appl.equ = equ;appl.var = {'lambda0','0.5'};fem.appl{1} = appl;fem.frame = {'ref'};fem.border = 1;clear units;units.basesystem = 'SI';fem.units = units;% ODE Settingsclear odeclear units;units.basesystem = 'SI';ode.units = units;fem.ode=ode;% Multiphysicsfem=multiphysics(fem);% Extend meshfem.xmesh=meshextend(fem);% Solve problemfem.sol=femstatic(fem, ...                  'complexfun','on', ...                  'solcomp',{'scEz'}, ...                  'outcomp',{'scEz'}, ...                  'blocksize','auto');% Save current fem structure for restart purposesfem0=fem;% Plot solutionpostplot(fem, ...         'tridata',{'scEz','cont','internal','unit','V/m'}, ...         'trimap','jet(1024)', ...         'title','Surface: Scattered electric field, z component [V/m]', ...         'axis',[-1.471475468307245,1.471475468307245,-1.3200000524520874,1.3200000524520874]);% Plot solutionpostplot(fem, ...         'tridata',{'scEz','cont','internal','unit','V/m'}, ...         'trimap','jet(1024)', ...         'title','Surface: Scattered electric field, z component [V/m]', ...         'axis',[-1.333770544802556,1.333770544802556,-1.3200000524520874,1.3200000524520874]);