blob: f25e7ebf20cbd365d36b2a47b59ffeb3164ab30f (
plain)
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
|
c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
c Copyright (C) ????-2008 - INRIA
c
c This file must be used under the terms of the CeCILL.
c This source file is licensed as described in the file COPYING, which
c you should have received as part of this distribution. The terms
c are also available at
c http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
subroutine dpodiv(a,b,na,nb )
implicit double precision (a-h,o-z)
dimension a(*),b(*)
c
c 1. a: avant execution c'est le tableau des coefficients du
c polynome dividende range suivant les puissances croissantes
c de la variable (ordre na+1).
c apres execution,contient dans les nb premiers elements
c le tableau des coefficients du reste ordonne suivant les
c puissances croissantes, et dans les (na-nb+1) derniers elements,
c le tableau des coefficients du polynome quotient range suivant
c les puissances croissantes de la variable.
c 2. b: tableau des coefficients du polynome diviseur range suivant
c les puissances croissantes de la variable (ordre nb+1).
c 3. na: degre du polynome a.
c 4. nb: degre du polynome b.
c
l=na-nb+1
2 if (l .le. 0) then
goto 5
else
goto 3
endif
3 n=l+nb
q=a(n)/b(nb+1)
nb1=nb+1
do 4 i=1,nb1
n1=nb-i+2
n2=n-i+1
4 a(n2)=a(n2)-b(n1)*q
a(n)=q
l=l-1
goto 2
5 continue
return
end
|
