I feel like I set expectations that I could just "do it for you", but we are clearly having difficulty. Also, it turns out, I'll be largely offline for the next week. So, I apologize for overpromising. If you can wait for awhile, and if you can find the actual syntactically correct running code, then maybe we can work this out later. But I'm starting to wonder if it wouldn't be easier for you to just implement the functionality directly in Python from scratch.

This line

M = Import[NotebookDirectory[] <> "file2.xlsx", {"Sheets", "Sheet1"}]

means this is dependent on a file, and obviously I don't have that file.

And there is still something wrong with the last function definition, the one for f[tarB_, tarG_, RhoBG_].

Have you actually successfully run any of this code? If so, something is getting lost between copying it from the notebook to pasting it in this forum.

If that's the entire code, then we have some problems. There are many undefined variables:

`dist`, `dB`, `dG`, `iJDP`, `K`, `t`, and `B`

The product t B is probably an error and was intended to be a single symboltB.

There is also a syntax problem in the definition of f. I'm not exactly sure what is being attempted, so I don't know if it's just missing brackets or missing lines of code.

So, I can't even get this to run properly.

The variable RhoLG was declared inside the Module. Mathematica uses a re-naming scheme to scope variables, so inside the Module the variable that looks to you like RhoLG is really RhoLG$19861474.

(The bit after the $ will be different each time the Module executes.

But, if you put the Print statement inside the Module it should work fine. So, you must have been putting the Print outside the Module. Since you want intermediate results, put the Print right next to each place that RhoLG is defined or updated.

If you're still having trouble, I'll try to put together a demo. But I don't have much time right now, so it'll have to wait until later (or maybe someone else can jump in and help you out).

Many thanks Eric, About the intermediate results, so for instance a step code that extract the RhoLG how should it be?' I am not very confident with the syntax of Mathematica and I would like to extract the intermediate results to compare them with the ones in Python

Many many thanks!!!!!!!!!! And if I would like to export as xlsx or csv or txt the RhoLG which step should I create? Still many thanks!!