Реферат: Моделирование ЭВМ

return (l*x)%k;

}

void main()

{ FILE *fout; //выходной файл случайных величин

int i; //параметр цикла

fout=fopen("vi_gpsc1.dat","w");

for(i=1; i<=1000; i++) fprintf(fout,"%f ",float((x=Rnd(x)))/65536);

fclose(fout);

}

Приложение № 3

uses crt;

var f1,f2,f3,f4:text;

i:integer;

x,z1,z2,y1,y2,a,y3,y4,y5,y6:real;

procedure norm1(var x1,x2:real);

begin

z1:=random;

z2:=random;

x1:=sqrt(-2*ln(z1))*cos(2*pi*z2);

x2:=sqrt(-2*ln(z1))*sin(2*pi*z2);

x1:=sqrt(2)*x1+16;

x2:=sqrt(2)*x2+16;

end;

procedure norm2(var x1,x2:real);

begin

z1:=random;

z2:=random;

x1:=sqrt(-2*ln(z1))*cos(2*pi*z2);

x2:=sqrt(-2*ln(z1))*sin(2*pi*z2);