MBN ex 2

MBNex2.m

pre = [0 1 0 0 1;
    1 0 0 0 0;
    0 1 0 0 0;
    0 0 1 0 0;
    0 0 0 2 0;
    0 0 0 0 1;
    0 0 1 0 0;
    0 0 0 1 0];

post = [1 0 0 0 0;
    0 1 0 0 1;
    0 1 1 0 0;
    0 0 1 1 0 ;
    0 0 0 0 1;
    0 0 0 2 0;
    0 0 0 0 0;
    0 0 0 0 0];


A = post - pre

S = A';

r = [5 0 4 10 0 0 0 3]'

M = [0 1 2 10 12]'

disp('Transition M* = M + S*r')
M + S*r

y = [1 1 0 0 0]'

disp('P-invariance: A*y=0')
A*y

x1 = [1 1 0 0 0 0 0 0]'
x2 = [0 0 0 0 1 1 0 0]'

disp('S*x1')
S*x1
disp('S*x2')
S*x2