clear all;close all;clc
mu=0.00004;
lamuta1=0.2;
lamuta2=0.9;
h1=0;
h2=0.5;
h12=h2-h1;
hs=(h2*h2-h1*h1)/2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k11=-0.932;
k22=-0.231;
I6=eye(6);
......................