eq1:h=sqrt(R^2-(ix-t*dx)^2-(iy-t*dy)^2)+t*dz;

eq2:diff(eq1,t);

eq3:eq2-dz;

eq4:eq3^2;

eq5:solve(eq4,t);

eq6:factor(eq5);