function [rz] = rz(x1,y1,z1,a,b,c)
rz = (b + (x1-c)*z1) ;