def draw_points(self, area, a, b):
        if len(a) < 3:
            return
        x1 = a[-1]
        y1 = b[-1]
        w, h = area.GetSize()
        a, b, ma, mb = self.scale_points(a, b, w, h)
        c = zip(a, b)
        c = filter_points.remove_dupes(c)
        
        area.Clear()
        area.SetPen(wx.Pen((0,0,0)))
       
        area.DrawPoints(c)

##        try:
##            rc = c#random.sample(c, 1000)
##            ra = [ x for x, y in rc ]
##            rb = [ y for x, y in rc ]
##
##            #a, b, c = matfunc.polyfit((ra, rb))
##            #lines = self.get_lines(w, a, b, c)
##            #a, b, c, d = matfunc.polyfit((ra, rb), 3)
##            a, b = matfunc.polyfit((ra, rb), 1)
##            lines = self.get_lines1(w, a, b)
##            area.DrawLines(lines)
##            #area.draw_lines(area.gc, lines)
##        except Exception, e:# ZeroDivisionError:
##            print e
##            pass


        u = SizedList(10)
        t = SizedList(10)
        points = []
        for (x,y) in c:
            u.append(x)
            t.append(y)

            #x2 = float(sum(u)) / float(len(u))
            x2 = x
            y2 = h - filter_points.standard_deviation(t)
            #y2 = float(sum(t)) / float(len(t))
            points.append((x2,y2))
                    
        #area.draw_lines(area.gc, points)
        area.SetPen(wx.Pen((255,0,0)))
        area.DrawPoints(points)

        def avg_atan2_p(a, b):
            atan = 0
            x = 0
            y = 0
            as = 0
            t = float(len(a))
            for p1, p2 in zip(a, b):
                atan += math.atan2(p1, p2) / t
                try:
                    as += math.sqrt(p1**2 + p2**2) / t
                except ZeroDivisionError, e:
                    print e
        def median(a, b):
            a = list(a)
            a.sort()
            b = list(b)
            b.sort()
            ap = min(len(a)-1, 5)
            bp = min(len(b)-1, 5)
            return a[ap], b[bp]           
            

        u = SizedList(10)
        t = SizedList(10)
        points = []
        for (x,y) in c:
            u.append(x)
            t.append(y)

            x2, y2 = avg_atan2_p(u, t)
            #x2, y2 = median(u,t)
            #x2 = x
            #y2 = float(sum(t)) / float(len(t))
            points.append((x2,y2))
                    
        #area.draw_lines(area.gc, points)
        area.SetPen(wx.Pen((0,0,255)))
        #area.DrawPoints(points)        

##        ##############################
##        u = SizedList(10)
##        t = SizedList(10)