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

[ScoRpiON]

Помогите пжлста. Мне надо затабулировать функцию р=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

[Selyd]

Я понимаю так.
Угол надо менять от 0 до 4pi, чтобы обойти один период.
Беру a=0, b=4pi=12.57, N=40(50), сколько хочешь.
К программе, набираешь таблицу данных и в Exel.