using Symbolics
using TyMath
using TySymbolicMath
R = 0.90
D1 = 238
D1a = 19.04
N1 = 120000
N1a = 12000
zr = quantile(Normal(), R)
eqs=[D1-N1/A-zr*sqrt(D1a^2+(N1a/A)^2)~0]
solve_for(eqs,[A])
r = sqrt(A / π)
println("r = ", r)

这个显示错误是哪里有问题