Index: mttroot/mtt/bin/trans/dae2daes_m ================================================================== --- mttroot/mtt/bin/trans/dae2daes_m +++ mttroot/mtt/bin/trans/dae2daes_m @@ -13,19 +13,22 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.1 1996/08/15 16:24:55 peter +## Initial revision +## ############################################################### echo Creating $1_daesol.m rm -f dae2daesol_m.log if [ "$2" = "" ]; then - PARAMS='T=[0:0.1:10];x0=zeros(10,1),dx0=zeros(10,1)' + PARAMS='T=[0:0.1:10];x0=zeros(nx+nz,1);dx0=zeros(nx+nz,1);' echo Using default parameter $PARAMS else PARAMS=$2; fi @@ -59,11 +62,11 @@ x = dassl('$1_dae', x0, dx0, T); i=0; for tt=T' i=i+1; - y(i) = $1_daeo(x(i,:),tt); + y(i,:) = $1_daeo(x(i,:),tt)'; end; write_matrix([T,y], '$1_daesol'); EOF