1*-----------------------------------------------------------------------
2*     Copyright (C) 2000-2004 GFD Dennou Club. All rights reserved.
3*-----------------------------------------------------------------------
4      PROGRAM USPK09
5
6      PARAMETER(N=200, M=5)
7      REAL X(N), Y(N), YC(N,M), A(M)
8
9      DT = 1./(N-1)
10      PI = 3.14159
11      DO 50 J=1, M
12        JJ = J*2-1
13        A(J) = (-1)**J *2./(JJ*PI)
14  50  CONTINUE
15
16      DO 100 I=1, N
17        T    = DT*(I-1)*2*PI
18        X(I) = DT*(I-1)
19        Y(I) = 0.
20        DO 150 J=1, M
21          JJ = J*2-1
22          YC(I,J) = A(J)*COS(JJ*T)
23          Y(I)  = Y(I) + YC(I,J)
24  150   CONTINUE
25  100 CONTINUE
26
27*--------------------------- 1ST PAGE ----------------------------------
28
29      WRITE(*,*) ' WORKSTATION ID (I)  ? ;'
30      CALL SGPWSN
31      READ (*,*) IWS
32
33      CALL GLRGET('RUNDEF',RUNDEF)
34
35      CALL GROPN(IWS)
36      CALL GRFRM
37      CALL GRSVPT(0.2, 0.8, 0.2, 0.6)
38      CALL USSPNT(N*M, RUNDEF, YC)
39
40      CALL USSTTL('X-AXIS', ' ', 'COMPONENTS', ' ')
41      CALL USGRPH(N, X, YC)
42
43      DO 200 J=2,M
44        IP = MOD(J-1,4) + 1
45        CALL UULINZ(N, X, YC(1,J), 1, IP)
46  200 CONTINUE
47
48*     --- NEW FIG ---
49      CALL GRFIG
50      CALL GRSVPT(0.2, 0.8, 0.62, 0.82)
51      CALL UZLSET('LABELXB', .FALSE.)
52
53      CALL USSTTL('X-AXIS', ' ', 'TOTAL', ' ')
54      CALL USGRPH(N, X, Y)
55
56*--------------------------- 2ND PAGE ----------------------------------
57
58      CALL GRFRM
59      CALL GRSVPT(0.2, 0.8, 0.2, 0.6)
60
61      CALL UZFACT(0.5)
62      CALL UZLSET('LABELXB', .TRUE.)
63
64      CALL USSPNT(N*M, RUNDEF, YC)
65      CALL USSTTL('X-AXIS', ' ', 'COMPONENTS', ' ')
66      CALL USGRPH(N, X, YC)
67
68      DO 300 J=2,M
69        IP = MOD(J-1,4) + 1
70        CALL UULINZ(N, X, YC(1,J), 1, IP)
71  300 CONTINUE
72
73*     --- NEW FIG ---
74      CALL GRFIG
75      CALL GRSVPT(0.2, 0.8, 0.62, 0.82)
76      CALL UZLSET('LABELXB', .FALSE.)
77
78      CALL USSTTL('X-AXIS', ' ', 'TOTAL', ' ')
79      CALL USGRPH(N, X, Y)
80
81      CALL GRCLS
82
83      END
84