function y=vmodel(x,b)
y=b(1) +b(2).*x(:,1) +b(3).*x(:,2);
return
