I have to solve a differential equation involving the error function
Erf, but my code is not able to evaluate it in the range I need.
I think I can formulate my problem in the following simplified terms: let's consider the following functions
VV[x_]:=Exp[-2 x^2] cc[x_]:=Erf[x]-1 f[r_]:=cc[r^-2]/Sqrt[cc[r^-2]^2 + (1 + r) VV[r^-2]^2]
I have to evaluate the function
f[r] for small values of
r. However, for
r<.2 I get the following message
Power::infy: "Infinite expression 1/0. encountered"
due to the presence of error functions. I tried to increase the precision but without success.
Do you have any suggestions to overcome this problem? Thank you very much for your help.