解等式约束信赖域子问题Matlab程序

%  clear;clc;
%  B=[510;130;00-1]
%  g=[111]‘   %min x‘Bx+2g‘x
%  delta= 0.6   %s.t. norm（s）^2=delta^2
 function[aabb]=testrs（Bgdelta）
epsi=10^（-10）;
[mn]=size（B）;
s=zeros（m1）;
I=diag（eye（m））;
k=0;

[UDB]=eig（B）;%将B对角化，那g也应该相应对角化吧？
Dg=U‘*g;

%变成向量
DDB=diag（DB）;

lambda=-DDB（1）;

%主要求解（B+lambda*I）*s=-g
DLag=DDB+lambda*I;

if abs（DLag（1））epsi
    k=1;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%无解
else
   for i=1:m
       if abs（DLag（i））            s（i）=0;
            c=i;
        else 
            break
        end
    end
    for i=c+1:m
        s（i）=-Dg（i）./DLag（i）;
    end
 
    n