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 -  "
print "Sobral - Ceará - Brasil " 
print "Aperte enter para continuar!"
pause -2
system "clear"

print   "Instruções para rodar o programa:";
print   "Este programa faz experiências com quatro tipos de funções";
print   "duas funções polinomiais, segundo e terceiro grau";
print   "a função seno";
print   "a função de Heaviside  - uma função que não é contínua.";
print   "As experiências estão divididas em quatro blocos separados";
print   "pela função quit - desta forma você pode repetir pelo menos";
print   "a primeira experiência várias vezes para se beneficiar do ";
print   "efeito visual da mesma. Também poderá cortar o programa em";
print   "vários blocos assim ter o efeito visual para cada um dos casos.";
print   "Chame o programa num terminal e o afaste do terminal gráfico";
print   "que gnuplot vai criar. Assim você poderá ler as mensagens ";
print   "enquanto roda o programa. Em Linux funciona!";
print "Aperte enter para continuar!"
pause -2
system "clear"



print "********  Os operadores diferença e quociente de diferenças  ******* "
print "Notação  D_f(x) = f(x + rho) - f(x)   --  Delta de f  "
print "Notação  Q_f(x) = (f(x + rho) - f(x))/rho    --  Quoc de f  "
print "Aperte enter para continuar!"
pause -2

print "************* Diferenças e a função  seno ************* "
print "Aperte enter para continuar!"
pause -2

print "rho = 0.5";
print "g(x) = sin(x)";
print "Q_g(x) = (g(x + rho) - g(x))/rho;";  ## operador quoc. de dif
print "D_g(x) = g(x + rho) - g(x) ";          ## operador difer.
print "Aperte enter para continuar!"
pause -2

print "************* Diferenças e a função  seno ************* ";
print  "rho = 0.5";
print  "g(x) = sin(x)";
print  "Q_g(x) = (g(x + rho) - g(x))/rho;"; ## operador quoc. de dif
print  "D_g(x) = g(x + rho) - g(x) ";         ## operador difer.
print  "plot g(x), Q_g(x), D_g(x),0  ";
print "Aperte enter para continuar!"
pause -2

rho = 0.5
g(x) = sin(x)
Q_g(x) = (g(x + rho) - g(x))/rho;           ## operador quoc. de dif
D_g(x) = g(x + rho) - g(x)                    ## operador difer.
plot g(x), Q_g(x), D_g(x),0   
print "Aperte enter para continuar!"
pause -2


print "************* Diferenças e a função  seno ************* ";
print  "rho = 0.2";
print  "g(x) = sin(x)";
print  "Q_g(x) = (g(x + rho) - g(x))/rho;";  ## operador quoc. de dif
print  "D_g(x) = g(x + rho) - g(x) ";          ## operador difer.
print  "plot g(x), Q_g(x), D_g(x),0  ";
print "Aperte enter para continuar!"
pause -2

rho = 0.2
g(x) = sin(x)
Q_g(x) = (g(x + rho) - g(x))/rho;            ## operador quoc. de dif
D_g(x) = g(x + rho) - g(x)                     ## operador difer.
plot g(x), Q_g(x), D_g(x),0   
print "Aperte enter para continuar!"
pause -2


print "************* Diferenças e a função  seno ************* ";
print  "rho = 0.1";
print  "g(x) = sin(x)";
print  "Q_g(x) = (g(x + rho) - g(x))/rho;";  ## operador quoc. de dif
print  "D_g(x) = g(x + rho) - g(x) ";          ## operador difer.
print  "plot g(x), Q_g(x), D_g(x),0  ";
print "Aperte enter para continuar!"
pause -2

rho = 0.1
g(x) = sin(x)
Q_g(x) = (g(x + rho) - g(x))/rho;            ## operador quoc. de dif
D_g(x) = g(x + rho) - g(x)                     ## operador difer.
plot g(x), Q_g(x), D_g(x),0   
print "Aperte enter para continuar!"
pause -2


print "************* Diferenças e a função  seno ************* ";
print  "rho = 0.05";
print  "g(x) = sin(x)";
print  "Q_g(x) = (g(x + rho) - g(x))/rho;";  ## operador quoc. de dif
print  "D_g(x) = g(x + rho) - g(x) ";          ## operador difer.
print  "plot g(x), Q_g(x), D_g(x),0  ";
print "Aperte enter para continuar!"
pause -2

rho = 0.05
g(x) = sin(x)
Q_g(x) = (g(x + rho) - g(x))/rho;            ## operador quoc. de dif
D_g(x) = g(x + rho) - g(x)                     ## operador difer.
plot g(x), Q_g(x), D_g(x),0   
print "Aperte enter para continuar!"
pause -2


print "************* Diferenças e a função  seno ************* ";
print  "rho = 0.01";
print  "g(x) = sin(x)";
print  "Q_g(x) = (g(x + rho) - g(x))/rho;";  ## operador quoc. de dif
print  "D_g(x) = g(x + rho) - g(x) ";          ## operador difer.
print  "plot g(x), Q_g(x), D_g(x),0  ";
print "Aperte enter para continuar!"
pause -2

rho = 0.01
g(x) = sin(x)
Q_g(x) = (g(x + rho) - g(x))/rho;            ## operador quoc. de dif
D_g(x) = g(x + rho) - g(x)                     ## operador difer.
plot g(x), Q_g(x), D_g(x),0   
print "Aperte enter para terminar !"
pause -2

quit