% graf_2var.m

% script en octave

% representa el graf d'una funcio en 2 variables sobre un rectangle

% Jaume Amoros, UPC, Barcelona
% 2003/11/11

% Cal posar aquest programa en el directori de treball d'Octave.
% Per Octave 2.1.36 en Windows 98 es 
% C:\Archivos de programa\GNU Octave 2.1.36\octave_files

% AVIS: aquest es el programa que dibuixa la funcio i els eixos
% amb el llistat mes facil de seguir, pero NO es el mes rapid.
% Podeu fer-lo mes rapid si calculeu la funcio amb operacions vectorials.

% inici del programa

% declara el rectangle [x0,xf]x[y0,yf] de dibuix del grafic
% i el nombre de punts de calcul de la funcio en cada eix nx,ny
% (escriure'ls en el propi llistat)


x0=-2;
xf=1;
y0=-1;
yf=1;

nx=21;
ny=21;

% declara els vectors de valors de x,y

x=linspace(x0,xf,nx);
y=linspace(y0,yf,ny);


% CALCUL PUNT A PUNT
% calcula la funcio en els punts (mes lent que fer-ho vectorialment!)

for j=1:ny,
  for i=1:nx,

% poseu aqui la vostra funcio
% es pot definir en mes d'una linea

    f(i,j)=x(i)^2-y(j)^2;
  end;
end;

% prepara els eixos pel dibuix:

% preparem les matrius de valors de x i de y en els punts de la malla rectangular

[xx,yy]=meshgrid(x,y);

% ajustem z

z0=min(min(f));
zf=max(max(f));

% Posa etiquetes en els eixos:

xlabel('eix x');
ylabel('eix y');
zlabel('z=f(x,y)');

% dibuixa la funcio

mesh(x,y,f);

% s'obren dues finestres del programa grafic gnuplot
% una de comandes i una amb la grafica.

% per variar el punt de vista sobre la grafica:
% en la finestra de comandes gnuplot que s'ha obert picar dues ordres, primer fer

% set view longitud,latitud

% (latitud i longitud en graus sexagesimals) i despres

% replot

% i el dibuix apareixera des del nou angle de vista