Changeset 485


Ignore:
Timestamp:
Mar 14, 2008, 10:56:48 AM (14 years ago)
Author:
Rodolfo Rodrigues
Message:

Updated sample/miscellaneous/sample_matrix

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/sample/miscellaneous/sample_matrix.mso

    r247 r485  
    1818* Author: Rodolfo Rodrigues
    1919* $Id: Matrix_Function.mso 89 2006-12-14 20:29:34Z bicca $
    20 *--------------------------------------------------------------------*#
     20*------------------------------------------------------------------*#
    2121#
    2222#               EMSO sample file of matrix functions
     
    3737* Given the matrix A:
    3838*
    39 *       |-1  -2   3  1  0|
    40 *       | 0  -1  -1  0  2|
    41 *       | 0   0   5  1  0|
     39*       |-2   3   1   0|
     40*       |-1  -1   0   2|
     41*       | 0   5   1   0|
    4242*
    4343*--------------------------------------------------*#
     
    4848        N                       as Integer;
    4949        A(M,N)          as Real;
     50       
    5051       
    5152        VARIABLES
     
    5859        Transp_(N,M)as Real;
    5960       
     61       
    6062        SET
    61         M = 5; # columns
    62         N = 3; # rows
     63        M = 3; # rows
     64        N = 4; # columns
    6365       
    64         A(:,1) = [-1, -2,  3, 1, 0];
    65         A(:,2) = [ 0, -1, -1, 0, 2];
    66         A(:,3) = [ 0,  0,  5, 1, 0];
     66        A(:,1) = [-2, -1,  0];
     67        A(:,2) = [ 3, -1,  5];
     68        A(:,3) = [ 1,  0,  1];
     69        A(:,4) = [ 0,  2,  0];
     70       
    6771       
    6872        EQUATIONS
     73#*---------------------------------------------------
     74*
     75*          |-2   3   1   0|   
     76*       sum|-1  -1   0   2| = |-3  7  2  2|
     77*          | 0   5   1   0|
     78*
     79*---------------------------------------------------*#
    6980        "Function sum"
    7081        Sum_    = sum(A);
    7182       
     83
     84#*---------------------------------------------------
     85*
     86*           |-2   3   1   0|   | 2|
     87*       sumt|-1  -1   0   2| = | 0|
     88*           | 0   5   1   0|   | 6|
     89*
     90*---------------------------------------------------*#
    7291        "Function sumt"
    7392        Sumt_   = sumt(A);
    7493       
     94       
     95#*---------------------------------------------------
     96*
     97*           |-2   3   1   0|
     98*       prod|-1  -1   0   2| = | 0  -15  0  0|
     99*           | 0   5   1   0|
     100*
     101*---------------------------------------------------*# 
    75102        "Function prod"
    76103        Prod_   = prod(A);
    77104       
     105       
     106#*---------------------------------------------------
     107*
     108*            |-2   3   1   0|   | 0|
     109*       prodt|-1  -1   0   2| = | 0|
     110*            | 0   5   1   0|   | 0|                                           
     111*
     112*---------------------------------------------------*# 
    78113        "Function prod"
    79114        Prodt_  = prodt(A);
    80 
     115       
     116       
     117#*---------------------------------------------------
     118*
     119*             |-2   3   1   0|   |-2 -1  0|
     120*       transp|-1  -1   0   2| = | 3 -1  5|
     121*             | 0   5   1   0|   | 1  0  1|
     122*                                                        | 0  2  0|
     123*
     124*---------------------------------------------------*#
    81125        "Function transp"
    82126        Transp_ = transp(A);
     127       
    83128       
    84129        OPTIONS
Note: See TracChangeset for help on using the changeset viewer.