function [fy] = ly(x1,y1,z1,r,b,sig)
fy =r*x1-y1-x1*z1;