برنامه متلب حل معادله موج یک بعدی مرتبه دو به روش صریح

این برنامه موج یک بعدی مرتبه دوم می باشد که فقط شکل موج را به صورت انیمیشن ترسیم می کند. در این برنامه جابجایی اولیه و سرعت اولیه رو به صورت تابع بیان کردم. برنامه فوق گزارش کار ندارد. 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)];...     ...