% orthog - Program to test if a pair of vectors % is orthogonal. Assumes vectors are in 3D space clear all; help orthog; % Clear the memory and print header %* Initialize the vectors a and b a = input('Enter the first vector: '); b = input('Enter the second vector: '); %* Evaluate the dot product as sum over products of elements a_dot_b = 0; for i=1:3 a_dot_b = a_dot_b + a(i)*b(i); end %* Print dot product and state whether vectors are orthogonal if( a_dot_b == 0 ) disp('Vectors are orthogonal'); else disp('Vectors are NOT orthogonal'); fprintf('Dot product = %g \n',a_dot_b); end