Добрый день.
Помогите исправить код для паскаля - тут по мелочи, но я не могу понять что имено исправить. Код вродебы написан правильно, но паскаль придирается к точкам с запятой и т.д.
uses crt;
const mf=500;
type
vector=array[l..mf] of real;
var
i, j, N, N1, N2, N3 : integer;
T, alfa, beta :vector;
ai, bi, ci, fi : real;
a, lamda, ro, c : real;
h, tau, t_end, time : real;
T0, L, kapa, Te, q :real;
f, tg : text;
begin
clrscr;
Writeln('Введите количество пространственных узлов в пластине, N’);
Readln(N).
Writeln('Введите окончание по времени, t_end’);
Readln(t_end).
Writeln('Введите толщину пластины, L’);
Readln(L).
Writeln('Введите коэффициент теплопроводности материала пластины,
lamda');
Readln(lamda).
Writeln('Введите плотность материала пластины, ro’);
Readln(ro);
Writeln('Введите теплоемкость материала пластины, c’);
Readln(c);
Writeln('Введите температуру окружающей среды, Te’);
Readln(Te);
Writeln('Введите коэффициент теплообмена, kapa’);
Readln(kapa);
Writeln('Введите начальную температуру, T0’);
Readln(T0);
Writeln('Введите мощность внутренних источников тепла, q’); Readln(q);
Nl:= N div 4;
N2:=2*N1;
N3:=3*N1;
h:=L/(N-l);
a:=lamda/(ro*c);
tau:=t_end/100.0;
for i:= 1 to N do
T[i]:=T0;
time:=0;
while time<t_end do
begin
time:=time+tau;
alfa[l]:=2.0*a*tau*lamda/(lamda*sqr(h)+2.0*a*tau*(lamda+kapa*h)); beta[l]:=(lamda*sqr(h)*T[1]+2.0*a*tau*kapa*h*Te)/(lamda*sqr(h)+2.0*a*tau*(lamda+kapa*h));
for i:= 2 to N-1 do
begin
ai:=lamda/sqr(h);
bi:=2.0*lamda/sqr(h)+ro*c/tau;
ci:=lamda/sqr(h):
fi:=-ro*c*T[i]/tau;
if (i=N1) or (i=N2) then fi:=-ro*c*T[i]/tau-h*(i-l)*q;
if i=N3 then fi:=-ro*c*T[i]/tau-h*(N1-1)*q;
alfa[i]:=ai/(bi-ci*alfa[i-1];
beta[i] :=(ci*beta [i-1]-fi)/(bi-ci*alfa[i-1]);
end;
T[N]:=(lamda*sqr(h)*T[N]+2.0*a*tau*(lamda*beta[N-l]+kapa*h*Te))/(lamda*sqr(h)+2.0*a*tau*(lamda*(l-alfa[N-1])+kapa*h));
for i:=N-l downto 1 do
T[i]:=alfa[i]*T[i+l]+beta[i];
end;
Assign(f,'res.txt');
Rewrite(f);
Writeln(f,'Толщина пластины L=’,L:6:4);
Writeln(f,'Число узлов по пространственной координате в пластине N = ‘,N);
Writeln(f,'Внутренние источники находятся в точках = ‘,N1,’ ‘,N2,’ ‘,N3); Writeln(f,’Внутренние источники находятся в точках = ',(N1-1)*h:6:4,' ‘,(N2-l)*h:6:4,’ ‘,(N3-l)*h:6:4);
Writeln(f,'Коэффициент теплопроводности материала пластины lamda = ',lamda:6:4);
Writeln(f,'Плотность материала пластины ro = ',ro:6:4);
Writeln(f,'Теплоемкость материала пластины c = ‘,c:6:4);
Writeln(f,'Начальная температура T0 = ‘,T0:6:4);
Writeln(f,'Коэффициент теплообмена kapa = ‘,kapa:6:4);
Writeln(f,'Температура окружающей среды Te = ',Te:6:4);
Writeln(f,'Составляющая мощности внутренних источников тепла q = ‘,q:6:4);
Writeln(f,'Результат получен с шагом по координате h = ',h:6:4);
Writeln(f,'Результат получен с шагом по времени tau = ',tau:6:4);
Writeln(f,'Температурное поле в момент времени t = ',t_end:6:4);
close(f);
Assign(g,'tempr.txt’);
Rewrite(g);
for i:=l to N do
Writeln(g,' ',h*(i-1):10:8,' ',T[i]:8:5);
close(g);
end.