3)
                 y[1]   y[2]
          ------------------------
          |      |       |       |  
          |      |       |       |
          |      | --->  |       |
      v[2]|------|-------|-------| o[2]
          |    P12 \    P22      |  
          |      |  \    |       |
          |      |    \  |       |
      v[1]|------|-------|-------| o[1]
          |    P11 --> P21       |  
          |      |       |       |
          |      |       |       |
          |------|-------|-------|
                a[1]    a[2]


 P11:   4*u[1,1]=a[1]+u[2,1]+u[1,2]+v[1]
 P21:   4*u[2,1]=a[2]+o[1]+u[2,2]+u[1,1]
 P12:   4*u[1,2]=u[1,1]+u[2,2]+y[1]+v[2]
 P22:   4*u[2,2]=u[2,1]+o[2]+y[2]+u[1,2]
eli (Maplen output-muoto):

                  4 u[1, 1] = a[1] + u[2, 1] + u[1, 2] + v[1]

                  4 u[2, 1] = a[2] + o[1] + u[2, 2] + u[1, 1]

                  4 u[1, 2] = u[1, 1] + u[2, 2] + y[1] + v[2]

                  4 u[2, 2] = u[2, 1] + o[2] + y[2] + u[1, 2]

Jotta matriisi näkyisi selvästi, merkitään: 
sijoitukset := {u[1, 1] = x1,u[2, 1] = x2, u[1, 2] = x3, u[2, 2] = x4}
yht1:=subs(sijoitukset,yht1); # ja samoin muut =>
     4 x1 = a[1] + x2 + x3 + v[1]
     4 x2 = a[2] + o[1] + x4 + x1
     4 x3 = x1 + x4 + y[1] + v[2]
     4 x4 = x2 + o[2] + y[2] + x3

Tästä näkyykin matriisi (siirretään termejä):

A:=matrix([[-4,1,1,0],[1,-4,0,1],[1,0,-4,1],[0,1,1,-4]]);
                               [-4     1     1     0]
                               [                    ]
                               [ 1    -4     0     1]
                          A := [                    ]
                               [ 1     0    -4     1]
                               [                    ]
                               [ 0     1     1    -4]

Samoin oikea puoli:



b:=vector([-a[1]-v[1],-a[2]-o[1],-y[1]-v[2],-o[2]-y[2]]);
       b := [-a[1] - v[1], -a[2] - o[1], -y[1] - v[2], -o[2] - y[2]]

b:=vector([-cos(Pi/3)-0,-cos(2*Pi/3)-0,-sin(Pi/3)-0,-0-sin(2*Pi/3)]);

                         [                  1/2         1/2]
                    b := [-1/2, 1/2, - 1/2 3   , - 1/2 3   ]

Nyt vain ratkaistaan yhtälösysteemi  A x = b

linsolve(A,b);
[             1/2                 1/2               1/2                 1/2]
  [5/48 + 1/16 3   , - 5/48 + 1/16 3   , 1/48 + 3/16 3   , - 1/48 + 3/16 3   ]

> evalf(");
              [.2124198422, .0040865088, .3455928598, .3039261932]
                  
> matrix(2,2,");
                          [.2124198422    .0040865088]
                          [                          ]
                          [.3455928598    .3039261932]

> swaprow(",1,2);
                          [.3455928598    .3039261932]
                          [                          ]
                          [.2124198422    .0040865088]

Tässä ovat samassa järjestyksessä, kuin yllä olevassa kuvassa.