Заштрихованный круг

clear, clc, close all
 
fprintf('\n\t Введите радиус круга:')
R = input(' ');
 
 
fprintf('\t Введите число штриховых линий:')
n0 = input(' ');
n = n0;
 
fprintf('\t Введите угол наклона штриховки круга в градусах угловых:')
alpha0 = input(' ');
alpha2 = alpha0*pi/180; %% перевод в радианы
 
alpha = (alpha2 - pi/2);
x = linspace(-R, R, 200); %% область определения, 200 точек
y = sqrt(R^2 - x.^2);    %% относительно центра координат
x = [x, x(end:-1:1)]; 
y = [y, -y]; 
 
n = n + 2;
 
Ax = linspace(-R, R, n);
%% Формирование двумерного массива из блоков
Lx = repmat(Ax, 2, 1); %% см. help repmat
 
Ly = sqrt(R^2 - Lx.^2);
Ly(1, :) = -Ly(1, :);
 
%%% Поворот осей координат, поворот штриховых линий
NLx = Lx*cos(alpha) - Ly*sin(alpha);
NLy = Lx*sin(alpha) + Ly*cos(alpha);
 
%%% Построение и оформление диаграммы
line(x, y, 'linew',2); %% окружность
line(NLx, NLy, 'color', 'r'); %% штриховые линии
grid on;
 
xlabel('\bf\fontsize{12} Ось абсцисс');
ylabel('\bf\fontsize{12} Ось ординат');
 
str = sprintf('%sРадиус R = %g; число штрихов N = %d; угол наклона штрихов %s = %g^o','\bf',...
R, n0, '\fontsize{14}\alpha\fontsize{10}',alpha0);
title(str);
 
set(gca,'fontsize',8, 'fontweight', 'bold');
set(gcf,'name', 'Штриховка круга', 'color','w'); %% белый цвет
axis equal;
ylim([-R - R/100,R + R/100])