Numerical model for low-frequency ion waves in magnetized dusty plasmas