坐标轮换法

function [xminf] = minconPS（fgx0deltauvareps1eps2）
if nargin == 7
    eps2 = 1.0e-6;
end
n = length（var）;
y = x0;
bmainCon = 1;

while bmainCon 
    yf = Funval（fvary）;
    yk_1 = y;
    for i=1:n
        tmpy = zeros（size（y））;
        tmpy（i） = delta（i）;
        tmpf = Funval（f vary+tmpy）;
        
        for j=1:length（g）
            cong（j） = Funval（g（j） vary+tmpy）;
        end
        if tmpf < yf && min（cong） >= 0
            bcon = 1;
            while bcon
                tmpy（i） = 2*tmpy（i）;
                tmpf_i = Funval（f vary+tmpy）;

                for j=1:length（g）
                    cong_i（j） = Funval（g（j） vary+tmpy）;
                end
                if tmpf_i =0
                    y_res = y  + tmpy;
                else
                    bcon = 0;
                end
            end
        else
            tmpy（i） = delta（i）;
            tmpf = Funval（f vary-tmpy）;