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