% Define a square Tip = []; x = 20; y = -20 z = 20; Q = InverseRhino([x,y,z]) Tip = Rhino(Q, zeros(4,1)); for y=-20:0.1:20 Tip = [Tip, [x, y, z,1]']; end for z = 20:0.1:60 Tip = [Tip, [x, y, z,1]']; end for y=20:-0.1:-20 Tip = [Tip, [x,y,z,1]']; end for z=60:-0.1:20 Tip = [Tip, [x,y,z,1]']; end Q = zeros(5,length(Tip)); for i=1:length(Q) Q(:,i) = InverseRhino(Tip(:,i)); end for i=1:length(Q) T = Rhino(Q(:,i), Tip); %TIP = [TIP,T]; end