برنامه فوق برای حل معادله موج یک بعدی درجه یک به روش اختلاف بالادست می باشد. این برنامه شکل موج را به صورت انیمیشن ترسیم می کند و دارای شرط مرزی دریکله در ابتدا و انتها می باشد. شرط اولیه را تابع sin(x) در نظر گرفتم. این برنامه فاقد گزارش کار می باشد.
clear all
close all
clc
dx=(pi/12);% interval
dt=0.1; % time interval
c=1; %velocity of wave
v=(c*dt/dx);% noo number
.
.
.
for i=1:Jx
u1(i)=sin(x(i));
end
.
.
.
برنامه متلب حل معادله موج یک بعدی درجه یک به روش اختلاف بالادست