First, a 'soapbox' about the quadrature method. It's perfectly fine and appropriate if the input uncertainties represent some number of standard deviations (the same number in each case) of roughly normally and independently distributed measurements. The answer then represents the same number of s.d. of the output.
But it's not unusual in engineering and in lab work that the input uncertainties would be better modeled as hard limits of uniform distributions. If there are many such inputs, the best approach is probably to turn each of those into a specific number of standard deviations, use quadrature, and understand that whatever number of s.d. you chose applies to the answer. But with only two inputs, you might be better served simply considering the four extreme values generated by plugging in the extremes of the two input ranges. The resulting range won't represent a uniform distribution, more triangular perhaps, but it will give you hard limits, and in engineering that may be more relevant.
Wrt your working so far, the formula you used for derivative of arctan assumes the angle is in radians. If it's in degrees you'll need to make a correction for that. It should be dimensionless - because angles are dimensionless.