1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
// function [A,degA] = putin(A,degA,B,degB,i,j)
// PUTS IN THE POLYNOMIAL B INTO THE MATRIX A AT (i,j)TH POSITION
// Modified by Kannan Moudgalya on 8 November 1992
function [A,degA] = putin(A,degA,B,degB,i,j)
[rA,cA] = polsize(A,degA);
if degB > degA
A = [A zeros(rA,(degB-degA)*cA)];
degA = degB;
end
for k = 0:degB
A(i,(k*cA)+j) = B(1,k+1);
end
if degA > degB
for k = degB+1:degA
A(i,(k*cA)+j) = 0;
end
[A,degA] = clcoef(A,degA);
end
endfunction;
|
