[routines] Implement take_along_axis method

[forfudan]

This will:

1. Implemented take_along_axis function in numojo/routines/indexing.mojo:

   • Takes values from an input array along a specified axis based on indices
   • Handles arrays of any dimension with axis specification
   • Supports broadcasting of indices when shapes don't exactly match
   • Includes proper error handling for invalid inputs

2. Created comprehensive test suite in tests/routines/test_indexing.mojo:

   • Tests for arrays of different dimensions (1D, 2D, 3D)
   • Tests with various axis specifications (positive, negative indices)
   • Tests with different index types (int32, int64)
   • Verification against NumPy's implementation for correctness