#! /usr/bin/gnuplot
system "clear"
print " "
print " "
print "============================================================"
print "			Os direitos do autor e a GPL  "
print "============================================================"
print "  "
print "Este programa é distribuido sob GPL. "
print "Se você não souber o que é a GPL, envie um e-mail ao autor "
print "                   tarcisio@member.ams.org "
print "mas em suma  quer dizer que você pode usar livremente o sistema"
print "de programas aqui mencionados, desde que esta informação seja mantida"
print "em todas as alterações que você fizer dos programas e das cópias que "
print "você deles fizer - é a proteção do direito autoral que não lhe impede"
print "de usar livremente o conhecimento mas que mantém na memória  quem "
print "participou de sua produção."
print "Este programa foi escrito por Tarcisio Praciano-Pereira "
print "professor da Universidade Estadual Vale do Acaraú - UeVA "
print "Laboratório de Matemática Computacional - Sobral Matematica " 
print "www.sobralmatematica.org "
print "Sobral - Ceará - Brasil " 
print "Aperte enter para continuar!"
pause -2
system "clear"


print "O programa pode ser tornado executável, depende que gnuplot esteja presente"
print "no sistema. Basta executar "
print "		chmod u+x sinal_kernel01_01.gnuplot "
print "e depois, num terminal, executar ./sinal_kernel01_01.gnuplot"
pause -2  "Aperte enter para continuar!"
system "clear"

print "Este programa é parte da solução do problema tratado em artigos sobre"
print "convolução de kernels, e serve para ilustrar pequenos detalhes neste"
print "artigos, mas também será  útil para o trabalho com funções características"
print "de retângulos em espaços de maior dimensão. "
print "  "
print "Os programas desta série tem a mesma raíz no nome \"sinal_kernel\"   "
print "e podem ser encontrados no mesmo diretório em que se encontra.   "
print "   "
pause -2 "aperte enter para continuar"
system "clear"


print "O programa contém a definição da convolução de duas funções"
print "características de retângulos centrados na origem."
print "O programa também faz o gráfico da convolução destas funções características"
print "Q(x,y) = chi_epsilon(x)*chi_epsilon(y) "
print "em que chi_epsilon(x) é a função característica de [-epsilon, epsilon]"
print "para obter a convolução de dois núcleos vou, posteriormente, "
print "multiplicar Q por 1/(4*power(epsilon,2)) "
print "ver préprint 2011.05 em www.sobralmatematica.org/preprints "
pause -2 "Aperte enter para continuar!"

print "   A definição em gnuplot da convolução:  "
print "ep = 1; ##  epsilon guardado  na variável ep"
print "set xrange [-5:5];"
print "set yrange [-5:5];"
print "Q(x,y)= (x<=-2*ep)?0:\"
print "	(x<=0)*(y<-2*ep)?0:\"
print "	(x<=0)*(y<0)?(y+2*ep)*(x+2*ep):\"
print "	(x<=0)*(y<2*ep)?(2*ep-y)*(x+2*ep):\"
print "	(x<=2*ep)*(y<-2*ep)?0:\"
print "	(x<=2*ep)*(y<0)?(y+2*ep)*(2*ep-x):\"
print "	(x<=2*ep)*(y<2*ep)?(2*ep-y)*(2*ep-x):0"

pause -2 "Aperte enter para continuar  " 

ep = 1; ##  epsilon na variável ep
set xrange [-5:5];
set yrange [-5:5];
Q(x,y)= (x<=-2*ep)?0:\
	(x<=0)*(y<-2*ep)?0:\
	(x<=0)*(y<0)?(y+2*ep)*(x+2*ep):\
	(x<=0)*(y<2*ep)?(2*ep-y)*(x+2*ep):\
	(x<=2*ep)*(y<-2*ep)?0:\
	(x<=2*ep)*(y<0)?(y+2*ep)*(2*ep-x):\
	(x<=2*ep)*(y<2*ep)?(2*ep-y)*(2*ep-x):0
splot Q(x,y)

pause -2 "Aperte enter para terminar "