Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

11. FEJEZET. PÉLDAPROGRAMOK 42
function kar_driver
kar driver.m
% adatok
D=0.05; A=D^2*pi/4; L=10; N=10; a=1000;
ro=1000; lambda=0.02; tmax=0.1;
% inicializalas
t=0; dt=L/N/a; x=0:L/N:L; dzdx=zeros(1:N);
% cell{} tip. vektorral több adattipust be lehet csomagolni
adatok{1}=[D L N a ro]; adatok{2}=dzdx;
% Inicializálás (t=0-hoz tartozó állapot)
p=2e5*ones(1,N+1); v=zeros(1,N+1); i=1;
% indulamandula
while t