function plot_series_cos2( a, b, nk, np )
%   funzione per confrontare il grafico
%   di cos(x^2) con x in [a,b]
%   con quello della funzione partsum_cos2( x, nk )
%   somma parziale nk-esima della serie di cos(x^2)
%   np: numero di punti di x 
%   (x = linspace(a,b,np))

x = linspace(a,b,np);
x2 = linspace(a, b, np);
cos2_app = partsum_cos2( x, nk );

hp = plot( x2, cos(x2.^2), 'k', x, cos2_app, 'r--');

hy = ylabel('\itf(x)', 'rotation', 0, 'FontSize', 12);

set( gca, 'FontSize', 12 );

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

titlestr1 = ...
    'Grafico della funzione {\itcos(x^2)} e della somma parziale della serie.';
titlestr2 = ['\it n_k = ', num2str(nk)];

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

hl = ...
    legend(hp, 'cos(x^2)', '\Sigma_0^{n_k}  a_k(x)',...
    'Location', 'SouthWest');

set( hl, 'FontSize', 12);

axis tight