层次分析法matlab程序

clc
a=[1,1/2,4,3,3
2,1,7,5,5
1/4,1/7,1,1/2,1/3
1/3,1/5,2,1,1
1/3,1/5,3,1,1];
[x,y]=eig(a);eigenvalue=diag(y);lamda=eigenvalue(1);
ci1=(lamda-5)/4;cr1=ci1/1.12
w1=x(:,1)/sum(x(:,1))
b1=[1  2  5;1/2  1  2;1/5  1/2  1];
[x,y]=eig(b1);eigenvalue=diag(y);lamda=eigenvalue(1);
ci21=(lamda-3)/2;cr21=ci21/0.58
w21=x(:,1)/sum(x(:,1))
b2=[1  1/3   1/8;3   1   1/3;8    3    1];
[x,y]=eig(b2);eigenvalue=diag(y);lamda=eigenvalue(1);
ci22=(lamda-3)/2;cr22=ci22/0.58
w22=x(:,1)/sum(x(:,1))
b3=[1    1   3;1   1    3;1/3   1/3   1];
[x,y]=eig(b3);eigenvalue=diag(y);lamda=eigenvalue(1);
ci23=(lamda-3)/2;cr23=ci23/0.58
w23=x(:,1)/sum(x(:,1))
b4=[1   3   4;1/3   1    1;1/4  1   1];
[x,y]=eig(b4);eigenvalue=diag(y);lamda=eigenvalue(1);
ci24=(lamda-3)/2;cr24=ci24/0.58
w24=x(:,1)/sum(x(:,1))
b5=[1   1    1/4;1   1    1/4;4   4   1];
[x,y]=eig(b5);eigenvalue=diag(y);lamda=eigenvalue(1);
ci25=(lamda-3)/2;cr25=ci25/0.58
w25=x(:,2)/sum(x(:,2))
w_sum=[w21,w22,w23,w24,w25]*w1
ci=[ci21,ci22,ci23,ci24,ci25];
cr=ci*w1/sum(0.58*w1)
运行结果是
cr1 =

    0.0161


w1 =

    0.2636
    0.4758
    0.0538
    0.0981
    0.1087


cr21 =

    0.0048


w21 =

    0.5954
    0.2764
    0.1283


cr22 =

    0.0013


w22 =

    0.0819
    0.2363
    0.6817


cr23 =

   -2.5862


w23 =

    3.0000
   -1.5000
   -0.5000


cr24 =

    0.0079


w24 =

    0.6337
    0.1919
    0.1744


cr25 =

   -2.5862


w25 =

    0.1667
    0.1667
    0.6667


w_sum =

    0.4376
    0.1415
    0.4208


cr =

   -0.4176
标准的结果是

权　   准则层
值
决策层	B1     B2     B3     B4    B5	组合权向量
 

	0.2636  0.4758  0.0538  0.0981  0.1087	
 
0.5954  0.0819   0.429   0.6337   0.166	 

 
0.2764  0.2363   0.429   0.1919   0.166	 

 
0.1283  0.6817   0.142   0.1744   0.668	 

 
3.0055   3.0015     3     3.0092    3	
 
0.00275  0.00075    0     0.0046    0	
 
0.58     0.58      0.58    0.58   0.58	
 
0.0047  0.0013      0      0.0079    0	
哪里错了？帮我改！