aaa.asv

高 斯 迭 代

load molecule;
[smallI,smallJ,smalldist]=find(smallS);
subplot(1,2,1);
for i=1:length(smalldist)
    line([0,smalldist(i)],[i i]);
        text(-0.07,i,sprintf('%d',smallJ(i)));
        text(smalldist(i)+.02,i,sprintf('%d',smallI(i)));
        
    end
    axis([-.1 1 0 10],'off');
    title('Measured distances between Pairs');
    subplot(1,2,2);
    plot(smallsoln(:,1),smallsoln(:,2),'r.','MarkerSize',15);
   for i = 1:length(smalldist)
        line([smallsoln(smallI(i),1) smallsoln(smallJ(i),1)],...
         [smallsoln(smallI(i),2) smallsoln(smallJ(i),2)]);
 end
textpos=smallsoln+[-.09 0;-.08 .08;.04 .01;.04 0;.03 -.07;.03 .07];
    for i=1:6;
        text(textpos(i,1),textpos(i,2),sprintf('%d',i));
    end
    axis([-0.1 1 0 1.1],'off','square');
    title('Atom Locations');