>> A = [1 1 0 0.5 ; 1 0 -1 0.5 ; 0 50 0 15 ; 0 0 40 15]

A =

    1.0000    1.0000         0    0.5000
    1.0000         0   -1.0000    0.5000
         0   50.0000         0   15.0000
         0         0   40.0000   15.0000

>> B = [-1 -1 0 0.5 ; 1 0 -1 0.5 ; 0 50 0 15 ; 0 0 40 15]

B =

   -1.0000   -1.0000         0    0.5000
    1.0000         0   -1.0000    0.5000
         0   50.0000         0   15.0000
         0         0   40.0000   15.0000

>> C = transpose([0,0,60,60])

C =

     0
     0
    60
    60

>> x1 = A ^ -1 * C

x1 =

    -2
     0
     0
     4

>> X2 = B ^ -1 * C

X2 =

    0.0896
    0.7164
    0.8955
    1.6119

>> 0.0896 .* 15

ans =

    1.3440

>> 1.6119 .* 15

ans =

   24.1785

>> 0.8955 .* 40

ans =

   35.8200

>> 0.7164 .* 50

ans =

   35.8200

>> 35.8200 .* 0.8955

ans =

   32.0768

>> 35.8200 .* 0.764

ans =

   27.3665

>> 0.896 .* 0.0896

ans =

    0.0803

>> 1.6119 .* 15

ans =

   24.1785

>> 24.1785 .* 1.6119

ans =

   38.9733

>> X8 = [1 -1 0 0 0 0 0 0 ; 0 1 -1 -1 0 0 0 0 ; 1 0 0 0 0 0 0 0 1 ; 0 0 0 0 1 0 -1 1 ; 0 0 0 1 -1 -1 0 0; -4 0 0 0 5 -5 20 0 ; -10 0 0 -5 -5 0 0 0 ; 4 0 -15 5 0 5 0 0]
Dimensions of arrays being concatenated are not consistent.

>> X8 = [1 -1 0 0 0 0 0 0 ; 0 1 -1 -1 0 0 0 0 ; 1 0 0 0 0 0 0 1 ; 0 0 0 0 1 0 -1 1 ; 0 0 0 1 -1 -1 0 0; -4 0 0 0 5 -5 20 0 ; -10 0 0 -5 -5 0 0 0 ; 4 0 -15 5 0 5 0 0]

X8 =

     1    -1     0     0     0     0     0     0
     0     1    -1    -1     0     0     0     0
     1     0     0     0     0     0     0     1
     0     0     0     0     1     0    -1     1
     0     0     0     1    -1    -1     0     0
    -4     0     0     0     5    -5    20     0
   -10     0     0    -5    -5     0     0     0
     4     0   -15     5     0     5     0     0

>> Y8 = transpose([-3 0 -3 0 0 -15 -10 0])

Y8 =

    -3
     0
    -3
     0
     0
   -15
   -10
     0

>> z8 = X8 ^ -1 * Y8

z8 =

   -0.5773
    2.4227
    0.6086
    1.8141
    1.3405
    0.4736
   -1.0822
   -2.4227

>> 15 + 20 .* -2.42

ans =

  -33.4000

>> 4 .* -0.58 + 5 .* -1.08

ans =

   -7.7200

>> -10 + 15 .* 1.81

ans =

   17.1500

>> function = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)));
 function = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)));
 ↑
Error: Function definition not supported in this context. Create functions in code
file.

>> function = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)))
 function = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)))
 ↑
Error: Function definition not supported in this context. Create functions in code
file.

>> function_shittot33 = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)));
>> integral (function_shittot33 , 1i , 1 ,"Waypoints",1+1i)

ans =

  -0.7103 + 3.0191i

>> fun2 = @(al) cos(al) + 1i * sin(al);
>> fun22 = @(al) -sin(al) + 1i * cos(al);
>> integral ( @(k) function_shittot33(fun2(t)) .* fun22(t) , 0.5 * pi , 0)
Undefined function or variable 't'.

Error in @(k)function_shittot33(fun2(t)).*fun22(t)

Error in integralCalc/iterateScalarValued (line 314)
                fx = FUN(t);

Error in integralCalc/vadapt (line 132)
            [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)
        [q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);

>> function_shittot33 = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)));
>> integral (function_shittot33 , 1i , 1 ,"Waypoints",1+1i)

ans =

  -0.7103 + 3.0191i

>> fun2 = @(al) cos(al) + 1i * sin(al);
>> fun22 = @(al) -sin(al) + 1i * cos(al);
>> integral ( @(k) function_shittot33(fun2(t)) .* fun22(t) , 0.5 * pi , 0)
Undefined function or variable 't'.

Error in @(k)function_shittot33(fun2(t)).*fun22(t)

Error in integralCalc/iterateScalarValued (line 314)
                fx = FUN(t);

Error in integralCalc/vadapt (line 132)
            [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)
        [q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);

>> function_shittot33 = @(z) log(1./ conj(z)) .* exp (sqrt (sin(z)));
>> integral (function_shittot33 , 1i , 1 ,"Waypoints",1+1i)

ans =

  -0.7103 + 3.0191i

>> fun2 = @(al) cos(al) + 1i * sin(al);
>> fun22 = @(al) -sin(al) + 1i * cos(al);
>> integral ( @(k) function_shittot33(fun2(k)) .* fun22(k) , 0.5 * pi , 0)

ans =

   0.1667 + 2.5673i

>>