; ; ; !y.range=1 close,1 openr,1,'moonlimbs' limbx=fltarr(79) limby=fltarr(79) x=findgen(79) y=findgen(79) for ii=0,78 do begin readf,1,xl,yl print,xl,yl limbx(ii)=xl limby(ii)=yl end plot,limbx,psym=1 xcoeff = poly_fit(x,limbx,1,cfit) oplot,cfit stop plot,limby,psym=1 ycoeff = poly_fit(y,limby,1,cfit) oplot,cfit stop end