function [x r]=multiLinearBDURLS(x0,A,La,b,Lb)
%用多步線性迭代法求解一個非線性最小二乘問題
% (Regularized least squares problem with bounded data uncertainty)
% BDU RLS:
% J3=min max[||[A+e_A]*x-[b+e_b]||^2+p||x||^2]
% where e_A,e_b:the uncertainty,||e_A||<=La,||e_b||<=Lb
% J3=>J4=min[[||Ax-b||+La*||x||+Lb]^2+p||x||^2]
% solution:x=inv(A'*A+r*I)*A'*b
%          r=La*||Ax-b||/||x||+p*||Ax-b||/(||Ax-b||+La*||x||+Lb)
% INPUT:
% x0: the initial estimate
% A,b: data
% La,Lb: the upper bound of the uncertainty of A and b
% 

Niter=500;
x=x0;
p=1;
for Iiter=1:Niter
    norm1=(A*x)'*(A*x);
    norm2=
    r=La*/(x'*x)+p*
end