Question: The output only "print" not results

Hi,

why the output for this program just only "print BChange" not a matrix?

 

Der:=proc(A,n)

local i,j,k,t,S1,S2,l,C,sols,eqns,BChange;

C:=matrix(n,n); BChange:=matrix(n,n); eqns:={};

for i to n-1 do

  for j from i+1 to n do

    for l to n do

      S1:=sum(A[i,j,k]*C[k,l],k=1..n);

      S2:=sum(A[k,j,l]*C[i,k]+A[i,k,l]*C[j,k],k=1..n);

      eqns:=union(eqns,{S1=S2})

    end do

  end do

end do;

 

sols:=[solve(eqns)];

t:=nops(sols);

for i to t do

  for j to n do

    for k to n do

      BChange[j,k]:=subs(sols[i],C[j,k])

    end do;

  end do;

end do;

print (BChange)

end proc:

 

> A1 := array(sparse, 1 .. 2, 1 .. 2, 1 .. 2, [(1, 1, 2) = 1]):
> Der(A1, 2);
                                                                   print BChange

 

Can someone please advice me on this?

thanks

witribm

 

Please Wait...