]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1996/04/01 15:02:29 mclareni | |
6 | * Mathlib gen | |
7 | * | |
8 | * | |
9 | #include "gen/pilot.h" | |
10 | FUNCTION RTRGSM(X,A,N,B,M,IOP) | |
11 | ||
12 | DIMENSION A(0:*),B(*) | |
13 | ||
14 | PARAMETER (R1 = 1, HF = R1/2) | |
15 | ||
16 | C=2*COS(X) | |
17 | H=0 | |
18 | IF(IOP .NE. 3) THEN | |
19 | A2=0 | |
20 | A1=0 | |
21 | A0=A(N) | |
22 | DO 1 K = N-1,0,-1 | |
23 | A2=A1 | |
24 | A1=A0 | |
25 | A0=A(K)+C*A1-A2 | |
26 | 1 CONTINUE | |
27 | H=HF*(A(0)+A0-A2) | |
28 | ENDIF | |
29 | IF(IOP .NE. 2) THEN | |
30 | A2=0 | |
31 | A1=0 | |
32 | A0=B(M) | |
33 | DO 2 K = M-1,1,-1 | |
34 | A2=A1 | |
35 | A1=A0 | |
36 | A0=B(K)+C*A1-A2 | |
37 | 2 CONTINUE | |
38 | H=H+A0*SIN(X) | |
39 | ENDIF | |
40 | RTRGSM=H | |
41 | RETURN | |
42 | END |