wglmxy3cf_k.m

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

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