code:
 class Cone(object):
    """docstring for Cone"""
    def __init__(self, v,d,theta):
        v.mustBePoint()
        d.mustBeVector()
        self.point = v
        self.vector = d
        self.theta = theta
    def intersectionTime(self,ray):
        """docstring for intersectionTime"""
        u=ray.point-self.point
        a=ray.vector.dot(self.vector)**2-self.vector.magnitude()**2*cos(self.theta)**2*ray.vector.magnitude()**2
        b=2*(ray.vector.dot(self.vector)*u.dot(self.vector) - u.dot(ray.vector)*self.vector.magnitude()**2*cos(self.theta)**2)
        c=u.dot(ray.vector)*self.vector.magnitude()**2*cos(self.theta)**2
        d=b**2-4*a*c
        if d<0:
            t=None
        else:
            t=[(-b-sqrt(d))/2./a,(-b+sqrt(d))/2./a]
            p=[self.intersectionPoint(x,ray) for x in t]
            pt=dict(zip(p,t))
            t=[pt[x] for x in filter(lambda x:self.inrange(x),p)]
            if len(t)==0:
                t=None
            else:
                t=min(t)
        return t
    def intersectionPoint(self,t,ray):
        if t==None:
            p=None
        else:
            # print t
            p=ray.point+ray.vector.scale(t)
        return p
    def inrange(self,p):
        """docstring for inrange"""
        u=p-self.point
        return u.magnitude()*cos(self.theta)<=self.vector.magnitude()
    def normalAt(self,p):
        """docstring for normalAt"""
        u=p-self.point
        return u.cross(u.cross(self.vector))
 

В ответ на:

---

те косяк алгебраический...
Поможите, - где косяк в алгебре?

---

Quote:

---

есть такой метод Рисовал, да - по шагам сложно. Вроде формулы правильные а получается фигня В смысле программинг тут наверное ни при чем - я походу на бумаге неправильно решил. Забыл уж этот линал нафиг.

---

code:
[pt[x] for x in filter(lambda x:self.inrange(x),p)]