Tips for converting this piece of f77 source code to Matlab

The code attached is a piece of from a publicly available FORTRAN source (the complete source code is found here ). I was wondering if anyone could direct me on how to convert it to Matlab code. I've tried to convert it using f2matlab, however I can't convert it to f90 (which f2matlab requires).

The Matlab file I've attached is the conversion I've got so far, but I get an error with "while(more)".

Any help would be greatly appreciated.

Answer 1

[Stanley Kubrick]

Source code

Answer 2

[Stanley Kubrick]

Matlab Conversion so far

Answer 3

[Star Strider]

The matrix exponential function expm and matrix left divide function (\) mldivide (to use instead of the explicit inverse) are your friends here. I’m not certain what you’re doing with your FORTRAN subroutines, so I can’t comment further.