We present a new Vlasov code for collisionless plasmas in the nonrelativistic regime. A Darwin approximation is used for suppressing electromagnetic vacuum modes. The spatial integration is based on an extension of the flux-conservative scheme, introduced by Filbet et al. [J. Comp. Phys. Vol. 172 (2001) 166]. Performance and accuracy is demonstrated by comparing it to a standard finite differences scheme for two test cases, including a Harris sheet magnetic reconnection scenario. This comparison suggests that the presented scheme is a promising alternative to finite difference schemes.