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 "Dados três pontos, não colineares, eles determinam um plano e a equação"
print "deste plano pode ser obtida com o seguinte algoritmo:"
print "1) Chame um dos pontos de P = (a,b,c) para considerar o ponto-origem."
print "1) Escreva dois vetores-diferença usando o ponto-origem  como "
print "origem da diferença e obtendo assim dois lados de um triângulo espacial"
print "que determina o plano;"
print "2) Escreva o determinante funcional usando os dois vetores-diferença como"
print "as duas primeiras linhas e na terceira linha escreva o vetor diferença"
print "(x-a, y-b, z-c)."
print "3) Calcule o determinante para obter a equação do plano."

print "Eu vou resolver com gnuplot a questão 4-a da lista 02."
print "Os demais itens da questão 4  podem ser resolvidos fazendo pequenas
print "modificações neste script." 

print "Aperte enter para continuar!"
pause -2
system "clear"


pow(x,n) = x**n;

F(x,y) = pow(x,2)*sin(x*y)*exp(x*y);
F_x(x,y) = (2*x*sin(x*y) + y*pow(x,2)*cos(x*y)  + y*pow(x,2)*sin(x*y))*exp(x*y);
F_y(x,y) = x*pow(x,2)*cos(x*y)*exp(x*y)+pow(x,2)*sin(x*y)*x)*exp(x*y);

set arrow from 0,0,0 to 1, 1, 1 
set arrow from 0,0,0 to 1, 2 ,3
set arrow from 0,0,0 to 3, -2, 3 

set xrange [-10:10]
set yrange [-10:10]
set zrange [-10:10]


G(x,y) = exp(pow(x,2))*(y+3)*(x+1);
G_x(x,y) = (2*x*(y+3)*(x+1) + y+3 )*exp(pow(x,2)) ;
G_y(x,y) = (x+1)*exp(pow(x,2));
a=1; b = -1; A=G_x(a,b); B = G_y(a,b);
print "vetor perpendicular em (",a,",",b,",",G(a,b),") = "," (",A,",",B,",",-1,") "
quit
splot G(x,y)

pause -2