Нарисовать шляпку подсолнуха с логарифмическими спиралями - Turbo Pascal

Формулировка задачи:

Программа должна рисовать шляпку подсолнуха с логарифмическими спиралями. Этот код работает, но нужно чтобы получалось что-то типа того что на картинке. Помогите, пожалуйста, кто сможет.

uses graph;
var x0,y0,x,y,i:integer;
    a,b,t,f:real;
begin
initgraph(x,y,'');
x0:=getmaxX div 2;
y0:=getmaxY div 2;
setcolor(2);
circle(x0,y0,153);
a:=1;
b:=0.5;
f:=0;
for i:=1 to 12 do
 begin
  t:=0;
  repeat
   t:=t+0.001;
   x:=round(a*exp(b*(t))*cos(t+f));
   y:=round(a*exp(b*(t))*sin(t+f));
   putpixel(x0+x,y0-y,14);
  until t>3.2*pi;
  f:=f+pi/6;
 end;
readln
end.

Решение задачи: «Нарисовать шляпку подсолнуха с логарифмическими спиралями»

Листинг программы

uses graph;
var x0,y0,x,y,i:integer;
    a,b,t,f:real;
begin
initgraph(x,y,'');
x0:=getmaxX div 2;
y0:=getmaxY div 2;
setbkcolor(15);
a:=1;
b:=0.5;
f:=0;
for i:=1 to 21 do
 begin
  t:=0;
  repeat
   t:=t+0.01;
   x:=round(a*exp(b*(t))*cos(t+f));
   y:=round(a*exp(b*(t))*sin(t+f));
   putpixel(x0+x,y0-y,8);
  until t>3.2*pi;
  f:=f+2*pi/21;
 end;
for i:=1 to 34 do
 begin
  t:=0;
  repeat
   t:=t+0.01;
   x:=round(a*exp(b*(t))*sin(t+f));
   y:=round(a*exp(b*(t))*cos(t+f));
   putpixel(x0+x,y0-y,8);
  until t>3.2*pi;
  f:=f+2*pi/34;
 end;
setcolor(2);
circle(x0,y0,155);
readln
end.

Объяснение кода листинга программы

Подключается библиотека graph.
Объявляются переменные x, y, i, a, b, t, f: integer; a, b, t, f: real.
Вызывается функция initgraph(x, y, ''), которая устанавливает графический контекст.
Вычисляются координаты центра шляпки подсолнуха x0, y0, разделив максимальные координаты по горизонтали и вертикали на 2.
Устанавливается цвет фона с помощью функции setbkcolor(15).
Инициализируются переменные a, b, t, f.
Задается начальное значение для переменной t в 0.
Задается начальное значение для переменной f в 0.
Запускается цикл for, который будет выполняться 21 раз.
В каждой итерации цикла вычисляются значения x и y по формуле x = a exp(b (t)) cos(t + f) и y = a exp(b (t)) sin(t + f).
Выводится пиксель на экран с помощью функции putpixel(x0 + x, y0 - y, 8), где x и y — вычисленные значения x и y.
Значение переменной t увеличивается на 0.01.
Значение переменной f увеличивается на 2 * pi / 21.
После завершения каждой итерации цикла выводится сообщение readln.
Заканчивается функция.
Вызывается функция circle(x0, y0, 155), которая рисует окружность с центром в точке (x0, y0) радиусом 155.
Завершается программа.