function vv = frestr(y)
global coor0, global vec0, global B, global vol0,
vv=zeros(2,1);
if (rows(y)==1)
 y = y';
endif
yd=diag(y);
coor=coor0+yd*vec0;
vv(1)=volume(coor)-vol0;
vv(2)=y(1);
end