Yes Prime can do a little bit better with units than Mathcad (e.g. in a solve block with odesolve). It also has some other advantages, like the usage of different dimensions in one matrix, matrix editing (matrix scrolling on the other hand is just a bad joke), matrix row selector operator, Ctrl-J for programming keywords.
But the handling, usability and performance speed of Prime and its user interface is inferior. And its still missing a lot of features we are used from Mathcad. So I don't consider it ready for use for real work.
> Why did you add "return" in the code?
For debug purposes. I just wanted to see which values are in use so I let the program return those values and stop (so we do not run in the error) and as you can see you run into a kind of divison by zero, or better division by a zero inside a matrix inside a matrix ...
In Mathcad we have the "trace" command inside a program and turn on debug mode and so we could see live the content of the various variables while the program runs. Prime does not offer any debug mode.
> The sum function should give as result a vector,
Your program never reaches that line because the line before you divide by this zero matrix (fugtv[i+3)
> The problem seems to be the function Z which I haven't looked at.
You had yourself evaluated it and if you turn on showing nested matrices you can see the problem - the same as I had seen when i used that return command.
The reason you have created that nasty nested matrices is because in the definition of f1 and f2 in your program you type square brackets which creates a matrix. Always use normal parentheses!
If you correct this you still have a division by zero as fug0 and fug1 result in zero.
Reason is the high values of f0 and f1.
E.g. f0=-2.0658207*10^6 and so exp(f0) is 2.848822*10^-897175 which for Mathcads numerics is zero.
Obviously there still is something else wrong.
