Табулирование функций
 

Есть функция abs(cos(3*x)). В Фортране мне её надо затабулировать. Проверьте пжлста коды. В екселе у меня получился график , но приподнятый по у вверх.

program lab1
real a,b,h,x,y
integer i,N
open(10,file='lab1.txt')
a=0
b=(2./3.)*(3.141593)
N=31
h=(b-a)/(N-1)
do i=1,N
x=a+(i-1)*h
y=abs(cos(3*x))
print*,x,y
write(10,*)x,y
end do
end

Слушайте, сударь - ну что Вы вцепились, как бульдог, в это N=31? Что за магическое число такое? У Вас явно просто шаг грубоват - ну так увеличьте своё N хоть до 200, и посмотрите, что будет.
Кстати, что за глупость, причем повторяемая во всех программах:
N=31
h=(b-a)/(N-1) ???
Почему бы не
N=30
h=(b-a)/N ?