How to express the following inequality in form of LMI in Matlab, and use the feasp function ?

This is my LMI constraint :

(P1) + (beta^2) * R * R' > 0

where

P1 = A'(I + FC)'P + P(I+F*C)*A + A'*C'G'(Py) + (Py)GCA - C'(Pk)' - (Pk)*C + 2I

R = P + PFC + (Py)GC

beta = some constant

A' means the transpose of A matrix

This is a LMI in variables P, Py and Pk.

The problem is to find P, Py and Pk such that beta is maximized.

How to express the above LMI in Matlab ? I need to find the solution by using the feasp function too, for finding a solution to the LMI feasibility problem ?

The usage requires lmivar to express the inequality and the feasp function thereafter to compute the solution.