Haskell & Scheme Code for Matrix Operations

I'm looking for a talented programmer with expertise in both Haskell and Scheme to develop code for performing matrix operations using the echelon method. The code should be able to handle determinant, rank, and inverse calculations for 1000x1000 matrices.

Key requirements:
- Proficient in both Haskell and Scheme
-Proficient in Math
- Experience with matrix operations and echelon method
- Able to write code handling 1000x1000 matrices
**In Scheme, external libraries are only allowed to be used in flomat and only for vectorization of code.
In Scheme, functions that modify the values ​​of objects are not allowed, in particular, functions with names that contain ! are not allowed.
**In Haskell, you can use Hmatrix for code vectorization only...What it means? you CANNOT use functions:
Numeric.LinearAlgebra.det
Numeric.LinearAlgebra.rank
Numeric.LinearAlgebra.inverse

Skills and experience:
- Strong programming skills
- Knowledge of linear algebra
- Experience with Haskell and Scheme programming
- Ability to write efficient code for matrix operations

Running time for both codes:
-Rank and Determinant should run at most 1 sec.
-Inverse at most 6 sec.