Please change in the following line, having permisible matrix dimention to works perfectly
Cep=k*kron(A,A)+2*M*(B-1/3*(kron(A,A))-kron(normal,normal));
Lets decompose to individual section (2nd Part with Bracket):
>> whos B
Name Size Bytes Class Attributes
B 6x6 288 double
Next
>> 1/3*(kron(A,A));
>> whos ans
Name Size Bytes Class Attributes
ans 36x1 288 double
Next
>> kron(normal,normal);
>> whos ans
Name Size Bytes Class Attributes
ans 9x9 648 double
Lets example matrix A(6x6)-B(36x1)-C(9x9), is that permisible? To addtion or subtraction, all must have same sizes. Please check the code and modify accordingly.
:)
