function plot_series_x3( a, b, nk, np )
%   funzione per confrontare il grafico
%   di 1/(1+x^3) con x in [a,b]
%   con quello della funzione partsum_x3( x, nk )
%   somma parziale nk-esima della serie di arctg(x)
%   np: numero di punti di x 
%   (x = linspace(a,b,np))

figure(nk+1)
clf

x = linspace(a,b,np);
%x2 = linspace(a, 2*b, 2*np);
frac_app = partsum_x3( x, nk );

hp = plot( x, 1./(1+x.^3), 'k', x, frac_app, 'r--');

hold on

plot([a b], [0 0], 'k')

hy = ylabel('\it f(x)', 'rotation', 0, 'FontSize', 12, ...
    'Position', [-1.1 3.5 0]);

set( gca, 'FontSize', 12 );

hx = xlabel( '{\it x}', 'FontSize', 12 );

titlestr1 = ...
    strcat('Grafico della funzione {\it 1/(1+x^3)} ',... 
        ' e della somma parziale della serie.');
titlestr2 = ['\it n_k = ', num2str(nk)];

title( {titlestr1;titlestr2}, 'FontSize', 12 );

hl = ...
    legend(hp, '1/(1+x^3)', '\Sigma_0^{n_k}  a_k(x)',...
    'Location', 'NorthEast');

set( hl, 'FontSize', 12);

axis tight