OberonCore

   PROCEDURE SetPoint(VAR ap: ARRAY OF Ports.Point; x,y : INTEGER; n: INTEGER);
   BEGIN
         ap[n].x := x; ap[n].y := y;
   END SetPoint;

   PROCEDURE Draw*;
   VAR p: ARRAY 10 OF Ports.Point;
   BEGIN
       SternCanvas.SetBgColor(100,100,1);
      (*подготовка массива точек*)
      SetPoint(p, 10,5, 0);
      SetPoint(p, 40,20, 1);
      SetPoint(p, 80,10, 2);
      SetPoint(p, 160,20, 3);
      SetPoint(p, 200,30, 4);
      SetPoint(p, 180,80, 5);
      SetPoint(p, 120,100, 6);
      SetPoint(p, 80,110, 7);
      SetPoint(p, 40,50, 8);
      SetPoint(p, 10,30, 9);
      SternCanvas.Path(p, 0, SternCanvas.red, SternCanvas.openBezier);
      SternCanvas.Path(p, 0, SternCanvas.blue, SternCanvas.openPoly);
   END Draw;

MODULE  СборGrm10; 
   IMPORT  Gr := Info21sysTPGraphics, Math;

   PROCEDURE Окр( x: REAL): INTEGER;
   BEGIN
      RETURN SHORT( ENTIER( x + 0.5 ) )
   END Окр;
         
   PROCEDURE Do*;
      VAR  i, k, x, y: INTEGER;
   BEGIN Gr.Clear;   
      Gr.SetColor( 0 ); Gr.SetBkColor( 15 );
      Gr.MoveTo( 0, 70 );
      FOR i := 0 TO 300 DO
         y := 70 - Окр( 15 * Math.Sin( i / 3 ) +
            40 * Math.Sin( i / 30 ) );
         x := i; Gr.LineTo( x, y )
      END;
      Gr.MoveTo( 0, 180 );
      FOR i := 0 TO 300 DO
         y := 180 - Окр( 10 * Math.Sin( i / 1.3 ) +
            30 * Math.Sin( i / 16 ) );
         x := i; Gr.LineTo( x, y )
      END;
      Gr.String(80, 15, "Сложение синусоид");
      Gr.Open
   END Do;

END СборGrm10.

(Q) СборGrm10.Do 

   Здесь показано сложение синусоид с разными периодами и амплитудами.

MODULE SternObx1;

   IMPORT SternCanvas;

   (* y=d*x^2+f  *)
   (* начало координат располагается в левом верхнем углу *)
   
   CONST
      d = 1; f = 1;
      scale = 500; (* масштабный коэффициент для ординаты *)
   
   PROCEDURE Y (x: INTEGER) : INTEGER;
   BEGIN
      RETURN d*(x*x)+f
   END Y;

   PROCEDURE SetPoint(VAR ap: ARRAY OF SternCanvas.Point; x: INTEGER; n: INTEGER);
   BEGIN
         ap[n].x := x; ap[n].y := Y(x) DIV scale;
   END SetPoint;

   PROCEDURE Draw*;
   VAR p: ARRAY 10 OF SternCanvas.Point;
   BEGIN
       SternCanvas.SetBgColor(200,200,200);
      (*подготовка массива точек*)
      SetPoint(p, 20, 0);
      SetPoint(p, 40, 1);
      SetPoint(p, 60, 2);
      SetPoint(p, 80, 3);
      SetPoint(p, 100, 4);
      SetPoint(p, 120, 5);
      SetPoint(p, 140, 6);
      SetPoint(p, 160, 7);
      SetPoint(p, 180, 8);
      SetPoint(p, 200, 9);
      SternCanvas.Path(p, 0, SternCanvas.blue, SternCanvas.openPoly);
   END Draw;

END SternObx1.

MODULE TestPlot;

   IMPORT LibPlotters, Ports, Log := StdLog;

      (*  y = a*x^2 + b*x + c  *)

   TYPE
      Plot = LibPlotters.Plotter;

   CONST
      SHIFT = 0.1; (*шаг смещения переменной X*)

   PROCEDURE DrawPlot (a, b, c, minX, maxX: INTEGER);
      VAR
         plot: Plot;
         minY, maxY, X, Y, shiftX, shiftY: REAL;
         n: INTEGER;
   BEGIN

      (*создаем новый график и открываем окно с ним*)
      plot := LibPlotters.dir.NewPlotter('');
      LibPlotters.OpenAux(plot, 150, 100, TRUE, 'График функции y = a*x^2 + b*x + c ');

      (*вычисляем максимальное и минимальное значения шкалы Y*)
      IF a > 0 THEN
         minY := - b / (2 * a);
         maxY := a * maxX * maxX + b * maxX + c;
      ELSE
         maxY := - b / (2 * a);
         minY := a * maxX * maxX + b * maxX + c
      END;

      plot.Scale(minX - 1, maxX + 1, minY - 1, maxY + 1);

      (*определяем размеры ячейки сетки координат*)
      shiftX := ENTIER((maxX - minX) / 10);
      shiftY := ENTIER((maxY - minY) / 10);

      (*рисуем сетку координат*)
      plot.Grid(minX - shiftX, maxX + shiftX, shiftX, minY - shiftY, maxY + shiftY, shiftY, - 2, Ports.grey12, FALSE);

      (*рисуем шкалу X*)
      plot.Xaxis(minX, maxX, 1, 5, 2, 0, {});

      (*рисуем шкалу Y*)
      plot.Yaxis(minY, maxY, 1, 5, 0, 0, TRUE, {});

      (*добавляем точки на график и их координаты значения в StdLog*)
      X := minX;
      Y := a * X * X + b * X + c;
      Log.Clear;

      WHILE X <= maxX DO
         Log.Ln;
         Log.RealForm(X, 5, 2, 0, 8FX); Log.String("  "); Log.RealForm(Y, 5, 2, 0, 8FX);
         plot.Point(X, Y, 0, 0, 1, 0, 1, Ports.blue);
         X := X + SHIFT;
         Y := a * (X * X) + b * X + c
      END;

      (*обнавляем график*)
      plot.Update
   END DrawPlot;

   PROCEDURE Do*;
   BEGIN
(*Нарисуем несколько графиков*)
      DrawPlot(1, 2, 0, - 14, 10);
      DrawPlot(-2, 3, 1, - 14, 10);
      DrawPlot(4, 5, 2, -5, 5);
   END Do;

END TestPlot.

 TestPlot.Do

from numpy import *
from pylab import *

x = arange(0, 10, 0.01)
y = abs(sin(x))

plot(x, y)
show()