lam0=0;
lam=0;
a0=[(42+lam0)  0
   0   (20+lam0)  ];
g=[-2;-20];
p0=a0\g;
delta0=norm(p0);
for i=1:1000,
a=[(42+lam(i))  0
   0   (20+lam(i))  ];
p(:,i)=a\g;
delta(i)=norm(p(:,i));
lam(i+1)=lam(i)+.1;
delta(i+1)=0;
end
X=p(1,:);Y=p(2,:);
plot(X,Y,'y-',X,Y,'go')
max(X.*X+Y.*Y)
pause
plot(lam,delta,'y-',lam,delta,'go')
