Цитата:
Сообщение от Dima007
Не могли бы Вы помочь мне с этим разобраться?)
|
Да куда я денусь...
Код:
Type
Vect=array[1..3] of real;
Matr=Array[1..3,1..3] of real;
Const
z:Vect=( -2,3, 2);
y:Vect=(0.5,3,-2);
a:Matr=((2,4,6),(3,5,2),(2,4,7));
Function Prod(PM:Matr;PV:Vect):Vect;
var
i,k:Byte;
M:Vect;
begin
for i:=1 to 3 do M[i]:=0;
for i:=1 to 3 do
for k:=1 to 3 do
M[i]:=M[i]+PM[i,k]*PV[k];
Prod:=M;
end;
Function Sums(V1,V2:Vect):Vect;
var
i:Byte;
begin
for i:=1 to 3 do Sums[i]:=V1[i]+V2[i];
end;
Function Scal_Prod(V1,V2:Vect):Real;
var
i:Byte;
SP:Real;
begin
SP:=0;
for i:=1 to 3 do SP:=SP+V1[i]*V2[i];
Scal_Prod:=SP;
end;
Var
k:Real;
Begin
k:=Scal_Prod(Sums(Prod(A,y),z),Prod(A,y));
Writeln('k = ',k:0:3);
Readln
End.