解一元二次方程

syms erfa H;

u0 = 961.97416;
v0 = 687.53441;
fx = 1857.481736;
fy = 1872.211946;

u = 487;
v = 992;
Y = 5;
X = -0.92;

%[erfa,H] = solve('H = Y * ( tan(erfa - atan( (v-v0)/fy ) ) )',' (u-u0)*H = X * sqrt(fx*fx+(v-v0)*(v-v0)) * sin(erfa - atan((v-v0)/fy)) ')

[erfa,H] = solve('H = 5 * ( tan(erfa - atan( (992-687.53441)/1872.211946 ) ) )',' (487-961.97416)*H = -0.92 * sqrt(1857.481736*1857.481736+(992-687.53441)*(992-687.53441)) * sin(erfa - atan((992-687.53441)/1872.211946)) ')