Su-34

Можно еще вопрос? То же самое надо сделать с р=(phi)^2. Ведь должна получиться парабола? Тогда надо писать:
program lab1
real a,b,h,x,y
integer i,N
open(10,file='lab1.txt')
a=-1
b=1
N=31
h=(b-a)/(N-1)
do i=1,N
f=a+h*i
r=f**2
x=r*cos(f)
y=r*sin(f)
print*,x,y
write(10,*)x,y
end do
end

Что не так?