setlmis( [] ) ;
p11 = lmivar( 1 , [ 1 , 0 ] ) ;
p22 = lmivar( 1 , [ 1 , 0 ] ) ;
p33 = lmivar( 1 , [ 1 , 0 ] ) ;
p44 = lmivar( 1 , [ 1 , 0 ] ) ;
p55 = lmivar( 1 , [ 1 , 0 ] ) ;
p66 = lmivar( 1 , [ 1 , 0 ] ) ;
lmiterm( [ 1 1 1 p11 ] , -1 , 1 ) ;
lmiterm( [ 1 1 4 p11] , -1 , 1 ) ;
lmiterm( [ 1 1 4 p22 ] , -1 , 1 ) ;
lmiterm( [ 1 2 2 p22 ] , -1 , 1 ) ;
lmiterm( [ 1 3 3 p33 ] , -1 , 1 ) ;
lmiterm( [ 1 4 4 p33 ] , -1 , 1 ) ;
lmiterm( [ 1 4 4 p22 ] , 1 , 1 ) ;
lmiterm( [ 1 5 5 p22 ] , -1 , 1 ) ;
lmiterm( [ 1 6 6 p33 ] , 1 , 1 ) ;
lmiterm( [ 1 6 6 p11 ] , -1 , 1 ) ;
lmiterm( [ 1 7 7 p22 ] , -1 , 1 ) ;
lmiterm( [ 1 7 7 p33 ] , -2 , 1 ) ;
lmiterm( [ 1 7 7 p11 ] , 1 , 1 ) ;
Lmis= getlmis ;
[ tmin feas ] = feasp( Lmis ) ;
p1