Заданы уравнения трех кривых на плоскости. Pascal

[LiVilman]

Заданы уравнения трех кривых на плоскости. Одно из уравнений содержит переменный коэффициент k. Кривые, пересекаясь, ограничивают замкнутую область D, размеры которой зависят от значения коэффициента k. задано начальное значение k нач , конечное значение k кон и шаг изменения k;

На той же плоскости заданы декартовы координаты точек, образующих множество М (координаты задаются произвольно, n не больше 20).

Для каждого из заданных значений k сформировать P - подмножество M, содержащее точки, которые находятся внутри замкнутой области D.



Уравниния линий:
y2=0,75x
y=ln(x-0,2)
y=0,2x-k

k нач 0,1
k кон 0,5
шаг изменения k 0,1

Прошу помощи в решении. Что то пытался сделать сам но не до конца.

PHP код:

uses crt, graph;

type
TPoint = record x,y : Real; end;

var
Pts : array [1..20] of TPoint;
drv, mode      : Integer;
RangeX, RangeY : Real;
ch, ch1        : Char;
ks             : Array[1..10] of Real;
curK           : Integer;

procedure InitK;
var x : Integer;
begin
for x := 1 to 10 do ks[x] := Random * (2.2 - 1) + 1;
end;

procedure InitPoints;
var x : Integer;
begin
for x := 1 to 20 do
        begin
        pts[x].x := RangeX*random/2;
        pts[x].y := RangeY*random/2;
        end;
end;

Function gFunc1(x:Real):Real;
begin
gFunc1 := x*x - x + 1;
end;

Function gFunc2(x, k:Real):Real;
begin
gFunc2 := k * cos(x + 0.25);
end;

Function gFunc3(x:Real):Real;
begin
if x>=0 then gFunc3 := sqrt(1.5*x)
    else gFunc3 := 0;
end;

procedure PointToScreen(X,Y:Real;color : Integer);
begin
PutPixel(trunc((X+RangeX/2)/RangeX*GetMaxX), trunc((RangeY/2-Y)/RangeY*GetMaxY), color);
end;

procedure Solve;
var x    : Integer;
    a, b : Integer;
begin
For x := 1 to 20 do
    begin
    if (gFunc1(pts[x].X) < pts[x].Y)
    and (gFunc2(pts[x].X, ks[CurK]) < pts[x].Y)
    and (gFunc3(pts[x].X) > pts[x].Y)
    and (pts[x].X >= 0) then
        begin
        a := trunc((pts[x].X+RangeX/2)/RangeX*GetMaxX);
        b := trunc((RangeY/2-pts[x].Y)/RangeY*GetMaxY);
        SetColor(13);
        Circle(a, b, 4);
        end;
    end;
end;

procedure Repaint;
var x   : Integer;
    r,a : Real;
begin
ClearViewPort;
SetColor(11);
Moveto(0, GetMaxY div 2);
LineTo(GetMaxX, GetMaxY div 2);
MoveTo(GetMaxX div 2, 0);
LineTo(GetMaxX div 2, GetMaxY);
for x := 0 to GetMaxx-1 do
    begin
    a := X/GetMaxX*RangeX - RangeX/2;
    PutPixel(x, trunc((RangeY/2 - gFunc1(a))/RangeY * GetMaxY), 15);
    PutPixel(x, trunc((RangeY/2 - gFunc2(a, ks[CurK]))/RangeY * GetMaxY), 15);
    PutPixel(x, trunc((RangeY/2 - gFunc3(a))/RangeY * GetMaxY), 15);
    end;
For x := 1 to 20 do
    PointToScreen(pts[x].X, pts[x].Y, 14);
Solve;
end;

begin
Randomize;
CurK := 1;
RangeX := 6;
RangeY := 4;
drv := detect;
mode := 2;
InitGraph(drv, mode, '');
InitPoints;
InitK;
Repaint;
repeat
ch := UpCase(Readkey);
ch1:= #0;
if ch=#0 then ch1 := Readkey;
if ch1='H' then
        begin
        RangeX := RangeX/1.1;
        RangeY := RangeY/1.1;
        Repaint;
        end;
if ch1='P' then
        begin
        RangeX := RangeX*1.1;
        RangeY := RangeY*1.1;
        Repaint;
        end;
if ch ='A' then
        begin
        if CurK < 10 then inc(CurK);
        Repaint;
        end;
if ch ='Z' then
        begin
        if CurK > 1 then dec(CurK);
        Repaint;
        end;
if ch ='Q' then
        begin
        InitPoints;
        Repaint;
        end;
until (ch=#27) or (ch=#13);
CloseGraph;
end.