Operaciones con matrices

Suma resta y multiplicación de matrices en pascal.

Para el uso de estos procedimientos es necesario declarar el tipo
tMatriz. Este tipo es un registro en el que f y c son las filas y columnas
de la matriz y datos es una matriz de numeros reales.

const
 N = 10;

type
 tMatriz= record
 datos: array [1..N,1..N] of real;
        f,c: 0..N;
    end;



Tarea: Suma de dos matrices
El procedimiento recibe dos matrices y devuelve la matriz sumada.


procedure Sumar(m1, m2: tMatriz; var resul: tMatriz);
var i,j:integer;
begin
      if (m1.c = m2.c) and (m1.f = m2.f) then begin
          resul.c:=m1.c;    { Crear Matriz resul del mismo rango que m1 y m2 }
          resul.f:=m2.c;
          for i:=1 to m1.c do
    for j:=1 to m1.f do
                        begin
                         resul.datos[i,j]:= m1.datos[i,j] + m2.datos[i,j];
                        end;
          end
       else writeln('No se pueden sumar matrices de diferente dimension');
end;



Tarea: Resta de dos matrices
El procedimiento recibe dos matrices y devuelve la matriz restada.


procedure Restar(m1, m2: tMatriz; var resul: tMatriz);
var i,j:integer;
begin
if (m1.c = m2.c) and (m1.f = m2.f) then begin
resul.c:=m1.c;    { Crear Matriz resul del mismo rango que m1 y m2 }
          resul.f:=m2.c;
          for i:=1 to m1.c do
                    for j:=1 to m1.f do
                     begin
                      resul.datos[i,j]:= m1.datos[i,j] - m2.datos[i,j];
                     end;
          end
       else writeln('No se pueden restar matrices de diferente dimension');
end;


Tarea: Multiplicación dos matrices
El procedimiento recibe dos matrices y devuelve la matriz multiplicada.


procedure Multiplicar(m1, m2: tMatriz; var resul: tMatriz);
var i,j:integer;
begin
        if (m1.c = m2.f) then begin
          resul.c:=m1.f;    { Crear Matriz resul de rango apropiado para multiplicar }
          resul.f:=m2.c;
          for i:=1 to m1.c do
                    for j:=1 to m1.f do
                     begin
                      resul.datos[i,j]:=resul.datos[i,j] + m1.datos[i,j] *
m2.datos[i,columna];
                     end;
end
       else writeln('No se pueden multiplicar matrices de esta dimension');
end;