I want to obtain the distance and angle between two coordinates and the arctan of the difference between these two coordinates

function [ vDist, vBearing] = vDistB( A , B )
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here

A = [E1,N1];
B = [E2,N2];

vDist = sqrt((E2-E1).^2)+(N2-N1).^2))

vBearing = atan((E2-E1)./(N2-N1))

[Walter Roberson]

We recommend atan2() instead of atan().

Your code uses E1, E2, N1, and N2, without defining them, and for unknown reasons it accepts A and B but ignores their values.