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

using namespace std;

int main() {
  const int POINTSX = 21, POINTSY = 21;  
  double z[POINTSX][POINTSX];  
  double x, y;
  double dx = 1, dy = 1;
  //tabulation de la fonction
  for (int i=0; i<POINTSX; i++)
    for (int j=0; j<POINTSY; j++) {
      x = -10. + i*dx;
      y = -10. + j*dy;
      z[i][j] = x*x + y*y;
    } 
  
  metafl("XWIN");
  disini();
  //dessin du graphique
//  titlin("f(x,y) = x*x + y*y",4);
  graf3d(-10.,10.,-10.,5.,-10.,10.,-10.,5.,0.,250.,0.,25.);
//  title();
  color("GREEN");

  //passage du tableau z par pointeur
  //&z[0][0] donne l'adresse memoire du premier element du tableau
  qplsur(&z[0][0],POINTSX,POINTSY);

  disfin();
   
  return 0;
}