Подпрограммы общего вида. (Fortran).

[Moustache]

Добрый вечер. Написал программу в фортране, но мне выдаются ошибки.

Само задание звучит так:
Составить подпрограмму нахождения индекса минимального элемента кратного заданному целому числу в целочисленном одномерном массиве. С помощью этой подпрограммы сформировать вектор из элементов целочисленной матрицы A(-3:0,-4:0)

Вот мое решение:

program LAB2
dimension A(-3:0,-4:0), B(5)
real min
real c
do i=-3,0
read *, (a(i,j),j=-4,0)
enddo
do i=-3,0
do j=-4,0
b(i)=a(i,j)
enddo
call pr(a,-3,0,-4,0,b,5,c)
min=b(i)
enddo
end
subroutine pr(a,i1,i2,j1,j2,b,k,c)
dimension A(i1:i2,j1:j2), b(k)
real min
real c
read *,c
print * , 'vvedite chislo'
min=b(1)
do i=1,n
if (b(i).LT.min) then
if ((b(i)/c).EQ.(b(i)/c)) then
min=b(i)
l=i
endif
endif
enddo
end subroutine

[Serebro]

Сообщение от Moustache

if ((b(i)/c).EQ.(b(i)/c)) then

1. Что это такое? (b(i)/c) всегда равно (b(i)/c).

2. Для чего в подпрограмме pr(...) требуется массив a(...)?

[Moustache]

Сообщение от Serebro

1. Что это такое? (b(i)/c) всегда равно (b(i)/c).

2. Для чего в подпрограмме pr(...) требуется массив a(...)?

1. Пытаюсь задать кратность. Если я пишу ((b(i)/c).EQ.(b(i)/c.)) , то программа вообще не запускается.

2. Хм. Ну если я не ошибаюсь, то нужно перечислить все элементы кода, которые в нем используются.

[Serebro]

Сообщение от Serebro

1. Что это такое? (b(i)/c) всегда равно (b(i)/c).

2. Для чего в подпрограмме pr(...) требуется массив a(...)?

1. Кратность проверяется так:

integer n,m,mmm

n=mmm - (mmm / m)*m                  // или n=mmm - int((mmm / m))*m 

if (n.eq.0) then
              //      mmm - кратно m
end if

2. В подпрограмме pr(...) массив a(...) НЕ ИСПОЛЬЗУЕТСЯ, поэтому и ставится вопрос:

2. Для чего в подпрограмме pr(...) требуется массив a(...)?

[Moustache]

Без понятия.