How to find Intersection point of plane and vector/line in 3d

Hello! Could you please help me why my fucntion does not calculate intersection point I. It gaves me my initial X0, Y0, Z0 coordinate

I used this example

https://www.mathworks.com/matlabcentral/fileexchange/103835-plane-and-line-intersection

[I ] = planeLocation6(1.5, 3.3371, 0.1820)
function [ I ] = planeLocation6(XBar, YBar, ZBar)
X0 = 1.5;      
Y0 = 1.5;
Z0 = 3.0;
Theta0 = 30;

planes(:,:,1) = [0 3 3; 0 0 3; 0 3 0; 0 0 0; 0 0 0];
planes(:,:,2) = [0 0 3; 3 0 3; 0 0 0; 3 0 0; 0 0 0];
planes(:,:,3) = [3 0 3; 3 3 3; 3 0 0; 3 3 0; 3 0 0];
planes(:,:,4) = [3 3 3; 0 3 3; 3 3 0; 0 3 0; 0 3 3];
planes(:,:,5) = [0 3 0; 3 3 0; 0 0 0; 3 0 0; 0 0 0];
planes(:,:,6) = [0 3 3; 3 3 3; 0 0 3; 3 0 3; 0 0 3];
location_plane = 6;
for j=1:6
    j==location_plane
    plane = planes(:,:,j);
    p0 = plane(1,:);
    p1 = plane(2,:);
    p2 = plane(3,:);
    p3 = plane(4,:);
    V0 = plane(5,:);
    Ri = [X0 Y0 Z0];  
    Rr = [XBar YBar ZBar];
    u = Rr-Ri;
    w=Ri-V0;
    A = p0-p2;
    B = p0-p3;
    n=cross(A,B);      % normal vector of the Plane 
    u = Rr-Ri;
    w=Ri-V0;
    D=dot(n,u);
    N=-dot(n,w);
    sI = N / D;
    I = Ri+ sI.*u;
    disp(I);
end

Thank you in advance, BR Aknur

[Rik]

Your code doesn't actually call the function you linked. Is there a reason why you don't want to use it?

[Rik]

I recovered the removed content from the Google cache (something which anyone can do). Editing away your question is very rude. Someone spent time reading your question, understanding your issue, figuring out the solution, and writing an answer. Now you repay that kindness by ensuring that the next person with a similar question can't benefit from this answer.

[Aknur]

Hello, dear @Rik I have removed because in general my code was not correct and it can lead to misunderstanding for other viewers. As I used this function with mistakes https://www.mathworks.com/matlabcentral/fileexchange/103835-plane-and-line-intersection

Answer 1

[Bhavana Ravirala]

Hi Ankur,

I understand that you are unable to calculate the intersection point ‘I’.

When I try to run the code at my end, I am getting a value of ‘I’ as Infinity and not the initial values.

You are not getting the intersection point because the dot product of ‘n’ and ‘u’ is zero with the inputs which you have provided to the function. Try giving different inputs to the function such that ‘D’ does not become zero and ‘I’ does not be Infinity.

Hope this helps!!

[Aknur]

Dear, @Bhavana Ravirala thank you for your feedback. Yes D become zero, Can I ask if I want extract special data in array which will met the requirements of if statement x > 0 & x < 3 & y > 0 & y < 3 & z > 0 & z < 3 how I can do it? I used isnan and isinf to clean data and take out Nan and Infinite values. Now I'm thining how to take only that I which will meet the requirement mentioned above.