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