%% reglage du format de l'image
300 setwidth
150 setheight

%% taille de la fenetre
-1 12 setxrange
-1 24 setyrange

%% trace du repere
traceaxes

/f {
   3 mul 4 sub
} def


%% on nomme le tableau de points
/les_points {
   %% le tableau des x
   [1 2 3 4 5 6 7 8 9]
   %% le tableau des y
   [0 2 3 10 14 13 13 20 23]
   fuz	   %% fusion des tableaux pour obtenir un tableau de points
} def

les_points points	   %% on represente les points
{f} courbe		   %% on trace la droite de regression

%% syntaxe : A D D' paralleloproj --> A', le projete de A sur D'
%% parallelement a D
/paralleloproj {
3 dict begin
   /d2 defdroite
   /d1 defdroite
   /A defpoint
   d2
   d1 A paral
   interdroite
end
} def

les_points {OX 0 4 -1 1 paralleloproj} papply
{times2} plot