I should stress that I know mathematica. I had tried all possibly variaions before wrote here. If anyone wants help me please first of all run the code. I think mathematica or the article has a bug. Please nobody reply the message without running the code.
You are right. But after I tried all variable ranges with correct code, nothing changed. Mathematica file is attached for anyone intersted in this plot. Thanks for help.
NascentOxygen and DrClaude, thank you for answers. I plotted this without elliptic integral. But this article referred to elliptic integrals. DrClaude, I am aware of your writings but mathematica has given error mesage: