t = 0:0.001*pi:pi; n = 9; A = zeros(n,1001); rad = t; A(1,:) = (sin(t)-t.*cos(t))/pi./(1-cos(t)); A(2,:) = (t-sin(t).*cos(t))/pi./(1-cos(t)); fori = 3:n A(i,:) = 2*(sin(i*t).*cos(t)-i*cos(i*t).* sin(t))/i/pi./(i*i-1)./(1-cos(t)); end rad = t/pi*180; fori = 1:n plot(rad,A(i,:)); hold on; end grid on; xlabel('θ/°'); ylabel('αn');