clear all;
RR=load('../obs_float/r_float_optim7.dat');
R1=RR(3:2:end,end);
R2=RR(4:2:end,end);
hold on
%figure;
plot(R1+1,R2+1,'k-');
hold on
plot(R1(1)+1,R2(1)+1,'*r');

obs=load('../obs_float/obs.dat')
plot(obs(1,end)+1,obs(2,end)+1,'*g');