xs = [ 0.0, 0.3, 0.9, 1.7, 2.2, 2.7, 3.7 ];
ys = [-3.0,-3.3,-4.1,-5.5,-6.6,-8.1,-10.9];
xm  =[ 0.0, 0.3, 0.6, 0.8, 1.2, 1.5, 1.75, 2.1, 2.4, 2.7, 3.0, 3.3 ];
ym1 =[-2.1,-2.5,-2.6,-3.1,-3.1,-3.4,-4.1, -3.3,-3.7 -3.4,-4.2,-4.4 ];
ym2 =[-2.5,-2.5,-2.5,-2.2,-3.0,-3.4,-3.8, -3.9,-3.6,-3.9,-4.0,-3.8 ];
ym3 =[-2.4,-2.4,-2.7,-3.2,-3.3,-3.9,-3.5, -3.4,-3.6,-4.1,-3.9,-4.3 ];
ym4 =[-2.3,-2.3,-2.4,-2.8,-3.3,-3.2,-3.2, -4.2,-3.1,-3.3,-4.4,-4.1 ];
figure ( 1 )
clf
xmc = [ 0.0,  4.0 ];
ymc = [ -2.0, -4.0 ];

figure ( 1 )
hold on
axis( [0,4,-11,-2] )
plot ( xm, ym1, 'c.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym2, 'g.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym3, 'b.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym4, 'm.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xmc, ymc, 'k-', 'LineWidth', 2 )
grid on
xlabel ( 'Log10(N), number of points' )
ylabel ( 'Log10(E), integral error' )
title ( 'Monte Carlo Estimates, Stochastic Diffusion Problem, Dimension 11' )

figure ( 2 )
hold on
axis( [0,4,-11,-2] )
plot ( xs, ys,  'r.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym1, 'c.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym2, 'g.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym3, 'b.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xm, ym4, 'm.-', 'LineWidth', 3, 'MarkerSize', 30 );
plot ( xmc, ymc, 'k-', 'LineWidth', 2 )
grid on
xlabel ( 'Log10(N), number of points' )
ylabel ( 'Log10(E), integral error' )
title ( 'Sparse Grid and Monte Carlo Estimates, Stochastic Diffusion Problem, Dimension 11' )
