1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
function [f1]=%r_m_p(f1,n2)
// r=%r_m_p(f1,n2) <=> r=f1*n2
// f1=rational n2=polynomial
//!
[n1,d1]=f1(["num","den"]),
[l1,m1]=size(n1);[l2,m2]=size(n2),
//
indef=%f
if l1==-1 then
n1=n1+0;d1=d1+0;l1=1;m1=1;
if l2*m2==1 then indef=%t,else error(14),end
end
if l2==-1 then
n2=n2+0;l2=1;m2=1;
if l1*m1==1 then indef=%t,else error(14),end
end
if min([l1*m1,l2*m2])==1 then,
num=n1*n2,
den=d1*ones(l2,m2),
else,
if m1<>l2 then error(10),end,
for i=1:l1,
n=n1(i,:);
[y,fact]=lcm(d1(i,:)),
den(i,1:m2)=ones(1,m2)*y,
for j=1:m2,
num(i,j)=n*(n2(:,j).*matrix(fact,l2,1)),
end,
end,
end,
if ndims(num)<=2 then
[num,den]=simp(num,den),
else
sz=size(num)
[num,den]=simp(num.entries,den.entries),
num=matrix(num,sz)
den=matrix(den,sz)
end
if indef then
num=num*eye()
den=den*eye()
end
f1=rlist(num,den,f1("dt"))
endfunction
|
