ADDED mttroot/mtt/bin/trans/m/mtt_implicit.m Index: mttroot/mtt/bin/trans/m/mtt_implicit.m ================================================================== --- /dev/null +++ mttroot/mtt/bin/trans/m/mtt_implicit.m @@ -0,0 +1,15 @@ +function [x] = mtt_implicit(x,dx,AA,AAx,dt,Nx,open); # Implicit update + +## ############################################################### +## ## Version control history +## ############################################################### +## ## $Id$ +## ## $Log$ +## ############################################################### + +## Copyright (C) 1999 by P.J. Gawthrop + + I = nozeros(!open.*[1:Nx]'); # Indices of states to update + x(I) = AA(I,I)\(AAx(I) + dx(I)*dt); # Implicit update (exept open switches); + +endfunction