program vartest;
uses dos,crt;
var tem1,tem2,us,dt,rr,vart1,vart2,varu,tem1a,tem2a,usa:real;
    tem10,tem20,us0:real;
    imax,anz:integer;
    tex:text;
    name:string;
    b1:char;

begin
clrscr;
writeln('Dieses Programm berechnet die Varianz der Temperatur- und');
writeln('Spannungsmeáwerte einer tkm-Datei.');
repeat
write('Name der Meáwertdatei (mit Extension): ');readln(name);
assign(tex,name);
reset(tex);
imax:=-1;vart1:=0;vart2:=0;varu:=0;
repeat
inc(imax);
readln(tex,anz);
readln(tex,tem1);
readln(tex,tem2);
readln(tex,us);
readln(tex,dt);
readln(tex,rr);
if imax>0 then begin
vart1:=vart1+sqr(tem1-tem1a);
vart2:=vart2+sqr(tem2-tem2a);
varu:=varu+sqr(us-usa);
end else begin
tem10:=tem1;tem20:=tem2;us0:=us;
end;
tem1a:=tem1;tem2a:=tem2;usa:=us;
until eof(tex);
close(tex);
vart1:=vart1-sqr(tem10-tem1a)/imax;
vart2:=vart2-sqr(tem20-tem2a)/imax;
varu:=varu-sqr(us0-usa)/imax;
if vart1>0 then
vart1:=sqrt(vart1/2/imax)
else vart1:=0;
if vart2>0 then
vart2:=sqrt(vart2/2/imax)
else vart2:=0;
if varu>0 then
varu:=sqrt(varu/2/imax)
else varu:=0;
writeln(name);
writeln('Varianz von U in æV/K: ',varu:8:4);
writeln('Varianz von T1 in K: ',vart1:8:4);
writeln('Varianz von T2 in K: ',vart2:8:4);
write('Neue Berechnung? (j/n) ');readln(b1);
until b1='n';
writeln;
writeln('Programmende.');
readln;
end.