Помогите составить программу вычисления матричного синуса, пожалуйста!

[Lunatik]

Здравствуйте!
Мне нужно разработать программу вычисления матричного синуса с помощью модульного программирования в Delphi.

Имеется формула:

SinX=X-(X^3)/(3!)+(X^5)/(5!)-...+((-1)^n)*(X^(2n+1))/((2n+1)!)

где Х - заданная матрица (nxn)

условие выхода ||(X^m)/(m!)||<eps

eps-заданная точность

||Z||=*сумма от i,j=1 до n*|Zij|

[Stilet]

Умеешь факториал считать?

[Lunatik]

Так?

procedure faktorial(n:integer; var a:real);
var i:integer;
begin
a:=1;
for i:=1 to n do
begin;
a:=a*i;
end;
end;

[Stilet]

Верно, только сделай это функцией:

Function faktorial(n:integer):LongInt;
var i:integer;a:LongInt;
begin
a:=1;
for i:=1 to n do
begin;
 a:=a*i;
end;
Faktorial:=a;
end;

Ну а теперь приблизительно так:

uses math,
...
Function aga(x:real,n:integer):Extended;
begin
 result:=power(-1,n)*power(x,2*n+1)/faktorial(2*n+1);
end;

И в цикле суммируй результы этой функи, передавая в нее X и счетчик цикла.

[Lunatik]

Спасибо, большое!

Но у нас так не прокатит. Надо процедурами написать:
умножение матриц, умножение матриц на число, сложение матриц. А потом все это собрать вместе. Через цикл.

[Stilet]

Ну как хочеш. Пиши процедурами, там переделывать немного.

[Lunatik]

Это не я хочу. Это требуют...