A = [0 1 0; 0 0 1; -3000 -650 -45];
P = sdpvar(3,3); #defining sdp variables
gamma = sdpvar(1,1);
M = [-P*gamma P; P M*A+A'P];
constraints = [M<0, P>0]
solvesdp(constraints)