A numerical investigation of coorbital stability and libration in three dimensions