A=[1-0.1-0.10-0.2-0.4;-0.31-0.20-0.30.2;0-0.21-0.5-0.10;-0.1-0.3-0.11-0.2-0.1;-0.2-0.3-0.2-0.11-0.1;-0.3-0.1-0.1-0.2-0.11];
S=zeros（66）;
S（16）=0.4;S（25）=0.3;S（34）=0.5;S（45）=0.2;S（56）=0.1;
I=eye（66）;
AA=（I+S）*A;
L=-tril（A）;
L=L+I;
U=-triu（A）;
U=U+I;