Dynamic analysis using numerical multi-body approach for quadruped robots