Построить график изотермы

program pro3;
uses graphABC;
var a,b,T,xn,xk,x,mx,my,dx:real;
    x0,y0,i,V1,V2,V3,V4,V5,V6,V7,V8:integer;
    s:string;
    V,P: array[1..8]of real;
 
begin
x0:=windowwidth div 2;{Центр экрана}
y0:=windowheight div 2;
xn:=-10;xk:=10;{Интервал по Х}
mx:=(x0-30)/xk;{Мастштаб по У}
my:=y0/10;{По У}
line(0,y0,windowwidth,y0);{оси}
line(x0,0,X0,windowheight);
for i:=1 to 10 do{Максимальное колво засечек на одну сторону}
  begin
    line(x0+round(i*mx),y0-3,x0+round(i*mx),Y0+3); {Засечки на оси Х}
    line(x0-round(i*mx),y0-3,x0-round(i*mx),Y0+3);
    line(x0+3,y0-round(i*my),x0-3,y0-round(i*my)); {Засечки на осиY}
    line(x0+3,y0+round(i*my),x0-3,y0+round(i*my));
    str(i,s);
    {подпись оси Х}
    textout(x0+round(i*mx),y0+10,s);
    textout(x0-round(i*mx),y0+10,'-'+s);
    {подпись оси Y}
    textout(x0-25,y0-round(i*my),s);
    textout(x0-25,y0+round(i*my),'-'+s);
  end;
{Центр}
textout(x0+5,y0+10,'0');
{Подписи концов осей}
textout(windowwidth-10,y0-10,'X');
textout(x0+5,10, 'Y');
{график}
x:=xn;
dx:=0.001;
begin
writeln('Введите объем газа V');
 for i:=1 to 8 do
begin
readln(V[i]);
 end;
  writeln('введите  a');
  readln(a);
   writeln('введите  b');
   readln(b);
    writeln('Введите температуру газа T');
    readln(T);
for i:=1 to 8 do
begin
P[i]:=((0.0635*T)/(V[i]-b))-((a/Power(V[i],2)));{Вычисление значения p}
writeln('Уравнение состояния расширения 1 кмоля углекислого газа = ',P[i]);
end;
end;
end.

uses graphABC;
const n=8;
type mas=array[1..n] of real;
const a=0.0288;
      b=0.032;
      t=15.7;
      v:mas=(0.05,0.1,0.15,0.2,0.4,0.6,0.8,1);
function f(vi:real ):real;
begin
f:=0.0635*t/(vi-b)-a/sqr(vi);
end;
var p:mas;
    x0,y0,m,i:integer;
    mx,my,max:real;
    s:string;
begin
//заполнение таблицы значений
max:=0;
for i:=1 to n do
 begin
  p[i]:=f(v[i]);
  if p[i]>max then max:=p[i];
 end;
setwindowsize(800,400);
centerwindow;
x0:=windowwidth div 4;//начало координат
y0:=windowheight-50;
m:=20;
mx:=(windowwidth-x0-20)/m; //масштабы по осям
my:=(y0-60)/max;
//координатная сетка
line(x0,y0,windowwidth,y0);
textout(windowwidth-20,y0-20,'X');
line(x0,y0,x0,20);
textout(x0+10,20,'Y');
for i:=0 to m do
 begin
  line(x0+round(mx*i),y0+3,x0+round(mx*i),y0-3);
  str(i*0.05:0:2,s);
  textout(x0+round(mx*i),y0+10,s);
  line(x0+3,y0-round(my*i*10),x0-3,y0-round(my*i*10));
  textout(x0-20,y0-round(my*i*10)-5,inttostr(i*10));
 end;
//график
setpencolor(clBlue);
setpenwidth(2);
circle(x0+round(v[1]*mx*20),y0-round(p[1]*my),2);
moveto(x0+round(v[1]*mx*20),y0-round(p[1]*my));
for i:=1 to n do
 begin
  setpencolor(clRed);
  lineto(x0+round(v[i]*mx*20),y0-round(p[i]*my));
  setpencolor(clBlue);
  circle(x0+round(v[i]*mx*20),y0-round(p[i]*my),2);
 end;
{построение таблицы}
setpencolor(clBlack);
setpenwidth(1);
textout(40,0,'Значения аргумента');
textout(60,20,'и функции');
line(40,50,40,50+(n+1)*18);
line(90,50,90,50+(n+1)*18);
line(140,50,140,50+(n+1)*18);
for i:=0 to n+1 do
line(40,50+i*18,140,50+i*18);
textout(65,55,'X');
textout(115,55,'Y');
for i:=1 to n do
 begin
  str(v[i]:0:2,s);
  textout(45,55+18*i,s);
  str(p[i]:0:3,s);
  textout(95,55+18*i,s);
 end;
end.