function btv_test06 ( ) %*****************************************************************************80 % %% BTV_TEST06 tests BURGERS_TIME_VISCOUS with the spline initial condition. % % Licensing: % % This code is distributed under the GNU LGPL license. % % Modified: % % 23 April 2012 % % Author: % % John Burkardt % fprintf ( 1, '\n' ); fprintf ( 1, 'BTV_TEST06\n' ); fprintf ( 1, ' Test BURGERS_TIME_VISCOUS with the spline initial condition.\n' ); fprintf ( 1, ' Use periodic boundaries.\n' ); nx = 81; nt = 200; t_max = 2.0; nu = 0.01; bc = 4; fprintf ( 1, '\n' ); fprintf ( 1, ' Initial condition: spline\n' ); fprintf ( 1, ' Number of space nodes = %d\n', nx ); fprintf ( 1, ' Number of time steps = %d\n', nt ); fprintf ( 1, ' Final time T_MAX = %g\n', t_max ); fprintf ( 1, ' Viscosity = %g\n', nu ); fprintf ( 1, ' Boundary condition = %d\n', bc ); U = burgers_time_viscous ( @ic_spline, nx, nt, t_max, nu, bc ); x = linspace ( -1.0, +1.0, nx ); figure ( 6 ) plot ( x, U(1:50:(nt+1),:), 'Linewidth', 3 ) grid on xlabel ( '<-- X -->' ) ylabel ( '<-- U(X,T) -->' ) title ( 'Burgers equation solutions over time, initial condition spline' ) filename = 'btv_test06.png'; print ( '-dpng', filename ) fprintf ( 1, '\n' ); fprintf ( 1, ' Saved plot as "%s"\n', filename ); return end