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
55
|
*-----------------------------------------------------------------------
PROGRAM USPK07
PARAMETER(N=200, M=10)
REAL X(N), Y(N), A(M)
*-----------------------------------------------------------------------
DT = 1./(N-1)
PI = 3.14159
DO 50 J=1, M
JJ = J*2-1
A(J) = (-1)**J *2./(JJ*PI)
50 CONTINUE
DO 100 I=1, N
T = DT*(I-1)*2*PI
X(I) = DT*(I-1)
YY = 0.
DO 150 J=1, M
JJ = J*2-1
YD = A(J)*COS(JJ*T)
YY = YY + YD
150 CONTINUE
Y(I) = YY*2.
100 CONTINUE
*-----------------------------------------------------------------------
WRITE(*,*) ' WORKSTATION ID (I) ? ;'
CALL SGPWSN
READ (*,*) IWS
CALL GROPN(IWS)
CALL GRFRM
* --- X AXIS ---
CALL UZRSET('UYUSER' , 0.)
CALL USCSET('CXSIDE' , 'U')
* --- Y AXIS ---
CALL UZISET('IROTLYL' , 1)
CALL UZISET('ICENTYL' , 0)
CALL USCSET('CYSIDE' , 'L')
* --- etc. ---
CALL UZISET('INNER' , -1)
CALL USCSET('CBLKT' , '[]')
CALL USSTTL('TIME', 'SEC', 'VOLTAGE', 'mV')
CALL USGRPH(N, X, Y)
CALL GRCLS
END
|
