source: trunk/sample/miscellaneous/sample_arrays2.mso @ 289

#*-------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
*
* This LIBRARY is free software; you can distribute it and/or modify
* it under the therms of the ALSOC FREE LICENSE as available at
* http://www.enq.ufrgs.br/alsoc.
*
* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
* from http://www.rps.eng.br Copyright (C) 2002-2004.
* All rights reserved.
*
* EMSO is distributed under the therms of the ALSOC LICENSE as
* available at http://www.enq.ufrgs.br/alsoc.
*
*----------------------------------------------------------------------
* Author: Rafael de Pelegrini Soares
* $Id: sample_arrays2.mso 176 2007-03-04 04:56:54Z arge $
*--------------------------------------------------------------------*#


Model Arrays
        PARAMETERS
        Dimension  as Integer(Brief="Integer dimension", Default=3);
        Scalar as Real(Default=1);

        VARIABLES
        Var1(Dimension,Dimension) as Real(Default=1);
        Var2(Dimension,Dimension) as Real(Default=1);

        EQUATIONS

        diff(Var1(1,1))*'s' = Var2(1,2);
        
        for i in [2 : Dimension-1]
                diff(Var1(1,i))*'s' = Var2(1,i-1)+Var2(1,i+1);
        end

        diff(Var1(1,Dimension))*'s' = Var2(1,Dimension);
        diff(Var1(2:Dimension,:))*'s' = Scalar*Var2(2:Dimension,:);

        exp(Scalar)*Var2 = sin(time*'rad/s');
end

FlowSheet SampleArrays

        DEVICES
        dev1 as Arrays;

        INITIAL
        dev1.Var1 = -1;

        OPTIONS
        TimeStep = 0.1;
        TimeEnd = 2*3.1415;
end