Consider the famous Van Der Pol system:
x' = y
y’ = ay(1-x^2)-x
Where a is for example 2.
Is there any way to find its limit cycle (attached picture) based on the equations and not simulating the system? Note that I don’t mean finding whether limit cycle exists or not. I mean determining its exact shape and location.