Перевод с С или Паскаля на Фортран

@Peperovich · Регистрация: 17.05.2011

Author24 — интернет-сервис помощи студентам

Всем добрый день! Помогите перевести программу на Фортран желательно версии 4. Я написал 1 программу на С++ и Паскале, мне надо еще на Фортране она. Но с Фортраном я вообще не знаком. Помогите кому не сложно.

C++

#include <iostream>
#include <conio.h>
#include <iomanip>
 
using namespace std;
 
int Matrix[3][3];
int TMatrix[3][3];
int MMatrix[3][3];
int A[9];
 
//ввод элементов матрицы
void interMatrix()
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            cout<<"Enter element's of Matrix["<<i<<"]["<<j<<"] = ";
            cin>>Matrix[i][j];
            cout<<endl;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//вывод матрицы в читабельном виде
void OutMatrix(int (&M)[3][3])
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            cout<<setw(4)<<M[i][j];
            j++;
        }
        while(j<3);
        cout<<endl;
        i++;
        j=0;
    }
    while(i<3);
}
 
//вычисляем определитель матрицы
int DeterMatrix()
{
    int i,j,D;
    i=j=D=0;
    D=Matrix[0][0]*Matrix[1][1]*Matrix[2][2]-Matrix[0][2]*Matrix[1][1]*Matrix[2][0];
    D=D+Matrix[0][1]*Matrix[1][2]*Matrix[2][0]+Matrix[0][2]*Matrix[1][0]*Matrix[2][1];
    D=D-Matrix[1][0]*Matrix[0][1]*Matrix[2][2]-Matrix[0][0]*Matrix[2][1]*Matrix[1][2];
    return D;
}
 
//умножение на число
void ChMatrix(int a)
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            Matrix[i][j]=Matrix[i][j]*a;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//получаем транспонированную матрицу
void getTMatrix()
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            TMatrix[i][j]=Matrix[j][i];
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//получаем диагональную матрицу
void getDMatrix()
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            if (i!=j) TMatrix[i][j]=0;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//суммируем матрицы
void sumMatrix()
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            Matrix[i][j]=Matrix[i][j]+TMatrix[i][j];
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//наибольший элемент матрицы
int MaxMatrix()
{
    int i,j,max;
    i=j=0;
    max=-999999;
    do
    {
        do
        {
            if (max<Matrix[i][j]) max=Matrix[i][j];
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
    return max;
}
 
//наименьший элемент матрицы
int MinMatrix()
{
    int i,j,min;
    i=j=0;
    min=999999;
    do
    {
        do
        {
            if (min>Matrix[i][j]) min=Matrix[i][j];
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
    return min;
}
 
//след матрицы
int SpMatrix()
{
    int i,j,Sp;
    i=j=Sp=0;
    do
    {
        do
        {
            if (i==j) Sp=Sp+Matrix[i][j];
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
    return Sp;
}
 
//инициализируем матрицу
void InMatrix()
{
    int i,j,a;
    i=j=a=0;
    do
    {
        do
        {
            Matrix[i][j]=0;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//перемножаем матрицы
void  MultMatrix()
{
    int i,j;
    i=j=0;
    do
    {
        do
        {
            MMatrix[i][j]=Matrix[i][0]*TMatrix[0][j]+Matrix[i][1]*TMatrix[1][j]+Matrix[i][2]*TMatrix[2][j];
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//сортировка по возрастанию
void AscMatrix()
{
    int i,j,c;
    i=j=c=0;
    do
    {
        do
        {
            A[c]=Matrix[i][j];
            c++;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
    i=j=c=0;
    do
    {
        do
        {
            if (A[j]>A[j+1]) 
            {
                c=A[j+1];
                A[j+1]=A[j];
                A[j]=c;
            }
            j++;
        }
        while(j<(8-i));
        i++;
        j=0;
    }
    while(i<8);
    i=j=c=0;
    do
    {
        do
        {
            Matrix[i][j]=A[c];
            c++;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//сортировка по убыванию
void DscMatrix()
{
    int i,j,c;
    i=j=c=0;
    do
    {
        do
        {
            A[c]=Matrix[i][j];
            c++;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
    i=j=c=0;
    do
    {
        do
        {
            if (A[j]<A[j+1]) 
            {
                c=A[j+1];
                A[j+1]=A[j];
                A[j]=c;
            }
            j++;
        }
        while(j<(8-i));
        i++;
        j=0;
    }
    while(i<8);
    i=j=c=0;
    do
    {
        do
        {
            Matrix[i][j]=A[c];
            c++;
            j++;
        }
        while(j<3);
        i++;
        j=0;
    }
    while(i<3);
}
 
//вычисление обратной матрицы
int BackMatrix()
{
    if (DeterMatrix()!=0)
    {
        TMatrix[0][0]=-1*(Matrix[1][1]*Matrix[2][2]-Matrix[2][1]*Matrix[1][2]);
        TMatrix[0][1]=(Matrix[1][0]*Matrix[2][2]-Matrix[1][2]*Matrix[2][0]);    
        TMatrix[0][2]=-1*(Matrix[1][0]*Matrix[2][1]-Matrix[1][1]*Matrix[2][0]);
        TMatrix[1][0]=(Matrix[0][1]*Matrix[2][2]-Matrix[0][2]*Matrix[2][1]);
        TMatrix[1][1]=-1*(Matrix[0][0]*Matrix[2][2]-Matrix[0][2]*Matrix[2][0]); 
        TMatrix[1][2]=(Matrix[0][0]*Matrix[2][1]-Matrix[0][1]*Matrix[2][0]);
        TMatrix[2][0]=-1*(Matrix[0][1]*Matrix[1][2]-Matrix[1][1]*Matrix[0][2]);
        TMatrix[2][1]=(Matrix[0][0]*Matrix[1][2]-Matrix[1][0]*Matrix[0][2]);
        TMatrix[2][2]=-1*(Matrix[0][0]*Matrix[1][1]-Matrix[1][0]*Matrix[0][1]);
        int i,j;
        i=j=0;
        do
        {
            do
            {
                MMatrix[i][j]=TMatrix[j][i];
                j++;
            }
            while(j<3);
            i++;
            j=0;
        }
        while(i<3);
        return 1;
    }
    else
    {
        return 0;
    }
}
 
int main()
{
    int a;
 
    InMatrix();
    interMatrix();
    system("cls");
 
    cout<<setw(10)<<"Matrix"<<endl;
    OutMatrix(Matrix);
    cout<<endl;
 
    cout<<"Max element of Matrix = "<<MaxMatrix();
    cout<<endl;
 
    cout<<"Min element of Matrix = "<<MinMatrix();
    cout<<endl;
 
    cout<<"Determinant Matrix = "<<DeterMatrix()<<endl;
 
    cout<<"Trace Matrix = "<<SpMatrix()<<endl;
    
    cout<<"Enter number = ";
    cin>>a;
    ChMatrix(a);
    cout<<endl;
    cout<<setw(10)<<"Matrix * number"<<endl;
    OutMatrix(Matrix);
 
    getTMatrix();
    cout<<endl;
    cout<<setw(10)<<"TMatrix"<<endl;
    OutMatrix(TMatrix);
 
    getDMatrix();
    cout<<endl;
    cout<<setw(10)<<"Diagonal Matrix from TMatrix"<<endl;
    OutMatrix(TMatrix);
 
    MultMatrix();
    cout<<endl;
    cout<<setw(10)<<"Matrix * TMatrix"<<endl;
    OutMatrix(MMatrix);
 
    sumMatrix();
    cout<<endl;
    cout<<setw(10)<<"Matrix + TMatrix"<<endl;
    OutMatrix(Matrix);
 
    AscMatrix();
    cout<<endl;
    cout<<setw(10)<<"Asceding Matrix"<<endl;
    OutMatrix(Matrix);
 
    DscMatrix();
    cout<<endl;
    cout<<setw(10)<<"Desceding Matrix"<<endl;
    OutMatrix(Matrix);
    
    cout<<endl;
    if (BackMatrix()==0) cout<<"LoopBack Matrix don't exist"<<endl;
    else
    {
        cout<<setw(10)<<"Loopback Matrix"<<endl;
        int i,j,f;
        i=j=f=0;
        do
        {
            if (i!=1)   cout<<"      ";
            do
            {
                if (i!=1)
                    cout<<setw(4)<<MMatrix[i][j];
                else
                {
                    if (f==0)
                    {
                    cout<<"1/"<<DeterMatrix()<<" *"<<setw(4)<<MMatrix[i][j];
                    f=1;
                    }
                    else cout<<setw(4)<<MMatrix[i][j];
                }
                j++;
            }
            while(j<3);
            cout<<endl;
            i++;
            j=0;
        }
        while(i<3);
    }
    getch();
    return 0;
}

Pascal

program Project2;
 
{$APPTYPE CONSOLE}
 
{$R *.res}
 
uses
  System.SysUtils,Windows;
 
Type X=array [1..3,1..3] of integer;
 
var
Matrix:x;
TMatrix:x;
MMatrix:x;
 
A:array [1..9] of integer;
a1,i,j,f:integer;
 
procedure InMatrix;
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    Matrix[i][j]:=0;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure interMatrix;
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    write('Enter elements of Matrix[',i,'][',j,'] = ');
    readln(Matrix[i][j]);
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure OutMatrix(M:x);
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    write('   ',M[i][j]);
    j:=j+1;
  until j>3;
  writeln;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure ChMatrix(a:integer);
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    Matrix[i][j]:=Matrix[i][j]*a;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure getTMatrix;
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    TMatrix[i][j]:=Matrix[j][i];
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
function DeterMatrix:integer;
var
i,j,D:integer;
begin
i:=1;
j:=1;
D:=0;
D:=Matrix[1][1]*Matrix[2][2]*Matrix[3][3]-Matrix[1][3]*Matrix[2][2]*Matrix[3][1];
D:=D+Matrix[1][2]*Matrix[2][3]*Matrix[3][1]+Matrix[1][3]*Matrix[2][1]*Matrix[3][2];
D:=D-Matrix[2][1]*Matrix[1][2]*Matrix[3][3]-Matrix[1][1]*Matrix[3][2]*Matrix[2][3];
DeterMatrix:=D;
end;
 
procedure getDMatrix;
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    if (i<>j) then TMatrix[i][j]:=0;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure sumMatrix;
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    Matrix[i][j]:=Matrix[i][j]+TMatrix[i][j];
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
function MaxMatrix:integer;
var
i,j,max:integer;
begin
i:=1;
j:=1;
max:=-99999;
repeat
  repeat
    if (max<Matrix[i][j]) then max:=Matrix[i][j];
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
MaxMatrix:=max;
end;
 
function MinMatrix:integer;
var
i,j,min:integer;
begin
i:=1;
j:=1;
min:=99999;
repeat
  repeat
    if (min>Matrix[i][j]) then min:=Matrix[i][j];
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
MinMatrix:=min;
end;
 
function SpMatrix:integer;
var
i,j,Sp:integer;
begin
i:=1;
j:=1;
Sp:=0;
repeat
  repeat
    if (i=j) then Sp:=Sp+Matrix[i][j];
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
SpMatrix:=Sp;
end;
 
procedure MultMatrix;
var
i,j:integer;
begin
i:=1;
j:=1;
repeat
  repeat
    MMatrix[i][j]:=Matrix[i][1]*TMatrix[1][j]+Matrix[i][2]*TMatrix[2][j]+Matrix[i][3]*TMatrix[3][j];
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure AscMatrix;
var
i,j,c:integer;
begin
i:=1;
j:=1;
c:=1;
repeat
  repeat
    A[c]:=Matrix[i][j];
    c:=c+1;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
i:=1;
j:=1;
c:=1;
repeat
  repeat
    if (A[j]>A[j+1]) then
            begin
                c:=A[j+1];
                A[j+1]:=A[j];
                A[j]:=c;
            end;
    j:=j+1;
  until (j>(9-i));
  i:=i+1;
  j:=1;
until i>8;
i:=1;
j:=1;
c:=1;
repeat
  repeat
    Matrix[i][j]:=A[c];
    c:=c+1;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
procedure DscMatrix;
var
i,j,c:integer;
begin
i:=1;
j:=1;
c:=1;
repeat
  repeat
    A[c]:=Matrix[i][j];
    c:=c+1;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
i:=1;
j:=1;
c:=1;
repeat
  repeat
    if (A[j]<A[j+1]) then
            begin
                c:=A[j+1];
                A[j+1]:=A[j];
                A[j]:=c;
            end;
    j:=j+1;
  until (j>(9-i));
  i:=i+1;
  j:=1;
until i>8;
i:=1;
j:=1;
c:=1;
repeat
  repeat
    Matrix[i][j]:=A[c];
    c:=c+1;
    j:=j+1;
  until j>3;
  i:=i+1;
  j:=1;
until i>3;
end;
 
function BackMatrix:integer;
var
i,j:integer;
begin
if (DeterMatrix<>0)   then
    begin
        TMatrix[1][1]:=-1*(Matrix[2][2]*Matrix[3][3]-Matrix[3][2]*Matrix[2][3]);
        TMatrix[1][2]:=(Matrix[2][1]*Matrix[3][3]-Matrix[2][3]*Matrix[3][1]);
        TMatrix[1][3]:=-1*(Matrix[2][1]*Matrix[3][2]-Matrix[2][2]*Matrix[3][1]);
        TMatrix[2][1]:=(Matrix[1][2]*Matrix[3][3]-Matrix[1][3]*Matrix[3][2]);
        TMatrix[2][2]:=-1*(Matrix[1][1]*Matrix[3][3]-Matrix[1][3]*Matrix[3][1]);
        TMatrix[2][3]:=(Matrix[1][1]*Matrix[3][2]-Matrix[1][2]*Matrix[3][1]);
        TMatrix[3][1]:=-1*(Matrix[1][2]*Matrix[2][3]-Matrix[2][2]*Matrix[1][3]);
        TMatrix[3][2]:=(Matrix[1][1]*Matrix[2][3]-Matrix[2][1]*Matrix[1][3]);
        TMatrix[3][3]:=-1*(Matrix[1][1]*Matrix[2][2]-Matrix[2][1]*Matrix[1][2]);
    i:=1;
    j:=1;
    repeat
      repeat
        MMatrix[i][j]:=TMatrix[j][i];
        j:=j+1;
      until j>3;
      i:=i+1;
      j:=1;
    until i>3;
    BackMatrix:=1;
   end
   else
   BackMatrix:=0
end;
 
begin
InMatrix;
interMatrix;
WinExec('cls',1);
writeln;
writeln('Matrix');
OutMatrix(Matrix);
writeln;
writeln('Max element of Matrix = ',MaxMatrix);
writeln('Min element of Matrix = ',MinMatrix);
writeln('Determinant Matrix = ',DeterMatrix);
writeln('Trace Matrix = ',SpMatrix);
write('Enter number = ');
readln(a1);
ChMatrix(a1);
writeln;
writeln('Matrix * number');
OutMatrix(Matrix);
getTMatrix;
writeln;
writeln('TMatrix');
OutMatrix(TMatrix);
getDMatrix;
writeln;
writeln('Diagonal Matrix from TMatrix');
OutMatrix(TMatrix);
MultMatrix;
writeln;
writeln('Matrix * TMatrix');
OutMatrix(MMatrix);
sumMatrix;
writeln;
writeln('Matrix + TMatrix');
OutMatrix(Matrix);
AscMatrix;
writeln;
writeln('Asceding Matrix');
OutMatrix(Matrix);
DscMatrix;
writeln;
writeln('Desceding Matrix');
OutMatrix(Matrix);
writeln;
if (BackMatrix=0) then writeln('LoopBack Matrix dont exist')
    else
    begin
    writeln('Loopback matrix');
   i:=1;
   j:=1;
   f:=0;
    repeat
      if (i<>2) then write('    ');
      repeat
        if (i<>2) then write(' ',MMatrix[i][j])
        else
        begin
          if (f=0) then
          begin
            write('1/',DeterMatrix,'* ',MMatrix[i][j],' ');
            f:=1;
          end
          else write(' ',MMatrix[i][j]);
        end;
        j:=j+1;
      until j>3;
      writeln;
      i:=i+1;
      j:=1;
    until i>3;
 
  end;
 
 
 
 
 
 
readln;
end.

@Kastaneda · 14.01.2013, 18:53

Там далеко не 5 строк кода, поэтому вряд ли кто-то возьмется переписывать. Могу посоветовать поискать конвертеры C to Fortran. Точно знаю, что трансляторов с Фортрана на С в сети полно.

@Peperovich · 15.01.2013, 09:06 **[ТС]**

Да там не 5 строк, но там процедуры все однотипны.

@Kastaneda 5231 / 3204 / 362 Регистрация: 12.12.2009 Сообщений: 8,116 Записей в блоге: 2
	14.01.2013, 18:53	2
	Там далеко не 5 строк кода, поэтому вряд ли кто-то возьмется переписывать. Могу посоветовать поискать конвертеры C to Fortran. Точно знаю, что трансляторов с Фортрана на С в сети полно. 0

@Peperovich 43 / 43 / 13 Регистрация: 17.05.2011 Сообщений: 162
	15.01.2013, 09:06 [ТС]	3
	Да там не 5 строк, но там процедуры все однотипны. 0