x'=x*(a-b*y)
y'=y*(-c+d*x)
par a=1,b=1,c=1,d=1
aux q=a*ln(abs(y))+c*ln(abs(x))-b*y-d*x
init x=.2,y=.1
done