Hello,

thanks a lot. I found a solution, probably not the most elegant one, but it works. I added:

...

result[i]:=fsolve({eq,eql},{mu,t},t=0..1,mu=-1..1):  

mures[i]:=rhs(result[i][1]):

 tres[i]:=rhs(result[i][2]):

 ires[i]:='MX[i,1]':

 end do:

And then I create a Matrix to store the all results and export the matrix to an excel file.

Thanks, Dominik.