Quel langage pour des applications numériques scientifiques

Bonjour,

Je tente ma question dans cette rubrique ...
Ingénieur retraité, j'ai beaucoup utilisé les langages Basic, Fortran (des milliers de lignes ...) et un peu C dans les années 75-90 pour des applications scientifiques (Fourier, Navier-Stokes, ...) essentiellement dans le domaine du Spatial.

A cet époque préhistorique, les Différences Finies (DF) se sont répandues, suivies des Eléments finis (EF) et des volumes finies (VF).

A titre pédagogique je voudrais refaire un (ou des) schéma simple (1D) aux DF pour la loi de Fourier, voire éventuellement ensuite pour d'autres applications, donc toujours dans le domaine des sciences physiques.

Question : quel langage me conseilleriez-vous aujourd'hui pour programmer ce genre de problème ?

Dois-je reprendre le Fortran et rechercher un (vieux) compilateur ? ou le Basic ?
Ou vers quel langage m'orienter (C++, python, ... que sais-je ce qu'il existe aujourd'hui), de préférence vers une solution disposant d'une bibliothèque de sous-programmes scientifiques assez riche ?
L'objectif pour moi est avant tout la mise en œuvre de l'application physique mais je suis prêt à m'investir ... (Je le fait déjà fait, dans un autre domaine, en créant des sites internet avec le logiciel libre OpenElement).
J'espère l'avis circonstancié d'un spécialiste dans ce domaine...
Merci d'avance