Included the NTV calculation in the matlab multiple species code /HS