coupled differential equations solver