From a math standpoint, the angle is arctan((y2-y1)/(x2-x1)). I think that to implement that on the robot, you'll need either the math option or to use Karel.
I've had some success using Taylor expansion to get around the math limitations on base fanucs. You have to scale the function though, to account for the fact that a real in the registers is cut off after 6 decimal places.