@Kitonum Thank you ever so much.  I have added interface(rtablesize=infinity): inside the main proc.  It worked but I am not sure if that is the right thing to do, is it?

In the old code linalg[vectdim](f) returned f's dimension.  I could not find anything similar in LinearAlgebra but nops(f) seems to work.  Would nops be ok?

Now that my old code is working I am trying to convert all deprecated functions to new ones.   The first setp was changing vector to Vector[row]. 

When trying to convert a list to Vector inside a proc, the outcome comes with a lot of info that wrecks the code.  Example:

aux:=[-x1*x2+u0,x1*x1-2*x2,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10];

f:=Vector[row](aux);

f returns [1..12 Vector rwo, Data Type: anything, ..., Order: fortran_order].   How to get a simple output as vector did?

Many thanks

Ed

@Carl Love Many many thanks.   Part of the old code is now working.   

@Carl Love Yes, but I found that there is nothing but (x,y) -> x-y.  At least it works for the cases I have tried. Many thanks.

@Carl Love Many thanks.  Would inner always be available?

@Kitonum Thank you ever so much.   There is another function that I could find in the version.  Here is an example:

(anglebracket(X-Y))(y2, -(-x1*x2+u0)*x2-(x1^2-2*x2)*x1+u1)

Unfortunately I don't have a foggiest idea where this function comes from.

I would rather not use LinearAlgebra right now.  If I do I have changing all the definitions from vector to Vector. Once I get the code working I will definitely change it.  Many thanks.