f(x,y,z,p,q)=p2+q2−2px−2qy+1=0
fx+pfzdp=fy+qfzdq=−pfp−qfqdz=−fpdx=−fqdy
pdp=qdq=p(p−x)+q(q−y)dz=p−xdx=q−ydy pdp=qdq⇒lnp=lnq+lnc⇒p=cq
f(x,y,z,p,q)=c2q2+q2−2cqx−2qy+1=q2(c2+1)−2q(cx+y)+1=0
q=2(c2+1)2(cx+y)±4(cx+y)2−4(c2+1) cx+y=t,c2+1=a
q=at±t2−a,p=cq=ac(t±t2−a)
dz=pdx+qdy=cqdx+qdy=q(cdx+dy)=qd(cx+y)=qdtz=∫at±t2−adt=2at2±a1∫t2−adt
z=2at2±21tt2−a∓21aln(t2−a+t)+C
z=2a(cx+y)2±21(cx+y)(cx+y)2−a∓21aln((cx+y)2−a+cx+y)+C