; __________________________________________________________________ ;
; MODEL3.PRO         ; The name of this program, should I forget it.
; __________________________________________________________________ ;
!P.NOERASE=1         ; Do not erase screen after each new PLOT command 
T=FINDGEN(101)/100.D0            ; t=0.00, 0.01, ...,  1.00
T=5.D0*T                         ; t=0.00, 0.05, ...,  5.00
P=1.5 & F=P^(-1.D0)              ; P fixed, f computed
W=COS(2.D0*!PI*T*F)              ; W(t)
SET_VIEWPORT,0.1,0.4,0.6,0.9     ; Panel location on screen/paper
PLOT,T,W                         ; Continuous line: PSYM=0 is default
OPLOT,T,W,PSYM=1                 ; Overplot with crosses: PSYM=1
U=FINDGEN(101)/100.D0            ; u=0.00, 0.01, ...,  1.00
U=4.D0*U-2.D0                    ; u=-2.0,-1.96, ...,1.96,2.0
G=EXP(-1.*U)                     ; g(u)
SET_VIEWPORT,0.6,0.9,0.6,0.9     ; Panel location on screen/paper
PLOT,U,G                         ; Continuous line: PSYM=0 is default
OPLOT,U,G,PSYM=4                 ; Overplot with diamonds: PSYM=4
; The X-axis and Y-axis limits are different in the two panels. 
X1=MIN([T,U]) & X2=MAX([T,U])    ; Suitable X-axis limits
Y1=MIN([W,G]) & Y2=MAX([W,G])    ; -"-      Y-axis limits
SET_XY,X1,X2,Y1,Y2               ; XY-axis limits fixed
SET_VIEWPORT,0.1,0.9,0.1,0.4     ; Panel location on screen/paper
 PLOT,T,W,PSYM=0, LINESTYLE=0    ; Continuous line
OPLOT,U,G,PSYM=0, LINESTYLE=1    ; Dotted continuous line
END
; __________________________________________________________________ ;