% Driver script for solving the 2D Maxwell's equations on TM form Globals2D; % Polynomial order used for approximation N = 8; % Initial fields load values1.dat; xa1 = values1(:,1); ya1 = values1(:,2); hxa1 = values1(:,3); hya1 = values1(:,4); eza1 = values1(:,5); %--------------- load values3.dat; xa3 = values3(:,1); ya3 = values3(:,2); hxa3 = values3(:,3); hya3 = values3(:,4); eza3 = values3(:,5); load values4.dat; xa4 = values4(:,1); ya4 = values4(:,2); hxa4 = values4(:,3); hya4 = values4(:,4); eza4 = values4(:,5); load values5.dat; xa5 = values5(:,1); ya5 = values5(:,2); hxa5 = values5(:,3); hya5 = values5(:,4); eza5 = values5(:,5); %--------------- load values7.dat; xa7 = values7(:,1); ya7 = values7(:,2); hxa7 = values7(:,3); hya7 = values7(:,4); eza7 = values7(:,5); load values8.dat; xa8 = values8(:,1); ya8 = values8(:,2); hxa8 = values8(:,3); hya8 = values8(:,4); eza8 = values8(:,5); load values9.dat; xa9 = values9(:,1); ya9 = values9(:,2); hxa9 = values9(:,3); hya9 = values9(:,4); eza9 = values9(:,5); load values10.dat; xa10 = values10(:,1); ya10 = values10(:,2); hxa10 = values10(:,3); hya10 = values10(:,4); eza10 = values10(:,5); load values11.dat; xa11 = values11(:,1); ya11 = values11(:,2); hxa11 = values11(:,3); hya11 = values11(:,4); eza11 = values11(:,5); load values12.dat; xa12 = values12(:,1); ya12 = values12(:,2); hxa12 = values12(:,3); hya12 = values12(:,4); eza12 = values12(:,5); load values13.dat; xa13 = values13(:,1); ya13 = values13(:,2); hxa13 = values13(:,3); hya13 = values13(:,4); eza13 = values13(:,5); load values14.dat; xa14 = values14(:,1); ya14 = values14(:,2); hxa14 = values14(:,3); hya14 = values14(:,4); eza14 = values14(:,5); load values15.dat; xa15 = values15(:,1); ya15 = values15(:,2); hxa15 = values15(:,3); hya15 = values15(:,4); eza15 = values15(:,5); load values16.dat; xa16 = values16(:,1); ya16 = values16(:,2); hxa16 = values16(:,3); hya16 = values16(:,4); eza16 = values16(:,5); load values17.dat; xa17 = values17(:,1); ya17 = values17(:,2); hxa17 = values17(:,3); hya17 = values17(:,4); eza17 = values17(:,5); load values18.dat; xa18 = values18(:,1); ya18 = values18(:,2); hxa18 = values18(:,3); hya18 = values18(:,4); eza18 = values18(:,5); load values19.dat; xa19 = values19(:,1); ya19 = values19(:,2); hxa19 = values19(:,3); hya19 = values19(:,4); eza19 = values19(:,5); load values20.dat; xa20 = values20(:,1); ya20 = values20(:,2); hxa20 = values20(:,3); hya20 = values20(:,4); eza20 = values20(:,5); load values21.dat; xa21 = values21(:,1); ya21 = values21(:,2); hxa21 = values21(:,3); hya21 = values21(:,4); eza21 = values21(:,5); load values22.dat; xa22 = values22(:,1); ya22 = values22(:,2); hxa22 = values22(:,3); hya22 = values22(:,4); eza22 = values22(:,5); load values23.dat; xa23 = values23(:,1); ya23 = values23(:,2); hxa23 = values23(:,3); hya23 = values23(:,4); eza23 = values23(:,5); load values24.dat; xa24 = values24(:,1); ya24 = values24(:,2); hxa24 = values24(:,3); hya24 = values24(:,4); eza24 = values24(:,5); load values25.dat; xa25 = values25(:,1); ya25 = values25(:,2); hxa25 = values25(:,3); hya25 = values25(:,4); eza25 = values25(:,5); load values26.dat; xa26 = values26(:,1); ya26 = values26(:,2); hxa26 = values26(:,3); hya26 = values26(:,4); eza26 = values26(:,5); load values27.dat; xa27 = values27(:,1); ya27 = values27(:,2); hxa27 = values27(:,3); hya27 = values27(:,4); eza27 = values27(:,5); load values28.dat; xa28 = values28(:,1); ya28 = values28(:,2); hxa28 = values28(:,3); hya28 = values28(:,4); eza28 = values28(:,5); load values29.dat; xa29 = values29(:,1); ya29 = values29(:,2); hxa29 = values29(:,3); hya29 = values29(:,4); eza29 = values29(:,5); load values30.dat; xa30 = values30(:,1); ya30 = values30(:,2); hxa30 = values30(:,3); hya30 = values30(:,4); eza30 = values30(:,5); load values31.dat; xa31 = values31(:,1); ya31 = values31(:,2); hxa31 = values31(:,3); hya31 = values31(:,4); eza31 = values31(:,5); load values32.dat; xa32 = values32(:,1); ya32 = values32(:,2); hxa32 = values32(:,3); hya32 = values32(:,4); eza32 = values32(:,5); load values33.dat; xa33 = values33(:,1); ya33 = values33(:,2); hxa33 = values33(:,3); hya33 = values33(:,4); eza33 = values33(:,5); load values34.dat; xa34 = values34(:,1); ya34 = values34(:,2); hxa34 = values34(:,3); hya34 = values34(:,4); eza34 = values34(:,5); load values35.dat; xa35 = values35(:,1); ya35 = values35(:,2); hxa35 = values35(:,3); hya35 = values35(:,4); eza35 = values35(:,5); load values36.dat; xa36 = values36(:,1); ya36 = values36(:,2); hxa36 = values36(:,3); hya36 = values36(:,4); eza36 = values36(:,5); load values37.dat; xa37 = values37(:,1); ya37 = values37(:,2); hxa37 = values37(:,3); hya37 = values37(:,4); eza37 = values37(:,5); load values38.dat; xa38 = values38(:,1); ya38 = values38(:,2); hxa38 = values38(:,3); hya38 = values38(:,4); eza38 = values38(:,5); load values39.dat; xa39 = values39(:,1); ya39 = values39(:,2); hxa39 = values39(:,3); hya39 = values39(:,4); eza39 = values39(:,5); load values40.dat; xa40 = values40(:,1); ya40 = values40(:,2); hxa40 = values40(:,3); hya40 = values40(:,4); eza40 = values40(:,5); load values41.dat; xa41 = values41(:,1); ya41 = values41(:,2); hxa41 = values41(:,3); hya41 = values41(:,4); eza41 = values41(:,5); load values42.dat; xa42 = values42(:,1); ya42 = values42(:,2); hxa42 = values42(:,3); hya42 = values42(:,4); eza42 = values42(:,5); load values43.dat; xa43 = values43(:,1); ya43 = values43(:,2); hxa43 = values43(:,3); hya43 = values43(:,4); eza43 = values43(:,5); load values44.dat; xa44 = values44(:,1); ya44 = values44(:,2); hxa44 = values44(:,3); hya44 = values44(:,4); eza44 = values44(:,5); load values45.dat; xa45 = values45(:,1); ya45 = values45(:,2); hxa45 = values45(:,3); hya45 = values45(:,4); eza45 = values45(:,5); %---------------- load values50.dat; xa50 = values50(:,1); ya50 = values50(:,2); hxa50 = values50(:,3); hya50 = values50(:,4); eza50 = values50(:,5); % Read in Mesh [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse1.dat'); % Initialize solver and construct grid and metric. StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx1,Hy1,Ez1] = magH2D(x,y,xa1,ya1,hxa1,hya1,eza1); PlotField2D(N, x, y, Hx1); pause(1);close; PlotField2D(N, x, y, Hy1); pause(1);close; K1 = K; % Read in Mesh [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse2.dat'); % Initialize solver and construct grid and metric. StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx2,Hy2,Ez2] = magH2D(x,y,xa3,ya3,hxa3,hya3,eza3); PlotField2D(N, x, y, Hx2); pause(1);close; PlotField2D(N, x, y, Hy2); pause(1);close; K2 = K; % Read in Mesh [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse3.dat'); % Initialize solver and construct grid and metric. StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx3,Hy3,Ez3] = magH2D(x,y,xa4,ya4,hxa4,hya4,eza4); PlotField2D(N, x, y, Hx3); pause(1);close; PlotField2D(N, x, y, Hy3); pause(1);close; K3 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse4.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx4,Hy4,Ez4] = magH2D(x,y,xa5,ya5,hxa5,hya5,eza5); K4 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse5.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx5,Hy5,Ez5] = magH2D(x,y,xa7,ya7,hxa7,hya7,eza7); K5 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse6.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx6,Hy6,Ez6] = magH2D(x,y,xa8,ya8,hxa8,hya8,eza8); K6 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse7.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx7,Hy7,Ez7] = magH2D(x,y,xa9,ya9,hxa9,hya9,eza9); K7 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse8.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx8,Hy8,Ez8] = magH2D(x,y,xa10,ya10,hxa10,hya10,eza10); K8 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse9.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx9,Hy9,Ez9] = magH2D(x,y,xa11,ya11,hxa11,hya11,eza11); K9 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse10.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx10,Hy10,Ez10] = magH2D(x,y,xa12,ya12,hxa12,hya12,eza12); K10 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse11.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx11,Hy11,Ez11] = magH2D(x,y,xa13,ya13,hxa13,hya13,eza13); K11 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse12.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx12,Hy12,Ez12] = magH2D(x,y,xa14,ya14,hxa14,hya14,eza14); K12 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse13.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx13,Hy13,Ez13] = magH2D(x,y,xa15,ya15,hxa15,hya15,eza15); K13 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse14.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx14,Hy14,Ez14] = magH2D(x,y,xa16,ya16,hxa16,hya16,eza16); K14 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse15.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx15,Hy15,Ez15] = magH2D(x,y,xa17,ya17,hxa17,hya17,eza17); K15 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse16.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx16,Hy16,Ez16] = magH2D(x,y,xa18,ya18,hxa18,hya18,eza18); K16 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse17.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx17,Hy17,Ez17] = magH2D(x,y,xa19,ya19,hxa19,hya19,eza19); K17 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse18.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx18,Hy18,Ez18] = magH2D(x,y,xa20,ya20,hxa20,hya20,eza20); K18 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse19.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx19,Hy19,Ez19] = magH2D(x,y,xa21,ya21,hxa21,hya21,eza21); K19 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse20.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx20,Hy20,Ez20] = magH2D(x,y,xa22,ya22,hxa22,hya22,eza22); K20 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse21.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx21,Hy21,Ez21] = magH2D(x,y,xa23,ya23,hxa23,hya23,eza23); K21 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse22.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx22,Hy22,Ez22] = magH2D(x,y,xa24,ya24,hxa24,hya24,eza24); K22 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse23.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx23,Hy23,Ez23] = magH2D(x,y,xa25,ya25,hxa25,hya25,eza25); K23 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse24.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx24,Hy24,Ez24] = magH2D(x,y,xa26,ya26,hxa26,hya26,eza26); K24 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse25.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx25,Hy25,Ez25] = magH2D(x,y,xa27,ya27,hxa27,hya27,eza27); K25 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse26.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx26,Hy26,Ez26] = magH2D(x,y,xa28,ya28,hxa28,hya28,eza28); K26 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse27.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx27,Hy27,Ez27] = magH2D(x,y,xa29,ya29,hxa29,hya29,eza29); K27 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse28.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx28,Hy28,Ez28] = magH2D(x,y,xa30,ya30,hxa30,hya30,eza30); K28 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse29.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx29,Hy29,Ez29] = magH2D(x,y,xa31,ya31,hxa31,hya31,eza31); K29 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse30.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx30,Hy30,Ez30] = magH2D(x,y,xa32,ya32,hxa32,hya32,eza32); K30 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse31.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx31,Hy31,Ez31] = magH2D(x,y,xa33,ya33,hxa33,hya33,eza33); K31 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse32.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx32,Hy32,Ez32] = magH2D(x,y,xa34,ya34,hxa34,hya34,eza34); K32 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse33.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx33,Hy33,Ez33] = magH2D(x,y,xa35,ya35,hxa35,hya35,eza35); K33 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse34.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx34,Hy34,Ez34] = magH2D(x,y,xa36,ya36,hxa36,hya36,eza36); K34 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse35.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx35,Hy35,Ez35] = magH2D(x,y,xa37,ya37,hxa37,hya37,eza37); K35 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse36.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx36,Hy36,Ez36] = magH2D(x,y,xa38,ya38,hxa38,hya38,eza38); K36 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse37.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx37,Hy37,Ez37] = magH2D(x,y,xa39,ya39,hxa39,hya39,eza39); K37 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse38.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx38,Hy38,Ez38] = magH2D(x,y,xa40,ya40,hxa40,hya40,eza40); K38 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse39.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx39,Hy39,Ez39] = magH2D(x,y,xa41,ya41,hxa41,hya41,eza41); K39 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse40.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx40,Hy40,Ez40] = magH2D(x,y,xa42,ya42,hxa42,hya42,eza42); K40 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse41.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx41,Hy41,Ez41] = magH2D(x,y,xa43,ya43,hxa43,hya43,eza43); K41 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse42.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx42,Hy42,Ez42] = magH2D(x,y,xa44,ya44,hxa44,hya44,eza44); K42 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse43.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx43,Hy43,Ez43] = magH2D(x,y,xa45,ya45,hxa45,hya45,eza45); K43 = K; [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarse44.dat'); StartUp2Dmichal; PlotMesh2D; pause(1); close; [Hx44,Hy44,Ez44] = magH2D(x,y,xa50,ya50,hxa50,hya50,eza50); K44 = K; % Globals2D; % Read in Mesh [Nv, Nn, VX, VY, K, EToV, mux, epx, sgx] = MeshReaderAnsys2D('EGcoarseAll.dat'); % Initialize solver and construct grid and metric. StartUp2D; PlotMesh2D; pause(10); % Set materials properties Epsilon = ones(Np,1)*epx; Mi = ones(Np,1)*mux; Sigma = ones(Np,1)*sgx; % Ez = zeros(Np, K); % ? 14,16,19,21 | 15,17,18,20 Hx = zeros(Np, K); Hy = zeros(Np, K); for i = 1:K for j = 1:Np if i <= K1 Hx(j,i) = Hx1(j,i); Hy(j,i) = Hy1(j,i); end; if i > K1 && i <= (K1 + K2) Hx(j,i) = Hx2(j,i-K1); Hy(j,i) = Hy2(j,i-K1); end; if i > (K1 + K2) && i <= (K2 + K3) Hx(j,i) = Hx3(j,i-K1-K2); Hy(j,i) = Hy3(j,i-K1-K2); end; if i > (K2 + K3) && i <= (K3 + K4) Hx(j,i) = Hx4(j,i-K2-K3); Hy(j,i) = Hy4(j,i-K2-K3); end; if i > (K3 + K4) && i <= (K4 + K5) Hx(j,i) = Hx5(j,i-K3-K4); Hy(j,i) = Hy5(j,i-K3-K4); end; if i > (K4 + K5) && i <= (K5 + K6) Hx(j,i) = Hx6(j,i-K4-K5); Hy(j,i) = Hy6(j,i-K4-K5); end; if i > (K5 + K6) && i <= (K6 + K7) Hx(j,i) = Hx7(j,i-K5-K6); Hy(j,i) = Hy7(j,i-K5-K6); end; if i > (K6 + K7) && i <= (K7 + K8) Hx(j,i) = Hx8(j,i-K6-K7); Hy(j,i) = Hy8(j,i-K6-K7); end; if i > (K7 + K8) && i <= (K8 + K9) Hx(j,i) = Hx9(j,i-K7-K8); Hy(j,i) = Hy9(j,i-K7-K8); end; if i > (K8 + K9) && i <= (K9 + K10) Hx(j,i) = Hx10(j,i-K8-K9); Hy(j,i) = Hy10(j,i-K8-K9); end; if i > (K9 + K10) && i <= (K10 + K11) Hx(j,i) = Hx11(j,i-K9-K10); Hy(j,i) = Hy11(j,i-K9-K10); end; if i > (K10 + K11) && i <= (K11 + K12) Hx(j,i) = Hx12(j,i-K10-K11); Hy(j,i) = Hy12(j,i-K10-K11); end; if i > (K11 + K12) && i <= (K12 + K13) Hx(j,i) = Hx13(j,i-K11-K12); Hy(j,i) = Hy13(j,i-K11-K12); end; if i > (K12 + K13) && i <= (K13 + K14) Hx(j,i) = Hx14(j,i-K12-K13); Hy(j,i) = Hy14(j,i-K12-K13); Ez(j,i) = 1; end; if i > (K13 + K14) && i <= (K14 + K15) Hx(j,i) = Hx15(j,i-K13-K14); Hy(j,i) = Hy15(j,i-K13-K14); Ez(j,i) = -1; end; if i > (K14 + K15) && i <= (K15 + K16) Hx(j,i) = Hx16(j,i-K14-K15); Hy(j,i) = Hy16(j,i-K14-K15); Ez(j,i) = 1; end; if i > (K15 + K16) && i <= (K16 + K17) Hx(j,i) = Hx17(j,i-K15-K16); Hy(j,i) = Hy17(j,i-K15-K16); Ez(j,i) = -1; end; if i > (K16 + K17) && i <= (K17 + K18) Hx(j,i) = Hx18(j,i-K16-K17); Hy(j,i) = Hy18(j,i-K16-K17); Ez(j,i) = -1; end; if i > (K17 + K18) && i <= (K18 + K19) Hx(j,i) = Hx19(j,i-K17-K18); Hy(j,i) = Hy19(j,i-K17-K18); Ez(j,i) = 1; end; if i > (K18 + K19) && i <= (K19 + K20) Hx(j,i) = Hx20(j,i-K18-K19); Hy(j,i) = Hy20(j,i-K18-K19); Ez(j,i) = -1; end; if i > (K19 + K20) && i <= (K20 + K21) Hx(j,i) = Hx21(j,i-K19-K20); Hy(j,i) = Hy21(j,i-K19-K20); Ez(j,i) = 1; end; if i > (K20 + K21) && i <= (K21 + K22) Hx(j,i) = Hx22(j,i-K20-K21); Hy(j,i) = Hy22(j,i-K20-K21); end; if i > (K21 + K22) && i <= (K22 + K23) Hx(j,i) = Hx23(j,i-K21-K22); Hy(j,i) = Hy23(j,i-K21-K22); end; if i > (K22 + K23) && i <= (K23 + K24) Hx(j,i) = Hx24(j,i-K22-K23); Hy(j,i) = Hy24(j,i-K22-K23); end; if i > (K23 + K24) && i <= (K24 + K25) Hx(j,i) = Hx25(j,i-K23-K24); Hy(j,i) = Hy25(j,i-K23-K24); end; if i > (K24 + K25) && i <= (K25 + K26) Hx(j,i) = Hx26(j,i-K24-K25); Hy(j,i) = Hy26(j,i-K24-K25); end; if i > (K25 + K26) && i <= (K26 + K27) Hx(j,i) = Hx27(j,i-K25-K26); Hy(j,i) = Hy27(j,i-K25-K26); end; if i > (K26 + K27) && i <= (K27 + K28) Hx(j,i) = Hx28(j,i-K26-K27); Hy(j,i) = Hy28(j,i-K26-K27); end; if i > (K27 + K28) && i <= (K28 + K29) Hx(j,i) = Hx29(j,i-K27-K28); Hy(j,i) = Hy29(j,i-K27-K28); end; if i > (K28 + K29) && i <= (K29 + K30) Hx(j,i) = Hx30(j,i-K28-K29); Hy(j,i) = Hy30(j,i-K28-K29); end; if i > (K29 + K30) && i <= (K30 + K31) Hx(j,i) = Hx31(j,i-K29-K30); Hy(j,i) = Hy31(j,i-K29-K30); end; if i > (K30 + K31) && i <= (K31 + K32) Hx(j,i) = Hx32(j,i-K30-K31); Hy(j,i) = Hy32(j,i-K30-K31); end; if i > (K31 + K32) && i <= (K32 + K33) Hx(j,i) = Hx33(j,i-K31-K32); Hy(j,i) = Hy33(j,i-K31-K32); end; if i > (K32 + K33) && i <= (K33 + K34) Hx(j,i) = Hx34(j,i-K32-K33); Hy(j,i) = Hy34(j,i-K32-K33); end; if i > (K33 + K34) && i <= (K34 + K35) Hx(j,i) = Hx35(j,i-K33-K34); Hy(j,i) = Hy35(j,i-K33-K34); end; if i > (K34 + K35) && i <= (K35 + K36) Hx(j,i) = Hx36(j,i-K34-K35); Hy(j,i) = Hy36(j,i-K34-K35); end; if i > (K35 + K36) && i <= (K36 + K37) Hx(j,i) = Hx37(j,i-K35-K36); Hy(j,i) = Hy37(j,i-K35-K36); end; if i > (K36 + K37) && i <= (K37 + K38) Hx(j,i) = Hx38(j,i-K36-K37); Hy(j,i) = Hy38(j,i-K36-K37); end; if i > (K37 + K38) && i <= (K38 + K39) Hx(j,i) = Hx39(j,i-K37-K38); Hy(j,i) = Hy39(j,i-K37-K38); end; if i > (K38 + K39) && i <= (K39 + K40) Hx(j,i) = Hx40(j,i-K38-K39); Hy(j,i) = Hy40(j,i-K38-K39); end; if i > (K39 + K40) && i <= (K40 + K41) Hx(j,i) = Hx41(j,i-K39-K40); Hy(j,i) = Hy41(j,i-K39-K40); end; if i > (K40 + K41) && i <= (K41 + K42) Hx(j,i) = Hx42(j,i-K40-K41); Hy(j,i) = Hy42(j,i-K40-K41); end; if i > (K41 + K42) && i <= (K42 + K43) Hx(j,i) = Hx43(j,i-K41-K42); Hy(j,i) = Hy43(j,i-K41-K42); end; if i > (K42 + K43) && i <= (K43 + K44) Hx(j,i) = Hx44(j,i-K42-K43); Hy(j,i) = Hy44(j,i-K42-K43); end; % end; end; % PlotMesh2D; pause(5); close; PlotField2D(N, x, y, Hx); pause(3); PlotField2D(N, x, y, Hy); pause(3); PlotField2D(N, x, y, Ez); pause(3); PlotField2D(N, x, y, Epsilon); pause(1); close; PlotField2D(N, x, y, Mi); pause(1); close; PlotField2D(N, x, y, Sigma); pause(1); close; % Solve Problem %[Hx,Hy,Ez,time] = Maxwell2D(Hx,Hy,Ez,FinalTime,Mi,Epsilon,Sigma);