Mercurial > repos > public > sbplib

comparison +scheme/Euler1d.m @ 78:80948a4084f3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added the factor or 1/2 to the wall bc term as well.
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 25 Nov 2015 15:41:10 +0100
parents 0b07ff8a0a12
children 2102af217134
comparison
equal deleted inserted replaced
77:0b07ff8a0a12 78:80948a4084f3
410 % w_s = T\q_s; 410 % w_s = T\q_s;
411 % w_s = Tinv * q_s; % Med analytisk matris 411 % w_s = Tinv * q_s; % Med analytisk matris
412 w_in = w_s(p_in); 412 w_in = w_s(p_in);
413 w_ot = w_s(p_ot); 413 w_ot = w_s(p_ot);
414 414
415 o = obj.Hi * tau * (w_in - R*w_ot); 415 o = 1/2*obj.Hi * tau * (w_in - R*w_ot);
416 end 416 end
417 417
418 closure = @closure_fun; 418 closure = @closure_fun;
419 end 419 end
420 420