След матрицы

[NauDS]

Даны квадратная матрица порядка m и натуральное число n. Вычислить следы матрицы А, А^2,..., А^n. Помогите, с процедурами и функциями не работал, а срочно нужно. язык паскаль.

[Sazary]

Смотрим вниз страницы на "Похожие темы" и видим...ничего знакомого?
Пользуйтесь поиском.
http://programmersforum.ru/showthread.php?t=39812

[NauDS]

Там решение на си, а для меня это еще больший лес

[Sazary]

А мою вторую фразу вы не увидели? Или не смотрели?
Есть поиск. Запрос "След матрицы". Одна из тем - http://programmersforum.ru/showthrea...F2%F0%E8%F6%FB
Для нахождения A^2 и т.д. введите "Умножение матриц".
Обсуждалось много раз.

[NauDS]

Вот моя программа, видимо я где-то с умножением ошибся, завтра сдавать, а я никак не могу разобраться до конца.
program sledmatr;
const
Nmax=10;
Type
tmatr=array[1..nmax,1..nmax] of real;
var
s:real;
Na,ia,ja:byte;
A,C:tmatr;
Procedure EnterMatr( var matr:tmatr; var n:byte; ch:char);
var
i,j:byte;
Begin
Writeln ('Введите размер матрицы ', ch);
write ('Количество строк=');
readln(n);
write ('Количество столбцов=');
readln(n);
for i:=1 to n do
for j:=1 to n do
begin
write(ch,'[',i,',',j,']=');
readln (matr[i,j]);
end;
End;
Procedure Ymn(a:tmatr; n:byte; var C:tmatr);
Var
i,j,k:byte;
s:real;
Begin
begin
for i:=1 to n do
for j:=1 to n do
c[i,j]:=0;
for k:=1 to n do
c[i,j]:=c[i,j]+a[i,k]*a[k,j];
End;
End;
Function Sled (const matr:tmatr;n:byte):integer;
var
i,j:byte;
s:real;
begin
s:=0;
for i:=1 to n do
s:=s+matr[i,i];
writeln(s:6:2);
End;
Procedure PrintMatr (const matr:tmatr; N:byte; ch:char);
var
i,j:byte;
s:real;
Begin
writeln('Матрица ',ch);
for i:=1 to n do
begin
for j:=1 to n do
write(matr[i,j]:6);
writeln;
end;
end;
Begin
EnterMatr(A,Na,'A');
Sled(A,Na);
Ymn(A,Na,C);
PrintMatr (A,Na,'A');
End.

[Sazary]

В функции умножения двойной begin:

Begin
begin

и с end'ом то же самое.
Про то, что матрицу нужно возвести в степень n (которую, кстати, надо ввести) вы, видимо, забыли.

program sledmatr;
const
Nmax=10;
Type
tmatr=array[1..nmax,1..nmax] of real;

var
Na,ia,ja:byte;
A,C:tmatr;
N,i : integer;

Procedure EnterMatr( var matr:tmatr; var n:byte; ch:char);
 var
 i,j:byte;
 Begin
 Writeln ('Введите размер матрицы ', ch);
 write ('Количество строк=');
 readln(n);
 write ('Количество столбцов=');
 readln(n);
 for i:=1 to n do
  for j:=1 to n do
   begin
   write(ch,'[',i,',',j,']=');
   readln (matr[i,j]);
   end;
 End;

Procedure Ymn(a:tmatr; n:byte; var C:tmatr);
 Var
 i,j,k:byte;
 Begin
 for i:=1 to n do
  for j:=1 to n do
   begin
   c[i,j]:=0;
   for k:=1 to n do
    c[i,j]:=c[i,j]+a[i,k]*a[k,j];
   End;
 End;

Function Sled (const matr:tmatr;n:byte):integer;
 var
 i,j:byte;
 s:real;
 begin
 s:=0;
 for i:=1 to n do
  s:=s+matr[i,i];
 writeln('Sled= ',s:6:2);
 End;

Procedure PrintMatr (const matr:tmatr; N:byte; ch:char);
 var
 i,j:byte;
 s:real;
 Begin
 writeln('Матрица ',ch);
 for i:=1 to n do
  begin
  for j:=1 to n do
   write(matr[i,j]:7:0);
  writeln;
  end;
 writeln;
 end;

Begin
EnterMatr(A,Na,'A');
write('Vvedite N: ');
readln(N);
Sled(A,Na);
for i:=1 to n-1 do
 begin
 writeln;
 Ymn(A,Na,C);
 A := C;
 PrintMatr (A,Na,'A');
 Sled(A,Na);
 end;
readln;
End.