Нарисовать 3 графика на плоскости координат разными цветами

uses graphABC;
const
x1=-4.5;//начало и конец графика в реальных значениях
x2=1.5;
function Y(x:real):real;
begin
Y:=x*x+3*x-4;
end;
function Z(x:real):real;
begin
Z:=abs(sin(3*x));
end;
function K(x:real):real;
begin
K:=ln(x*x+3*x+1);
end;
 
var m,max,min,x,dx:real;
    x0,y0,i,n:integer;
    s:string;
begin
setwindowsize(500,500);
x0:=round(-(windowwidth-40)*x1/(x2-x1));//начало координат по оси Х
max:=Y(x1); //найдем макс и мин для масштаба по оси У для параболы,
min:=Y(x1); // синус от 0 до 1, логарифм от минус бесконечности до плюс бесконечности
dx:=0.05;
x:=x1;
while x<x2 do
 begin
  if Y(x)>max then max:=Y(x)
  else if Y(x)<min then min:=Y(x);
  x:=x+dx;
 end;
y0:=round((windowwidth-40)*max/(max-min));//начало координат по оси Y
if max-min>x2-x1  then m:=(windowwidth-40)/(max-min)//масштаб
else m:=(windowwidth-40)/(x2-x1);
n:=trunc(x2)+1;//количество делений по осям в 1 сторону
if -min>x2 then n:=trunc(-min)+1;
if max>-min then n:=trunc(max)+1;
line(20,y0,windowwidth-20,y0); //ось Х
line(x0,20,x0,windowheight-20); //ось У
for i:=1 to n do
 begin
  line(x0-3,y0-round(i*m),x0+3,y0-round(i*m));//засечки на оси У
  line(x0-3,y0+round(i*m),x0+3,y0+round(i*m));
  line(x0+round(i*m),y0+3,x0+round(i*m),Y0-3); //засечки на оси Х
  line(x0-round(i*m),y0+3,x0-round(i*m),Y0-3);
  str(i,s);
  //подпись оси У
  textout(x0-20,y0-round(i*m),s);{соответственно засечкам}
  textout(x0-25,y0+round(i*m),'-'+s);
  //подпись оси Х
  textout(x0+round(i*m),y0+10,s);
  textout(x0-round(i*m),y0+10,'-'+s);
  end;
//центр
textout(x0+5,y0+10,'0');
//подписи концов осей
textout(windowwidth-30,y0-20,'X');
textout(x0+10,10, 'Y');
//график
x:=x1;
dx:=0.001;
while x<=x2 do
 begin
  setpixel(x0+round(x*m),y0-round(Y(x)*m),clRed);
  setpixel(x0+round(x*m),y0-round(Z(x)*m),clBlue);
  if x*x+3*x+1>0 then setpixel(x0+round(x*m),y0-round(K(x)*m),clGreen);
  x:=x+dx; //наращиваем х
 end;
textout(30,y0+100,'Графики функций');
setfontcolor(clRed);
textout(30,y0+120,'Y=x^2+3x-4');
setfontcolor(clBlue);
textout(30,y0+140,'Z=|sin3x|');
setfontcolor(clGreen);
textout(30,y0+160,'Y=ln(x^2+3x+1)');
end.