avaliações numéricas de chaminés de equilíbrio - ppgerha ...

for i := 1 to 2 do
begin
v[i,1]^ := v[i,2]^;
h[i,1]^ := h[i,2]^;
end;
lax;
reservatorio;
chamine;
valvula;
idt := idt + 1;
if idt = ndt
then
begin
idt := 0;
writeln(t,' ',ta);
writeln(arq,t:10:5,' ',h[1,2]^[0],' ',hch[2],
' ',h[2,2]^[nx[2]]);
end;
t := t + dt;
until t > ta;
close(arq);
end.
146