function [w,y]=diedai(d,m,k,n,u)
format long;
for j=1:250
    u=d*u;
    ma=u(1,1);
    for l=1:1:n
        if ma<=u(l,1)
        ma=u(l,1);
        end
    end
    u=u/ma;
    w=sqrt((u'*k*u)/(u'*m*u));
end
w=w;
y=u/sqrt(u'*m*u);