/* Neznáme koeficienty sú v nasledovnom poradí: */
/* a_0,0  a_0,1  a_0,2  a_0,3  a_1,0  a_1,1  a_1,2  a_1,3 */
A: matrix([1, 0, 0, 0, 0, 0, 0, 0],     /* funkcia prechádza P_0 */ 
          [0, 0, 0, 0, 1, 0, 0, 0],     /* funkcia prechádza P_1 */ 
          [0, 0, 0, 0, 1, 2, 4, 8],     /* funkcia prechádza P_2 */ 
          [1, 3, 9,27,-1, 0, 0, 0],     /* C^0 spojitosť */ 
          [0, 1, 6,27, 0,-1, 0, 0],     /* C^1 spojitosť */ 
          [0, 0, 2,18, 0, 0,-2, 0],     /* C^2 spojitosť */ 
          [0, 0, 2, 0, 0, 0, 0, 0],     /* prirodzený splajn -- p''_0(1) = 0 */ 
          [0, 0, 0, 0, 0, 0, 2,12]);    /* prirodzený splajn -- p''_1(6) = 0 */

/* Pravá strana sústavy */
b: transpose(matrix([1, 3, 1, 0, 0, 0, 0, 0]));

/* Riešenie sústavy */
x: invert(A).b;