ANEW -fill_demo

NEEDS -gaussj
NEEDS -gnuplot

DOC
(*
	#!/usr/bin/env python
	import numpy as np
	import matplotlib.pyplot as plt

	t = np.arange(0.0, 1.01, 0.01)
	s = np.sin(2*2*np.pi*t)

	plt.fill(t, s*np.exp(-5*t), 'r')
	plt.grid(True)
	plt.show()
*)
ENDDOC

: fill_demo ( -- )
	#101 DOUBLE LARRAY x{
	#101 DOUBLE LARRAY y{
	0 FIGURE CLF TRUE GRID
	x{ 0e 1.01e #101 LINSPACE
	#101 0 ?DO  	I S>F 0.01e F*  FDUP  x{ I } DF!
			FDUP  PI*2 F* F2* FSIN    
			FSWAP -5e F* FEXP F*  y{ I } DF!  
	      LOOP
	x{ y{ S" red" PFILL ;