;pro plot_runge

;NOTE: THIS PLOTTING PROGRAM ASSUMES THAT THE FINAL 
;FOURTH COLUMN IN THE FILE "runge.out" IS Y(X).

a=fltarr(4,100000) & b=fltarr(4) & stuff=strarr(1)
openr,2,'runge.out' & readf,2,stuff & readf,2,stuff
ic=0
while (not eof(2)) do begin
  readf,2,b & a(0:3,ic)=b(0:3) & ic=ic+1
endwhile
close,2

set_plot,'ps'
device,/landscape & !p.font=0
plot,a(0,0:ic-1),a(3,0:ic-1),/xs,xthick=3,ythick=3,thick=5,$
     xcharsize=1.2,ycharsize=1.2,charsize=1.2,$
     title='Numerical Solution From Program "RUNGE.F"',$
     xtitle='X',ytitle='Y(X)'
if max(a(3,*)) gt 0 and min(a(3,*)) lt 0 then begin
  oplot,[min(a(0,*)),max(a(0,*))],[0,0],thick=3
endif

device,/close

end