Как найти площадь треугольника методом Монте-Карло?

# -*- coding: utf-8 -*- 
import math
import random
 
e = 0
n = 7 # Степень числа N
t = 0.01 #  Точность сравнения вещественных чисел
N = 0 # Колличество точек
M = 0 # Колличество точек попавщих в треугольник
 
cx = input("Введите координату точки c по X:")
cy = input("Введите координату точки c по Y:")
bx = input("Введите координату точки b по X:")
by = input("Введите координату точки b по Y:")
ex = input("Введите координату точки e по X:")
ey = input("Введите координату точки e по Y:")
 
 
def PloshPr(ex,ey): #Вычисление площади прямоугольника 
    PloshPr = math.sqrt(ey**2) * math.sqrt(ex**2) 
    print "Ploshad Prjam", PloshPr
    return PloshPr
    
    
def PloshTr(cx,cy,bx,by,ex,ey): #Вычисление площади прямоугольника
    PloshTr = abs((cx-bx)*(ey-by)-(cy-by)*(ex-bx))/2.0
    return PloshTr
    
    
def  Proverka(M,N,n,ex,ey,t,e):
        for i in range(3, n+1):     
            N = 10**i           
            touch = [(random.uniform(0,ex+1),random.uniform(0,ey+1)) for i in range(10**i)] # Заполнение списка рандомными координатами   
            for (x,y) in touch:
                Stouch =  PloshTr(x,y,bx,by,ex,ey)+PloshTr(cx,cy,x,y,ex,ey)+PloshTr(cx,cy,bx,by,x,y) #Проверка попала ли точка в треугольник 
                if abs(PloshTr(cx,cy,bx,by,ex,ey) - Stouch) < t:
                    M = M+1
                else:
                    e = e+1
        PlTr = PloshPr(ex,ey)*M/N # Вычисление площади треугольника
        print "N=",N
        print "E", e
        print "M", M
        print PlTr
 
PloshPr(ex,ey)
PloshTr(cx,cy,bx,by,ex,ey)      
Proverka(M,N,n,ex,ey,t,e)