программа на фортране с массивами

[ЮлияБаля]

Я понимаю что это нечто но помогите пожалуйста я вообще не соображаю(
своими словами
есть уравнение его решение я предусмотрела но суть в том что теперь нужно подготовить программу, чтобы функцию С(х,р) при изменении значения аргумента X от значения xl до значения х2 с шагом Нх и аргумента Р - от p1 до р2 с шагом Нр. Результаты напечатайте в виде таблицы-матрицы: т.е. мы задает начальные знач х и р, шаг и кол-во шагов а потом когда все расчитается это должно быть выведено ввиде табл. Если возник ошибка (ноль в знаменателе) то в табл в этой графе будет набор звездочек.

При вводе данных организуйте контроль на достоверность: при положительном шаге Нx недопустимо, чтобы xl было боль¬ше х2, а при отрицательной величине ша-га Нx нельзя задавать xl < х2 (аналогично для Нp );

Вот я накидала немного но с выводом я не знаю что делать((

Program massiv
	N=5
	M=5
	Real p(N), x(N), c(N,M) 
	Character S(N,M)
	write(*,*) 'vvedite x1'
	read(*,*) x1
	write(*,*) 'vvedite x2'
	read(*,*) x2
	write(*,*) 'vvedite Hx'
	read(*,*) Hx
	write(*,*) 'vvedite p1'
	read(*,*) p1
	write(*,*) 'vvedite p2'
	read(*,*) p2
	write(*,*) 'vvedite Hp'
	read(*,*) Hp

	Nx=abs((x2-x1)/Hx)+1;
 	Np=abs((p2-p1)/Hp)+1;

	x(1)=x1
	p(1)=p1
 	do 1,i=2,Nx
1	x(i)=x1+(i-1)*Hx
	do 2,i=2,Np
2	p(i)=p1+(i-1)*Hp
	do 3,i=1,N
	do 3,j=1,M

3	buf1=x(i)*x(i)*x(i)-1.5
	if (buf1) 10,20,30
	

10	write(*,*) S(i,j)='*'
	goto 1000


20	write(*,*) S(i,j)='***'
	goto 1000
	
30	buf2=abs(x(i)*x(i)-(p(j)/(x(i)*x(i)-1.5)))-0.2
	
	if (buf2) 40,40,50

40	write(*,*) S(i,j)='**'
	goto 1000

50	buf3=log((abs(x(i)*x(i)-(p(j)/(x(i)*x(i)-1.5)))-0.2))

	if (buf3) 10,60,60
	
	
60	buf4=sqrt(log((abs(x(i)*x(i)-(p(j)/(x(i)*x(i)-1.5)))-0.2)))

	if (buf4) 70,80,70

80	write(*,*) S(i,j)='***'
	goto 1000
	
70	C(i,j)= (p(j)*x(i)-sqrt(buf1)-sin(x(i)))/(buf4)
	write(*,*) C
	
1000	end ?
	

	write(*,*)  'x/p    '
	do 4,j=2,M
4	write(*,*) p(j)

write(*,*) '*-net kornya iz otr, **-logarifma net, ***-0 v znamenatele'

[s-andriano]

Интересно, какой Фортран позволяет так описывать константы?

N=5
M=5

А вот это:

Nx:=abs((x2-x1)/Hx)+1;
Np:=abs((p2-p1)/Hp)+1;

больше похоже на Паскаль, чем на Фортран.

И еще пару советов:
1. Оформляйте код, как это принято.
2. Не забывайте про комментарии.
3. Напишите точно, что Вам нужно.

[ЮлияБаля]

Я не знаю как правильно( я привыкла к паскалю и мне сложно переключится
А задание.. У меня было первое задание написать решение уравнения относит х и р.
А теперь нужно это задание первое преобразовать. т.е. теперь нужно решить уравнение на наборе значений х и р. Но еще есть шаг. т.е. например я ввожу х1=0.9 х2=3 а р1=0.7 р2=5 и по этим формулам рассчитывается шаг Nx:=abs((x2-x1)/Hx)+1;
Np:=abs((p2-p1)/Hp)+1;
а главная задача вывести результаты в виде таблицы когда х(1)=0.9 а p(1)=0.7, потом х(2)=0.9+шаг и р(2)=0,7+шаг и т.д.

[s-andriano]

Давайте по порядку.
Какой стандарт Фортрана Вы используете?
Какой компилятор Вы используете?
Компилируется ли Ваша программа? Если нет - какова диагностика ошибок?

И по поводу первых двух строчек.
Я не помню навскидку, как точно записывается определение констант в Фортране, но, как мне кажется, в зависимости от диалекта, должно быть что-то вроде

integer (parameter N=5)

или

integer, parameter :: N=5

[ЮлияБаля]

Я пишу на 77 Фортране
Компилятор 95
Сегодня была у преподавателя и в таком виде программа компилируется
Program massiv

Real p(5), x(5), c(5,5)
Character S(5,5)
N=5
M=5
write(*,*) 'vvedite x1'
read(*,*) x1
write(*,*) 'vvedite x2'
read(*,*) x2
write(*,*) 'vvedite Hx'
read(*,*) Hx
write(*,*) 'vvedite p1'
read(*,*) p1
write(*,*) 'vvedite p2'
read(*,*) p2
write(*,*) 'vvedite Hp'
read(*,*) Hp

Nx=abs((x2-x1)/Hx)+1;
Np=abs((p2-p1)/Hp)+1;

x(1)=x1
p(1)=p1
do 1,i=2,Nx
1 x(i)=x1+(i-1)*Hx
do 2,i=2,Np
2 p(i)=p1+(i-1)*Hp
do 3,i=1,N
do 3,j=1,M

3 buf1=x(i)*x(i)*x(i)-1.5
if (buf1) 10,20,30


10 write(*,*) S(i,j)='*'
goto 1000


20 write(*,*) S(i,j)='***'
goto 1000

30 buf2=abs(x(i)*x(i)-(p(j)/(x(i)*x(i)-1.5)))-0.2

if (buf2) 40,40,50

40 write(*,*) S(i,j)='**'
goto 1000

50 buf3=log((abs(x(i)*x(i)-(p(j)/(x(i)*x(i)-1.5)))-0.2))

if (buf3) 10,60,60


60 buf4=sqrt(log((abs(x(i)*x(i)-(p(j)/(x(i)*x(i)-1.5)))-0.2)))

if (buf4) 70,80,70

80 write(*,*) S(i,j)='***'
goto 1000

70 C(i,j)= (p(j)*x(i)-sqrt(buf1)-sin(x(i)))/(buf4)
write(*,*) C

1000 end ?


write(*,*) 'x/p '
do 4,j=2,M
4 write(*,*) p(j)







'*-net kornya iz otr, **-logarifma net, ***-0 v znamenatele'

Последняя строчка записана через format правда не помню точно как(

[s-andriano]

Повторюсь:
1. Оформляйте код, как это принято.
2. Не забывайте про комментарии.

1 - очень трудно читать.
2 - до меня не доходит логика, которой Вы руководствовались, когда писали программу. Поэтому комментарии нужны чтобы:
- понять, что же Вы хотите сделать,
- на Фортране пишете именно то, что хотели.

И еще: все это компилироваться точно не может, но, возможно, Ваш компилятор (кстати, Вы так и не написали о нем) что-то пропускает.
Но последний оператор в программе - только end и никогда не write.

[ЮлияБаля]

Ладно ясно.. буду пытаться сама как нибудь! А то мы так все равно не разберемся(