ANEW -subplot_demo

NEEDS -gnuplot

DOC
(*
	def f(t):
	    s1 = cos(2*pi*t)
	    e1 = exp(-t)
	    return multiply(s1,e1)

	t1 = arange(0.0, 5.0, 0.1)
	t2 = arange(0.0, 5.0, 0.02)
	t3 = arange(0.0, 2.0, 0.01)

	subplot(211)
	l = plot(t1, f(t1), 'bo', t2, f(t2), 'k--', markerfacecolor='green')
	grid(True)
	title('A tale of 2 subplots')
	ylabel('Damped oscillation')

	subplot(212)
	plot(t3, cos(2*pi*t3), 'r.')
	grid(True)
	xlabel('time (s)')
	ylabel('Undamped')
	show()

*)
ENDDOC

: f(t) ( F: t -- r ) FDUP FNEGATE FEXP   FSWAP PI*2 F* FCOS  F* ;
: g(t) ( F: t -- r ) PI*2 F* FCOS ;

: subplot_demo ( -- )
	0 FIGURE CLF  TRUE GRID
	S" A tale of 2 subplots" TITLE
	2 1 
	SUBPLOT
	  S" Damped oscillation" YLABEL
	  [PHOLD 
	    0e 5e S" blue--"   ['] f(t) #250 FPLOT
	    0e 5e S" !greeno"  ['] f(t)  #50 FPLOT
	  PHOLD]

	  S" time (s)" XLABEL
	  S" Undamped" YLABEL
	  [PHOLD
	    0e 2e S" @!reds"   ['] g(t) #50  FPLOT ( '.' not supported )
	    0e 2e S" red"      ['] g(t) #100 FPLOT 
	  PHOLD]
	SPEND ;