Курсовая работа: Вероятностные процессы и математическая статистика в автоматизированных системах

type arr=array [1..3,1..3] of real;

var d,d1,d2,d3:arr; i:integer; bd:mas;

function delta(a:arr):real;

begin delta:=a[1,1]*(a[2,2]*a[3,3]-a[2,3]*a[3,2])-

a[1,2]*(a[2,1]*a[3,3]-a[2,3]*a[3,1])+

a[1,3]*(a[2,1]*a[3,2]-a[2,2]*a[3,1]); end;

begin d[1,1]:=2*b[7]; d[1,2]:=b[4]; d[1,3]:=b[5];

d[2,1]:=b[4]; d[2,2]:=2*b[8]; d[2,3]:=b[6];

d[3,1]:=b[5]; d[3,2]:=b[6]; d[3,3]:=2*b[9];

bd[1]:=-b[1]; bd[2]:=-b[2]; bd[3]:=-b[3];

for i:=1 to 3 do

begin d1[i,2]:=d[i,2]; d1[i,3]:=d[i,3];

d2[i,1]:=d[i,1]; d2[i,3]:=d[i,3];

d3[i,2]:=d[i,2]; d3[i,1]:=d[i,1];

d1[i,1]:=bd[i];d2[i,2]:=bd[i];d3[i,3]:=bd[i]; end;

xc[1]:=delta(d1)/delta(d);

xc[2]:=delta(d2)/delta(d);

xc[3]:=delta(d3)/delta(d);

ycen:=fy(xc[1],xc[2],xc[3]);

la[1]:=findl(-6,-2);

la[2]:=findl(-2,2);

la[3]:=findl(2,6);

for i:=1 to 3 do begin

ml[i]:=(b[6]*(b[7]-la[i])-0.5*b[4]*b[5])/(b[5]*(b[8]-la[i])-0.5*b[4]*b[6]);

nl[i]:=(-2*(b[7]-la[i])*(b[8]-la[i])+0.5*b[4]*b[4])/(b[5]*(b[8]-la[i])-0.5*b[4]*b[6]);

l[i]:=1/sqrt(1+sqr(ml[i])+sqr(nl[i]));

m[i]:=ml[i]*l[i];

n[i]:=nl[i]*l[i];

maximminim; end; end;