If the tilt is a rotation about the yaxis, then I suggest you first calculate the (x,y) coordinates of one of the satellites as above, then calculate the new (x',y',z') coordinates after a tilt of theta degrees about the yaxis as follows:
x' = x cos(theta)
y' = y
z' = x sin(theta)
This way you can have one satellite orbiting in the xy plane and the other orbiting at an angle theta relative to that.
