duffing_lyap.c