C R.C.
      PROGRAM OSCIATION
      PARAMETER( NMAX=1000 )
      REAL Y(NMAX), T(NMAX)
C 配列を宣言するとき、配列の長さは変数で指定できない。ここでの変数は NMAX.
C そのため、配列の長さを変数で指定したい場合は、PARAMETER宣言を行う.
*-- 出力ファイルをダンプ (保存) するための命令 ---
*      CALL SWpSET( 'LDUMP', .TRUE. )
*      CALL SWpSET( 'LWAIT', .FALSE. )
*      CALL SWpSET( 'LWAIT1', .FALSE. )
*-- 時間および空間分解 ---
      DT=0.01
*-- 初期条件 ---
      Y(1)=1.0
      T(1)=0.0
*-- 支配方程式 ---
      DO 10 I=1,NMAX
         T(I+1)=T(I)+DT
         Y(I+1)=(1.0-DT)*Y(I)
  10  CONTINUE

*-- グラフの出力装置の指定 (IWS=4 は GTK 指定) ----
      WRITE(*,*) ' WORKSTATION ID (I)  ? ;'
      CALL SGPWSN
*      IWS=4
      READ (*,*) IWS

*-- 出力装置の呼び出し
      CALL GROPN( IWS )

      CALL GRFRM
      CALL USSTTL( 'TIME', '', 'DEPARTMENT', '' )
      CALL USGRPH( NMAX, T, Y )
      CALL GRSWND(  0.2,  0.8,  0.2,  0.8 )

      CALL GRCLS
      STOP
      END