function DrawPlane(n,p,scale)

hold on
v1 = [1 0 0]';
v1 = v1 - (v1'*n)*n;
v1 = v1/norm(v1);
v2 = Vec2Skew(v1)*n;


for i = -10 : 10
    x1 = -10 * scale * v1 + i * scale * v2 + p;
    x2 = 10 * scale *v1 + i * scale * v2  + p;
    hold on
    plot3([x1(1) x2(1)], [x1(2) x2(2)], [x1(3) x2(3)], '-', 'Color', [0.5 0.5 0.5]);
end

for i = -10 : 10
    x1 = -10 * scale * v2 + i * scale * v1 + p;
    x2 = 10 * scale *v2 + i * scale * v1  + p;
    hold on
    plot3([x1(1) x2(1)], [x1(2) x2(2)], [x1(3) x2(3)], '-', 'Color', [0.5 0.5 0.5]);
end