makemat(n)

[compsystems]

Hi BP, just to bring the matlab style in the makemat function, add please the following

makemat (n) returns a zeros square matrix of dimension (n * n)


makemat(0,s,s); == makemat(s);

Code : Tout sélectionner

hilbert_( s ):={
local h, c, r;
h := makemat(0,s,s);
for c from 1 to s do
    for r from 1 to s do
        h[[r,c]] := 1/(r+c-1);
    end;
end;
return h;
}

or

Code : Tout sélectionner

hilbert_( s ):={
local h, c, r;
h := makemat(s);
for c from 1 to s do
    for r from 1 to s do
        h[[r,c]] := 1/(r+c-1);
    end;
end;
return h;
}

hilbert_( 3 )
[[1,1/2,1/3],
[1/2,1/3,1/4],
[1/3,1/4,1/5]]

[parisse]

Ok, adding makemat(n) as a synonym for matrix(n,n).

[compsystems]

also for the matrix(n) equal to the matrix(n, n) the square matrices are the most used =)