#include <iostream>
#include "dislin.h"     //bibliotheque DISLIN

using namespace std;

//declaration la fonction a dessiner
double fonct(double x);

int main() {
  cout << "Dessin d'une fonction" << endl;

  //tabulation de la faonction  
  const int POINTS = 100;
  double x[POINTS+1];
  double y[POINTS+1];
  double dx = 1.;
  //boucle pour remplir les tableaux x et f(x)
  for (int i=0; i<=POINTS; i++) {
    x[i] = i * dx;
    y[i] = fonct(x[i]); 
  }
  
  //affichage sur l'ecran
  metafl("XWIN");
  qplot(x,y,POINTS);
 
  return 0;
}

//definition de la fonction a dessiner
double fonct(double x) {
  double y = 1./(1.+x*x);
  return y;       
}