Ввод векторов и матриц, умножение матрицы на вектор и вычисление скалярного призведения двух векторов

[zverushka]

Задача: Даны 6-элементные вещественные векторы х и у и квадратные матрицы А,В,С 6-ого порядка. Надо найти величину (Ах,Ву)+(Сх,у)/(х,Ву)ю Ввод векторов и матриц, умножение матрицы на вектор и вычисление скалярного призведения двух векторов оформить ввиде подпрограмм.

Я не знаю, где ошибка.
program zad6;
uses SysUtils;
const n=6;
type
TMatrix=array[1..n,1..n] of integer;
TVekt=array[1..n] of integer;
var
A,B,C:TMatrix;
X,Y,PAX,PBY,PCX, Res:TVekt;
SP1,SP2,SP3,R: integer;

i,j,sum: integer;

procedure VvodA;
begin
writeln ('Vvod matr A porjdka ',n);
for j:=1 to n do
for i:=1 to n do
begin
write ('A[',i,',',j,'] : ');
readln (A[i,j]);
end;
end; {VvodA}

procedure VvodB;
begin
writeln ('Vvod matr B porjdka ',n);
for j:=1 to n do
for i:=1 to n do
begin
write ('B[',i,',',j,'] : ');
readln (B[i,j]);
end;
end; {VvodB}

procedure VvodC;
begin
writeln ('Vvod matr C porjdka ',n);
for j:=1 to n do
for i:=1 to n do
begin
write ('C[',i,',',j,'] : ');
readln (C[i,j]);
end;
end; {VvodC}

procedure VvodX;
begin
writeln ('Vvod vekt X porjdka ',n);
for i:=1 to n do
begin
write ('X[',i,'] : ');
readln (X[i]);
end;
end; {VvodX}

procedure VvodY;
begin
writeln ('Vvod vekt Y porjdka ',n);
for i:=1 to n do
begin
write ('Y[',i,'] : ');
readln (Y[i]);
end;
end; {VvodY}

procedure ProMV ( M : TMatrix; V : TVekt; var MV: TVekt);
var i,j: integer;
begin
SetLEngth (MV, Legth(V));
for i:=0 to Length(V)-1 do
begin
MV[i]:=0;
for j:=0 to Length(V)-1 do
begin
MV[i]:=MV[i]+M[i,j]*V[j];
end;
end;
end;{ProMV}

function SkalProVekt (var a,b: array of real; n:integer):real;
var i:integer; e: real;
begin
e:=0;
for i:=1 to n-1 do
e:=e+a[i]*b[i];
SkalProVekt:=e;
writeln;
end;{SkalProVekt}

begin
VvodA;
VvodB;
VvodC;
VvodX;
VvodY;

ProMV(A,X,PAX);
ProMV(B,Y,PBY);
ProMV(C,X,PCX);

SP1:=SkalProVekt(PAX,PBY,n);
SP2:=SkalProVekt(PCX,Y,n);
SP3:=SkalProVekt(X,PBY,n);

R:=SP1+SP2/SP3;
writeln ('formula =',R);
end.

[Sazary]

Вам компилятор говорит где ошибка. Как можно не знать где она?
Вот несколько:
в процедуре proMV:

SetLength (MV, Length(V));

SetLength изменяет размер динамического массива, а у вас - статический.
К тому же, еще и опечатка:

Legth(V)

Вот здесь:

SP1:=SkalProVekt(PAX,PBY,n);
SP2:=SkalProVekt(PCX,Y,n);
SP3:=SkalProVekt(X,PBY,n);

Функция принимает динамический массив array of real, а вы ей даете

TVekt=array[1..n] of integer;

Функция SkalProVekt возвращает real, а вы присваиваете ее значение целочисленной переменной:

SP1:=SkalProVekt(PAX,PBY,n);

-----------
Дальше не углублялся.

ps Для кода используйте тег CODE.

[Mikhail Bakurov]

Дофига там еще косяков, в основном все несоответствие типов

заметил кроме вышеперечисленных

SP1,SP2,SP3,R: integer;
................
R:=SP1+SP2/SP3;

Если деление, сразу R:Real;

Да, чуть не забыл, ввод можно сделать одной процедурой для каждого типа

[zverushka]

теперь выдает ошибку, что не может найти sysutils.tpu
Помогите, пожалуйста, с типами
Я не знаю, как сделать ввод одной процедурой для каждого типа

[Sazary]

Подправил самое необходимое. Программа компилируется и даже чего-то выводит.. Проверять, естественно, не стал, т.к. мне стало тоскливо вбивать вручную элементы уже на второй матрице.

program zad6;
uses SysUtils;
const n=6;
type
TMatrix=array[1..n,1..n] of integer;
TVekt=array[1..n] of integer;
var
A,B,C:TMatrix;
X,Y,PAX,PBY,PCX, Res:TVekt;
SP1,SP2,SP3,R: real;

i,j,sum: integer;

procedure VvodABC(var MAS : TMatrix);
begin
for j:=1 to n do
for i:=1 to n do
begin
write ('el[',i,',',j,'] : ');
readln (MAS[i,j]);
end;
end; {VvodA}

procedure VvodX;
begin
writeln ('Vvod vekt X porjdka ',n);
for i:=1 to n do
begin
write ('X[',i,'] : ');
readln (X[i]);
end;
end; {VvodX}

procedure VvodY;
begin
writeln ('Vvod vekt Y porjdka ',n);
for i:=1 to n do
begin
write ('Y[',i,'] : ');
readln (Y[i]);
end;
end; {VvodY}

procedure ProMV ( M : TMatrix; V : TVekt; var MV: TVekt);
var i,j: integer;
begin
for i:=0 to Length(V)-1 do
begin
MV[i]:=0;
for j:=0 to Length(V)-1 do
begin
MV[i]:=MV[i]+M[i,j]*V[j];
end;
end;
end;{ProMV}

function SkalProVekt (var a,b: TVekt; n:integer):real;
var i:integer; e: real;
begin
e:=0;
for i:=1 to n-1 do
e:=e+a[i]*b[i];
SkalProVekt:=e;
writeln;
end;{SkalProVekt}

begin
writeln('matrix A:');
VvodABC(A);
writeln('matrix B:');
VvodABC(B);
writeln('matrix C:');
VvodABC(C);

VvodX;
VvodY;

ProMV(A,X,PAX);
ProMV(B,Y,PBY);
ProMV(C,X,PCX);

SP1:=SkalProVekt(PAX,PBY,n);
SP2:=SkalProVekt(PCX,Y,n);
SP3:=SkalProVekt(X,PBY,n);

R:=SP1+SP2/SP3;
writeln ('formula =',R);
readln;
end.

Пример вывода:

matrix A:
........здесь куча чисел........
Y[6] : 2



formula =-1.40617124792769E+009

ps А не проще случайными заполнять?

[zverushka]

спасибо большое,
пошла проверять

[zverushka]

У меня выдает ошибку
Runtime error 205 at 0000:031D
Подскажите, что не так?

[Sazary]

zverushka
Мало информации.
Какие числа вводите?
Runtime error 205, насколько я помню, возникает при переполнении операций с плавающей точкой.

[zverushka]

я вводила разные,
ошибка даже , если матрицы единичные, а вектора х=1.2.3.4.5.6 и у=1.1.1.1.1.1 ответ должел быть 22

[Sazary]

zverushka Вы так и вводите что-ли, с точками? о_О
Если да, то естественно будет ошибка.
После ввода каждого элемента - <enter>