Программа по вычислению объёма функции z=x^2+y^2

[ivan.tiran]

Прошу прощения, что задаю вопрос в разделе "Паскаль". Мне нужен раздел "Fortran", а его здесь нет. Но вы прогеры умные, думаю, мне поможете.
Вот код:

		program prog
		real*8 s(0:10000,0:10000)
		real*8 x1(10000),x(10000),y1(10000),y(10000)
		real*8 dlina,shag,v
		integer*4 n,i,j
		x1(0)=0
		y1(0)=x1(0)
		v=0
		print*,'Vvedite dlinu uchastka'
		read(*,*)dlina
		print*,'Vvedite shag'
		read(*,*)shag
		n=dlina/shag
		
		
		do i=1,n
		x1(i)=x1(i-1)+shag/2
		x(i)=x1(i)
		print*,'x[',i,']=',x(i)
		x1(i)=x1(i)+shag/2
		enddo
		
		do j=1,n
		y1(j)=y1(j-1)+shag/2
		y(j)=y1(j)
		print*,'y[',j,']=',y(j)
		y1(j)=y1(j)+shag/2
		enddo
		
		
		do i=1,n
		do j=1,n
		s(i,j)=(y(i)**2)+(x(j)**2)*(shag**2)
		v=v+s(i,j)
		enddo 
		enddo
		print*,v
      end program prog

Может я ослеп и не вижу простейшей ошибки. Но в чём причина не могу понять. Когда заполняю массив с иксами, то шаг начинается с 0,5, как мне и нужно. А игрек не с 0,5, а 1,5. Если меняю местами циклы, то такая же ерунда происходит с точностью наоборот. Пробовал в один цикл записывать - там такая же ерунда.



________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)
Не забывайте об этом!
Модератор.

[ivan.tiran]

Да, для примера беру длину отрезка за 10, а шаг за 1.

[NetSpace]

В вузовских старых учебниках по МатАнализу (автор Демидович) есть в конце примеры программ на фортране. Думаю, можете там поискать.

[Anatole]

Я с фортраном работал где-то в конце 80-х. И как помнится компилятор не контролирует выход индексов за пределы массива и индексация массива начинается с 1. Проверте не выходят ли у вас индексы за пределы массива при i или j=1

[TinMan]

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

компилятор не контролирует выход индексов за пределы массива и индексация массива начинается с 1. Проверте не выходят ли у вас индексы за пределы массива при i или j=1

Именно так.
У автора темы без тени сомнения использует x1(0) и y1(0) - а их не существует. В цикле же следующий элемент считается через предыдущий, и цикл начинается с первого. Нужно сменить описания массивов с
x1(10000)
на
x1(0:10000)

[offtop]С Фортраном вообще хохма.. Как-то раз я испортил себе единицу, передавая ее в параметрах, и дня два искал причину совершенно непонятного поведения программы [/offtop]

[ivan.tiran]

Спасибо, всё работает.