求解数值x,y，方程及各参数如下，请问如何写程序求出x,y值？
c0=0.4
c01=0.295
c1=0.225
c12=0.095
c2=0.05

D0=2e-14
D1=2e-12
D2=2e-10
r0=sqrt(D0/D1)
r1=sqrt(D1/D2)

f1=(c0-c01)/(c01-c1)
f2=exp(-x^2)/(2-erfc(y))
f3=(c1-c12)/(c01-c1)
f4=exp(-(x*r0)^2)/(erfc(x*r0)-erfc(y))
f5=(c1-c12)/(c12-c2)
f6=exp(-y^2)/(erfc(x*r0)-erfc(y))
f7=c2/(c12-c2)
f8=exp(-(y*r1)^2)/erfc(y*r1)

f1*f2-1/r0*f3*f4-x*sqrt(pi)=0
f5*f6-1/r1*f7*f8-y*sqrt(pi)=0

求解数值x,y，方程及各参数如下，请问如何写程序求出x,y值？
c0=0.4
c01=0.295
c1=0.225
c12=0.095
c2=0.05

D0=2e-14
D1=2e-12
D2=2e-10
r0=sqrt(D0/D1)
r1=sqrt(D1/D2)

f1=(c0-c01)/(c01-c1)
f2=exp(-x^2)/(2-erfc(y))
f3=(c1-c12)/(c01-c1)
f4=exp(-(x*r0)^2)/(erfc(x*r0)-erfc(y))
f5=(c1-c12)/(c12-c2)
f6=exp(-y^2)/(erfc(x*r0)-erfc(y))
f7=c2/(c12-c2)
f8=exp(-(y*r1)^2)/erfc(y*r1)

f1*f2-1/r0*f3*f4-x*sqrt(pi)=0
f5*f6-1/r1*f7*f8-y*sqrt(pi)=0