Табулирование функции (полярные коор.)

[Su-34]

Помогите пжлста. Мне надо затабулировать функцию р=sin(f/2), заданную полярными координатами с помощью Visual Fortran.
Что я не так сделал, что нужно исправить чтобы в екселе получился график 2). Пока он мне рисует только правую часть 1)

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