Bird plot by mathematical equation

Sobhan Bamdad writes this code in Policlinico and It plots a free Bird using some mathematical equations from Hamid Naderi Yegane.

clc;clear all;close all;

%%% This code is written by Sobhan Bamdad in Policlinico

% to plot a free Bird by some mathematical equations from Hamid Naderi Yegane

fc_circle_plot = @(xs,ys,rs,varargin)...

arrayfun(@(x,y,r)plot(x,y,'o','MarkerSize',r,varargin{:}),xs,ys,rs,'UniformOutput',false);

k=-20000:1:20000;

x=(3*k/45000)+ sin((17*pi/20)*(k/20000).^5) .* (cos(41*pi*k/20000)).^6 +((1/3)* ((cos(41*pi*k/20000)).^16 ) +...

(1/3)* (cos(41*pi*k/20000)).^80) .* ((cos(pi*k/40000)).^12) .*(sin(6*pi*k/20000));

y=((15/30)*(k/20000).^4) - cos((17*pi/20)*(k/20000).^5).* ( (11/10)+(45/20)*((cos(pi*k/40000)).^8) ...

.*((cos(3*pi*k/40000)).^6) ) .* ((cos(41*pi*k/20000)).^6) + (12/20)*((cos(3*pi*k/200000)).^10).* ...

((cos(9*pi*k/200000)).^10) .* ((cos(18*pi*k/200000)).^10);

R=(1/2)+ (1/2)*((sin(41*pi*k/20000)).^2).* ((sin(9*pi*k/200000)).^2) + ...

(1/2)*((cos(41*pi*k/20000)).^2).*((cos(pi*k/40000)).^10) ;

figure;hold on;

fc_circle_plot(x,y,R,'Color','m','LineWidth',2);hold off;