function W = weight(D,V,L,alpha,b)
Vp=pinv(V);
r=size(V,2);
B=bigmat(b,r);
W=(D*V-B)*Vp+L*(eye(size(D))-V*Vp);