function [f,dfdx,dfdy] = lin(x,y)

f = y - x - 1;
dfdx = -1;
dfdy = 1;