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