# Python script
from math import *
def f(x):
	return 1.0+floor(sin(10*x**2-0.3))
l=[ f(x/300.0) for x in range(2000) ]



def zz(l):
	o=[]
	for i in range(len(l)):
		zs=0.0
		s=0.0
		for k in range(20):
			j=k-10
			z=exp(-0.03*j**2))
			try:
				n=l[i+j]
				zs+=z
				s+=n*z
			except:
				None
		o.append(s/(1.0*zs))
		#except:	#	o.append(0.0)
	return o
out=zz(l)
out=zz(out)
out=zz(out)
out=zz(out)
out=zz(out)
out=zz(out)
out=zz(out)

for i in range(len(out)-700):
	print ""
	for j in range(100):
		print i,j,out[i+100],l[i+100] 
 # plot the result unsing gnuplot with pm3d
 # f(x)= (a*x+b)>1? 1: (a*x+b)
 # g(x)=f(x)>0.55? f(x): 0.55
 # h(x,y)=((x-0.5)/(2*(g(y)-0.5)))+0.5
 # a = -0.0006
 # b = 1.1
 # set pm3d map
 # splot "data" u 1:2:(h($3,$1))