Multiplication de polynômes de matrices Friday 6 April 2007, by Cayrel Ce programme effectue la multiplication de polynômes de matrices mult_pol := proc(A,B) local E,F,P,i,j,k,C,p,n; n := rowdim(lcoeff(B)); k := degree(A) + egree(B); E := array(1..k+1,1,..k+1): for i from 1 to k+1 do for j from 1 to k+1 do E[j,i] := coeff(B,x,i-j); od; od; C := array(1..k+1): for i from 1 to k+1 do C[i] := coeff(A,x,i-1); od; P := array(1..k+1): for i from 1 to k+1 do p := sum(’evalm(C[j]&*E[j,i])’,’j’=1..k+1); P[i] := map(x-> x mod 2, evalm(p)); if op(convert(evalm(P[i]),set))=0 then P[i] := 0;fi; od; C := sum(’evalm(P[i+1])*x^i’,’i’=0..k); C; end;
RSS 2.0