MATLAB外点惩罚函数法不知道那错了

function f=obj(x)
f=(x(1)-5)^2+(x(2)-6)^2+7;
function fai=objfunout(x)
%辅助函数
global r
fai=obj(x)+r*([2x1-x2-1;0])^2+r([x1+x2-2;0])^2;
function y=fun2pen(x)
%罚项函数
y=(max(0,1-x(1)))^2;
function penitentout
% 外点法罚函数方法
% fun2pen 约束函数
% objfun 目标函数
clc
global r
x0=[0;1];
r=0.3; c=5; ep=1e-5; k=1;
while rfun2pen(x0)>=ep
[x0,f]=fminsearch(@objfunout,x0)
r=cr; k=k+1;
end
disp('最优解x0='),disp(x0)
disp('k='),disp(k)
disp('目标函数f(x)='),disp(obj(x0))