این برنامه موج یک بعدی مرتبه دوم می باشد که فقط شکل موج را به صورت انیمیشن ترسیم می کند. در این برنامه جابجایی اولیه و سرعت اولیه رو به صورت تابع بیان کردم. برنامه فوق گزارش کار ندارد. n = 101; % number of space gridpoints without boundaries dt = 1.6e-2; % time step tf = 8e0; % final time x = linspace(-1,1,n+2)'; x = x(2:end-1); h = x(2)-x(1); r = dt/h; disp(sprintf( 'CFL number: %0.2f' ,r)) u0 = f(x); U = [u0,u0+dt*g(x)];... ...