//dessin d'une fonction bi-dimensionnelle
#include <iostream>
#include <cmath>
#include "dislin.h"     //bibliotheque DISLIN

using namespace std;

const double DEG2RAD = M_PI/180.;

int main() {
  const int POINTSX = 19, POINTSY = 19;  
  double z[POINTSX][POINTSX];  
  double x, y;
  double dx = 20., dy = 20.;
  //tabulation de la fonction
  for (int i=0; i<POINTSX; i++)
    for (int j=0; j<POINTSY; j++) {
      x = i*dx * DEG2RAD;
      y = j*dy * DEG2RAD;
      z[i][j] = 2. * sin(x) * cos(y);
    } 
  
  //initialisation de DISLIN
  metafl("PDF");   //sortie fichier PDF
  setfil("fxy.pdf");   //nom du fichier
  filmod("VERSION");
  //metafl("XWIN");   //sortie sur l'ecran
  disini();
  name("axe-X","X");
  name("axe-Y","Y");
  name("axe-Z","Z");
  titlin("f(x,y) = 2 * sin(x) * cos(y)",4);

  //dessin du graphique
  graf3d(0.,360.,0.,90.,0.,360.,0.,90.,-3.,3.,-3.,1.);
  title();
  
  color("GREEN");
  //passage du tableau z par pointeur
  //&z[0][0] donne l'adresse memoire du premier element du tableau
  surmat(&z[0][0],POINTSX,POINTSY,1,1);
  
  disfin();
 
  return 0;
}