% Program: sunflower.m
% Author:  Michael Goldwasser
% This is a solution for CSCI-145, asgn03

d = 137.51;
n = 1:1000;

r = sqrt(n);
theta = pi * d .* n / 180;
x = r .* sin(theta);
y = r .* cos(theta);

plot(x,y,'o');
axis off;          % do not display axis markers or white background
axis equal;        % plot both dimensions with equal scale
