Const x0=1; xk=2; h1=0.1; y0=2; yk=3; h2=0.2; var x,y,z:real; i,j:integer; begin x:=1-h1; for i:=1 to round((xk-x0)/h1)+1 do begin x:=x+h1; y:=2-h2; for j:=1 to round((yk-y0)/h2)+1 do begin y:=y+h2; z:=x*x+y+y; writeln(x:3:1,y:5:1,z:7:3); end; end; end.
// PascalABC.NET 3.3, сборка 1627 от 27.01.2018 // Внимание! Если программа не работает, обновите версию!
begin var n:=4; var xp:=new real[n]; var yp:=new real[n]; Writeln('Введите парами координаты четырех вершин'); for var i:=0 to n-1 do Read(xp[i],yp[i]); var (x,y):=ReadReal2('Введите координаты точки:'); var (i,j,c):=(0,n-1,False); while i<n do begin if ((yp[i]<=y) and (y<yp[j]) or (yp[j]<=y) and (y<yp [i])) and (x<(xp[j]-xp[i])*(y-yp[i])/(yp[j]-yp[i])+xp[i]) then c:=not c; j:=i; i+=1; end; if c then Writeln('Внутри') else Writeln('Не внутри') end.
Примеры Введите парами координаты четырех вершин 1 4 6 9 9 6 4 1 Введите координаты точки: 1.5 1.6 Не внутри
Введите парами координаты четырех вершин 1 4 6 9 9 6 4 1 Введите координаты точки: 4 7 Внутри
Function sq(var x1,y1,x2,y2,x3,y3: real):real; begin // вычисляем площадь треугольника с // векторного произведения sq:=abs((x2-x1)*(y3-y1)-(x3-x1)*(y2-y1))/2; end; // основная программа var x,y,x1,y1,x2,y2,x3,y3,x4,y4,st,sp:real; begin write('координаты 4 вершин x1 y1 x2 y2 x3 y3 x4 y4 через пробел: '); readln(x1,y1,x2,y2,x3,y3,x4,y4); sp:=2*sq(x1,y1,x2,y2,x3,y3); write('координаты точки x y через пробел: '); readln(x,y); st:=sq(x1,y1,x2,y2,x,y) + sq(x2,y2,x3,y3,x,y) + sq(x3,y3,x4,y4,x,y) + sq(x4,y4,x1,y1,x,y); if abs(sp-st)<1.0e-9 then writeln('точка принадлежит прямоугольнику') else writeln('точка не принадлежит прямоугольнику'); end.
координаты 4 вершин x1 y1 x2 y2 x3 y3 x4 y4 через пробел: 1 4 6 9 9 6 4 1 координаты точки x y через пробел: 4 7 точка принадлежит прямоугольнику
координаты 4 вершин x1 y1 x2 y2 x3 y3 x4 y4 через пробел: 1 4 6 9 9 6 4 1 координаты точки x y через пробел: 1.5 1.6 точка не принадлежит прямоугольнику
var x,y,z:real; i,j:integer;
begin
x:=1-h1;
for i:=1 to round((xk-x0)/h1)+1 do
begin
x:=x+h1;
y:=2-h2;
for j:=1 to round((yk-y0)/h2)+1 do
begin
y:=y+h2;
z:=x*x+y+y;
writeln(x:3:1,y:5:1,z:7:3);
end;
end;
end.
Результат:
1.0 2.0 5.000
1.0 2.2 5.400
1.0 2.4 5.800
1.0 2.6 6.200
1.0 2.8 6.600
1.0 3.0 7.000
1.1 2.0 5.210
...
2.0 2.6 9.200
2.0 2.8 9.600
2.0 3.0 10.000