wglmxy3cf_k.m

集中质量法求解二维弹性波动散射问题的函数

%wglm function calculate stiff matrix K of each cell  3c   xy-mode

function [k111,k112,k113,k114,k115,k116,k117,k118,k121,k122,k123,k124,k125,k126,k127,k128,...
    k131,k132,k133,k134,k135,k136,k137,k138,k141,k142,k143,k144,k145,k146,k147,k148,...
    k151,k152,k153,k154,k155,k156,k157,k158,k161,k162,k163,k164,k165,k166,k167,k168,...
    k171,k172,k173,k174,k175,k176,k177,k178,k181,k182,k183,k184,k185,k186,k187,k188,...
    k211,k212,k213,k214,k215,k216,k217,k218,k221,k222,k223,k224,k225,k226,k227,k228,...
    k231,k232,k233,k234,k235,k236,k237,k238,k241,k242,k243,k244,k245,k246,k247,k248,...
    k251,k252,k253,k254,k255,k256,k257,k258,k261,k262,k263,k264,k265,k266,k267,k268,...
    k271,k272,k273,k274,k275,k276,k277,k278,k281,k282,k283,k284,k285,k286,k287,k288,...
    k311,k312,k313,k314,k315,k316,k317,k318,k321,k322,k323,k324,k325,k326,k327,k328,...
    k331,k332,k333,k334,k335,k336,k337,k338,k341,k342,k343,k344,k345,k346,k347,k348,...
    k351,k352,k353,k354,k355,k356,k357,k358,k361,k362,k363,k364,k365,k366,k367,k368,...
    k371,k372,k373,k374,k375,k376,k377,k378,k381,k382,k383,k384,k385,k386,k387,k388,...
    k411,k412,k413,k414,k415,k416,k417,k418,k421,k422,k423,k424,k425,k426,k427,k428,...
    k431,k432,k433,k434,k435,k436,k437,k438,k441,k442,k443,k444,k445,k446,k447,k448,...
    k451,k452,k453,k454,k455,k456,k457,k458,k461,k462,k463,k464,k465,k466,k467,k468,...
    k471,k472,k473,k474,k475,k476,k477,k478,k481,k482,k483,k484,k485,k486,k487,k488]=wglmxy3cf_k(kmat,kinc,kinc3,N,ib,ie,i3b,i3e);
%K of each 4 adjacent cell
k111(1:N,1:N)=kmat(1,1);k121(1:N,1:N)=kmat(2,1);k131(1:N,1:N)=kmat(3,1);k141(1:N,1:N)=kmat(4,1);
k112(1:N,1:N)=kmat(1,2);k122(1:N,1:N)=kmat(2,2);k132(1:N,1:N)=kmat(3,2);k142(1:N,1:N)=kmat(4,2);
k113(1:N,1:N)=kmat(1,3);k123(1:N,1:N)=kmat(2,3);k133(1:N,1:N)=kmat(3,3);k143(1:N,1:N)=kmat(4,3);
k114(1:N,1:N)=kmat(1,4);k124(1:N,1:N)=kmat(2,4);k134(1:N,1:N)=kmat(3,4);k144(1:N,1:N)=kmat(4,4);
k115(1:N,1:N)=kmat(1,5);k125(1:N,1:N)=kmat(2,5);k135(1:N,1:N)=kmat(3,5);k145(1:N,1:N)=kmat(4,5);
k116(1:N,1:N)=kmat(1,6);k126(1:N,1:N)=kmat(2,6);k136(1:N,1:N)=kmat(3,6);k146(1:N,1:N)=kmat(4,6);
k117(1:N,1:N)=kmat(1,7);k127(1:N,1:N)=kmat(2,7);k137(1:N,1:N)=kmat(3,7);k147(1:N,1:N)=kmat(4,7);
k118(1:N,1:N)=kmat(1,8);k128(1:N,1:N)=kmat(2,8);k138(1:N,1:N)=kmat(3,8);k148(1:N,1:N)=kmat(4,8);
k151(1:N,1:N)=kmat(5,1);k161(1:N,1:N)=kmat(6,1);k171(1:N,1:N)=kmat(7,1);k181(1:N,1:N)=kmat(8,1);
k152(1:N,1:N)=kmat(5,2);k162(1:N,1:N)=kmat(6,2);k172(1:N,1:N)=kmat(7,2);k182(1:N,1:N)=kmat(8,2);
k153(1:N,1:N)=kmat(5,3);k163(1:N,1:N)=kmat(6,3);k173(1:N,1:N)=kmat(7,3);k183(1:N,1:N)=kmat(8,3);
k154(1:N,1:N)=kmat(5,4);k164(1:N,1:N)=kmat(6,4);k174(1:N,1:N)=kmat(7,4);k184(1:N,1:N)=kmat(8,4);
k155(1:N,1:N)=kmat(5,5);k165(1:N,1:N)=kmat(6,5);k175(1:N,1:N)=kmat(7,5);k185(1:N,1:N)=kmat(8,5);
k156(1:N,1:N)=kmat(5,6);k166(1:N,1:N)=kmat(6,6);k176(1:N,1:N)=kmat(7,6);k186(1:N,1:N)=kmat(8,6);
k157(1:N,1:N)=kmat(5,7);k167(1:N,1:N)=kmat(6,7);k177(1:N,1:N)=kmat(7,7);k187(1:N,1:N)=kmat(8,7);
k158(1:N,1:N)=kmat(5,8);k168(1:N,1:N)=kmat(6,8);k178(1:N,1:N)=kmat(7,8);k188(1:N,1:N)=kmat(8,8);
%%%%2
k211(1:N,1:N)=kmat(1,1);k221(1:N,1:N)=kmat(2,1);k231(1:N,1:N)=kmat(3,1);k241(1:N,1:N)=kmat(4,1);
k212(1:N,1:N)=kmat(1,2);k222(1:N,1:N)=kmat(2,2);k232(1:N,1:N)=kmat(3,2);k242(1:N,1:N)=kmat(4,2);
k213(1:N,1:N)=kmat(1,3);k223(1:N,1:N)=kmat(2,3);k233(1:N,1:N)=kmat(3,3);k243(1:N,1:N)=kmat(4,3);
k214(1:N,1:N)=kmat(1,4);k224(1:N,1:N)=kmat(2,4);k234(1:N,1:N)=kmat(3,4);k244(1:N,1:N)=kmat(4,4);
k215(1:N,1:N)=kmat(1,5);k225(1:N,1:N)=kmat(2,5);k235(1:N,1:N)=kmat(3,5);k245(1:N,1:N)=kmat(4,5);
k216(1:N,1:N)=kmat(1,6);k226(1:N,1:N)=kmat(2,6);k236(1:N,1:N)=kmat(3,6);k246(1:N,1:N)=kmat(4,6);
k217(1:N,1:N)=kmat(1,7);k227(1:N,1:N)=kmat(2,7);k237(1:N,1:N)=kmat(3,7);k247(1:N,1:N)=kmat(4,7);
k218(1:N,1:N)=kmat(1,8);k228(1:N,1:N)=kmat(2,8);k238(1:N,1:N)=kmat(3,8);k248(1:N,1:N)=kmat(4,8);
k251(1:N,1:N)=kmat(5,1);k261(1:N,1:N)=kmat(6,1);k271(1:N,1:N)=kmat(7,1);k281(1:N,1:N)=kmat(8,1);
k252(1:N,1:N)=kmat(5,2);k262(1:N,1:N)=kmat(6,2);k272(1:N,1:N)=kmat(7,2);k282(1:N,1:N)=kmat(8,2);
k253(1:N,1:N)=kmat(5,3);k263(1:N,1:N)=kmat(6,3);k273(1:N,1:N)=kmat(7,3);k283(1:N,1:N)=kmat(8,3);
k254(1:N,1:N)=kmat(5,4);k264(1:N,1:N)=kmat(6,4);k274(1:N,1:N)=kmat(7,4);k284(1:N,1:N)=kmat(8,4);
k255(1:N,1:N)=kmat(5,5);k265(1:N,1:N)=kmat(6,5);k275(1:N,1:N)=kmat(7,5);k285(1:N,1:N)=kmat(8,5);
k256(1:N,1:N)=kmat(5,6);k266(1:N,1:N)=kmat(6,6);k276(1:N,1:N)=kmat(7,6);k286(1:N,1:N)=kmat(8,6);
k257(1:N,1:N)=kmat(5,7);k267(1:N,1:N)=kmat(6,7);k277(1:N,1:N)=kmat(7,7);k287(1:N,1:N)=kmat(8,7);
k258(1:N,1:N)=kmat(5,8);k268(1:N,1:N)=kmat(6,8);k278(1:N,1:N)=kmat(7,8);k288(1:N,1:N)=kmat(8,8);
%%%%3
k311(1:N,1:N)=kmat(1,1);k321(1:N,1:N)=kmat(2,1);k331(1:N,1:N)=kmat(3,1);k341(1:N,1:N)=kmat(4,1);
k312(1:N,1:N)=kmat(1,2);k322(1:N,1:N)=kmat(2,2);k332(1:N,1:N)=kmat(3,2);k342(1:N,1:N)=kmat(4,2);
k313(1:N,1:N)=kmat(1,3);k323(1:N,1:N)=kmat(2,3);k333(1:N,1:N)=kmat(3,3);k343(1:N,1:N)=kmat(4,3);
k314(1:N,1:N)=kmat(1,4);k324(1:N,1:N)=kmat(2,4);k334(1:N,1:N)=kmat(3,4);k344(1:N,1:N)=kmat(4,4);
k315(1:N,1:N)=kmat(1,5);k325(1:N,1:N)=kmat(2,5);k335(1:N,1:N)=kmat(3,5);k345(1:N,1:N)=kmat(4,5);
k316(1:N,1:N)=kmat(1,6);k326(1:N,1:N)=kmat(2,6);k336(1:N,1:N)=kmat(3,6);k346(1:N,1:N)=kmat(4,6);
k317(1:N,1:N)=kmat(1,7);k327(1:N,1:N)=kmat(2,7);k337(1:N,1:N)=kmat(3,7);k347(1:N,1:N)=kmat(4,7);
k318(1:N,1:N)=kmat(1,8);k328(1:N,1:N)=kmat(2,8);k338(1:N,1:N)=kmat(3,8);k348(1:N,1:N)=kmat(4,8);
k351(1:N,1:N)=kmat(5,1);k361(1:N,1:N)=kmat(6,1);k371(1:N,1:N)=kmat(7,1);k381(1:N,1:N)=kmat(8,1);
k352(1:N,1:N)=kmat(5,2);k362(1:N,1:N)=kmat(6,2);k372(1:N,1:N)=kmat(7,2);k382(1:N,1:N)=kmat(8,2);
k353(1:N,1:N)=kmat(5,3);k363(1:N,1:N)=kmat(6,3);k373(1:N,1:N)=kmat(7,3);k383(1:N,1:N)=kmat(8,3);
k354(1:N,1:N)=kmat(5,4);k364(1:N,1:N)=kmat(6,4);k374(1:N,1:N)=kmat(7,4);k384(1:N,1:N)=kmat(8,4);
k355(1:N,1:N)=kmat(5,5);k365(1:N,1:N)=kmat(6,5);k375(1:N,1:N)=kmat(7,5);k385(1:N,1:N)=kmat(8,5);
k356(1:N,1:N)=kmat(5,6);k366(1:N,1:N)=kmat(6,6);k376(1:N,1:N)=kmat(7,6);k386(1:N,1:N)=kmat(8,6);
k357(1:N,1:N)=kmat(5,7);k367(1:N,1:N)=kmat(6,7);k377(1:N,1:N)=kmat(7,7);k387(1:N,1:N)=kmat(8,7);
k358(1:N,1:N)=kmat(5,8);k368(1:N,1:N)=kmat(6,8);k378(1:N,1:N)=kmat(7,8);k388(1:N,1:N)=kmat(8,8);
%%%%4
k411(1:N,1:N)=kmat(1,1);k421(1:N,1:N)=kmat(2,1);k431(1:N,1:N)=kmat(3,1);k441(1:N,1:N)=kmat(4,1);
k412(1:N,1:N)=kmat(1,2);k422(1:N,1:N)=kmat(2,2);k432(1:N,1:N)=kmat(3,2);k442(1:N,1:N)=kmat(4,2);
k413(1:N,1:N)=kmat(1,3);k423(1:N,1:N)=kmat(2,3);k433(1:N,1:N)=kmat(3,3);k443(1:N,1:N)=kmat(4,3);
k414(1:N,1:N)=kmat(1,4);k424(1:N,1:N)=kmat(2,4);k434(1:N,1:N)=kmat(3,4);k444(1:N,1:N)=kmat(4,4);
k415(1:N,1:N)=kmat(1,5);k425(1:N,1:N)=kmat(2,5);k435(1:N,1:N)=kmat(3,5);k445(1:N,1:N)=kmat(4,5);
k416(1:N,1:N)=kmat(1,6);k426(1:N,1:N)=kmat(2,6);k436(1:N,1:N)=kmat(3,6);k446(1:N,1:N)=kmat(4,6);
k417(1:N,1:N)=kmat(1,7);k427(1:N,1:N)=kmat(2,7);k437(1:N,1:N)=kmat(3,7);k447(1:N,1:N)=kmat(4,7);
k418(1:N,1:N)=kmat(1,8);k428(1:N,1:N)=kmat(2,8);k438(1:N,1:N)=kmat(3,8);k448(1:N,1:N)=kmat(4,8);
k451(1:N,1:N)=kmat(5,1);k461(1:N,1:N)=kmat(6,1);k471(1:N,1:N)=kmat(7,1);k481(1:N,1:N)=kmat(8,1);
k452(1:N,1:N)=kmat(5,2);k462(1:N,1:N)=kmat(6,2);k472(1:N,1:N)=kmat(7,2);k482(1:N,1:N)=kmat(8,2);
k453(1:N,1:N)=kmat(5,3);k463(1:N,1:N)=kmat(6,3);k473(1:N,1:N)=kmat(7,3);k483(1:N,1:N)=kmat(8,3);
k454(1:N,1:N)=kmat(5,4);k464(1:N,1:N)=kmat(6,4);k474(1:N,1:N)=kmat(7,4);k484(1:N,1:N)=kmat(8,4);
k455(1:N,1:N)=kmat(5,5);k465(1:N,1:N)=kmat(6,5);k475(1:N,1:N)=kmat(7,5);k485(1:N,1:N)=kmat(8,5);
k456(1:N,1:N)=kmat(5,6);k466(1:N,1:N)=kmat(6,6);k476(1:N,1:N)=kmat(7,6);k486(1:N,1:N)=kmat(8,6);
k457(1:N,1:N)=kmat(5,7);k467(1:N,1:N)=kmat(6,7);k477(1:N,1:N)=kmat(7,7);k487(1:N,1:N)=kmat(8,7);
k458(1:N,1:N)=kmat(5,8);k468(1:N,1:N)=kmat(6,8);k478(1:N,1:N)=kmat(7,8);k488(1:N,1:N)=kmat(8,8);
%%%%%inc3
%%%1
k111(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,1);k121(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,1);
k112(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,2);k122(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,2);
k113(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,3);k123(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,3);
k114(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,4);k124(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,4);
k115(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,5);k125(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,5);
k116(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,6);k126(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,6);
k117(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,7);k127(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,7);
k118(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(1,8);k128(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(2,8);
k131(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,1);k141(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,1);
k132(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,2);k142(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,2);
k133(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,3);k143(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,3);
k134(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,4);k144(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,4);
k135(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,5);k145(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,5);
k136(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,6);k146(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,6);
k137(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,7);k147(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,7);
k138(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(3,8);k148(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(4,8);
k151(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,1);k161(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,1);
k152(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,2);k162(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,2);
k153(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,3);k163(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,3);
k154(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,4);k164(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,4);
k155(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,5);k165(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,5);
k156(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,6);k166(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,6);
k157(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,7);k167(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,7);
k158(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(5,8);k168(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(6,8);
k171(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,1);k181(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,1);
k172(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,2);k182(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,2);
k173(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,3);k183(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,3);
k174(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,4);k184(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,4);
k175(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,5);k185(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,5);
k176(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,6);k186(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,6);
k177(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,7);k187(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,7);
k178(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(7,8);k188(i3b*N:i3e*N,i3b*N:i3e*N)=kinc3(8,8);
%%%2
k211(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,1);k221(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,1);
k212(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,2);k222(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,2);
k213(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,3);k223(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,3);
k214(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,4);k224(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,4);
k215(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,5);k225(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,5);
k216(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,6);k226(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,6);
k217(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,7);k227(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,7);
k218(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(1,8);k228(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(2,8);
k231(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,1);k241(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,1);
k232(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,2);k242(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,2);
k233(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,3);k243(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,3);
k234(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,4);k244(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,4);
k235(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,5);k245(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,5);
k236(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,6);k246(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,6);
k237(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,7);k247(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,7);
k238(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(3,8);k248(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(4,8);
k251(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,1);k261(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,1);
k252(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,2);k262(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,2);
k253(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,3);k263(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,3);
k254(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,4);k264(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,4);
k255(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,5);k265(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,5);
k256(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,6);k266(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,6);
k257(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,7);k267(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,7);
k258(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(5,8);k268(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(6,8);
k271(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,1);k281(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,1);
k272(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,2);k282(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,2);
k273(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,3);k283(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,3);
k274(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,4);k284(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,4);
k275(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,5);k285(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,5);
k276(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,6);k286(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,6);
k277(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,7);k287(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,7);
k278(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(7,8);k288(i3b*N+1:i3e*N+1,i3b*N:i3e*N)=kinc3(8,8);
%%%3
k311(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,1);k321(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,1);
k312(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,2);k322(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,2);
k313(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,3);k323(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,3);
k314(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,4);k324(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,4);
k315(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,5);k325(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,5);
k316(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,6);k326(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,6);
k317(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,7);k327(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,7);
k318(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(1,8);k328(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(2,8);
k331(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,1);k341(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,1);
k332(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,2);k342(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,2);
k333(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,3);k343(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,3);
k334(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,4);k344(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,4);
k335(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,5);k345(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,5);
k336(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,6);k346(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,6);
k337(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,7);k347(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,7);
k338(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(3,8);k348(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(4,8);
k351(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,1);k361(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,1);
k352(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,2);k362(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,2);
k353(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,3);k363(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,3);
k354(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,4);k364(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,4);
k355(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,5);k365(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,5);
k356(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,6);k366(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,6);
k357(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(5,7);k367(i3b*N+1:i3e*N+1,i3b*N+1:i3e*N+1)=kinc3(6,7);