zu1=-θu1-(z(2v0-z)/(2v1)),