To solve equation (2) using the value of P obtained from equation (3), you can find the maximum of the singular value decomposition of the P and E matrices and then check whether the condition in equation (2) is satisfied. Fo this you can use the svd function in MATLAB. You can use the below code snippet for your reference:
sigma_max_P = max(svd(P)); % maximum singular value of P
sigma_max_E = max(svd(E)); % maximum singular value of E
% Check if the condition in equation (2) is satisfied
if sigma_max_E < 1/sigma_max_P
disp('System is stable');
else
disp('System is not stable');
end
This snippet will determine whether the system remains stable based on the condition in equation (2). To know more about the svd function you can refer to its documentation here: https://in.mathworks.com/help/matlab/ref/double.svd.html