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