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
|
C/MEMBR ADD NAME=SCAPOL,SSI=0
c Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
c Copyright (C) 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 scapol(na,a,nb,b,y)
c!but
c cette subroutine a pour but de calculer le produit
c scalaire de deux polynomes
c!liste d'appel
c subroutine scapol(na,a,nb,b,y)
c Entree :
c a. est le premier polynome de degre na
c b. est le second polynome du produit, et est de degre nb
c
c Sortie :
c y. est le resultat du produit scalaire <a,b>
c!
implicit double precision (a-h,o-y)
dimension a(0:*),b(0:*)
c
if (na.ge.nb) then
nmax=nb
else
nmax=na
endif
c
aux=0.0d+0
do 20 k=0,nmax
aux=aux + a(k)*b(k)
20 continue
y=aux
c
end
|
