[debateshang]

clear all
clc
syms a b c d m n p q
A=[a b;c d];
C=[m n;p q];
syms x
% A=[1.2 0;0 1.3];
% C=[0.5 0;0 0.6];
k1=C*[exp(4*x-4) 0;0 exp(-x+1)]*A;
K=[5 0;0 -3]*k1.'*k1;
L=k1.'*k1*[5 0;0 -3];
E=K+L;
F=[exp(x) 0;0 exp(-2*x)]*[10;10];
y=F.'*E*F;
% solve(y)