function ee = fenergy(y)
global coor0, global vec0, global B
if (rows(y)==1)
 y = y';
endif
yd=diag(y);
coor=coor0+yd*vec0;
ee=energy(coor);
endfunction