Skip to main content
Tweeted twitter.com/StackMma/status/1177191301000892417
edited tags
Link
xzczd
  • 68.4k
  • 9
  • 174
  • 489
edited tags
Link
bbgodfrey
  • 62.1k
  • 18
  • 92
  • 160
deleted 69 characters in body
Source Link
xzczd
  • 68.4k
  • 9
  • 174
  • 489

This is the first time I ask a question. I have seen many solutions to this error and tried but they are not working.

Here is the code.

Needs["DifferentialEquations`NDSolveProblems`"];
 
Needs["DifferentialEquations`NDSolveUtilities`"];
 
Needs["FunctionApproximations`"];

Clear["Global`*"];

e = 1.60217662*10^-19;
 
c = 299792458;
 
m = 9.10938356*10^-31;
 
ϵ0 = 8.854187817*10^-12;
 
A = 10^17;
 
h = 6.62607004*10^-34; hbar = 1.0545718*10^-34;
 k = 
 1/ϵ0 (2 e/h^3) (4 Pi/3);
 
 x0 = (3/(4 *Pi)*
     A*1*A*1.60439572*10^-24/(4*10^14))^(1/3)/100;

 xf = 10^0;

s = NDSolve[{1/x D[x y[x], {x, 2}] == 
    k ((e y[x]/c)^2 + 2 m*e y[x])^1.5, y[x0] == 10^6, y[xf] == 10^-5},
   y[x], {x, x0, xf}]

xfxf should be infinity but NDSolveNDSolve doesn't allow infinity so I put a relatively large number.

y[xf] is also needed to be zero at infinity. I put a relatively small number since I expected that y[x] is close to zero when x is approaching infinity.

Here is the error:

At x == 0.112455460174063`, step size is effectively zero; singularity or stiff system suspected. >>

This is the first time I ask a question. I have seen many solutions to this error and tried but they are not working.

Here is the code.

Needs["DifferentialEquations`NDSolveProblems`"];
 
Needs["DifferentialEquations`NDSolveUtilities`"];
 
Needs["FunctionApproximations`"];

Clear["Global`*"];

e = 1.60217662*10^-19;
 
c = 299792458;
 
m = 9.10938356*10^-31;
 
ϵ0 = 8.854187817*10^-12;
 
A = 10^17;
 
h = 6.62607004*10^-34; hbar = 1.0545718*10^-34;
 k = 
 1/ϵ0 (2 e/h^3) (4 Pi/3);
 
 x0 = (3/(4 *Pi)*
     A*1.60439572*10^-24/(4*10^14))^(1/3)/100;

 xf = 10^0;

s = NDSolve[{1/x D[x y[x], {x, 2}] == 
    k ((e y[x]/c)^2 + 2 m*e y[x])^1.5, y[x0] == 10^6, y[xf] == 10^-5},
   y[x], {x, x0, xf}]

xf should be infinity but NDSolve doesn't allow infinity so I put a relatively large number.

y[xf] is also needed to be zero at infinity. I put a relatively small number since I expected that y[x] is close to zero when x is approaching infinity.

Here is the error:

At x == 0.112455460174063`, step size is effectively zero; singularity or stiff system suspected. >>

This is the first time I ask a question. I have seen many solutions to this error and tried but they are not working.

Here is the code.

Needs["DifferentialEquations`NDSolveProblems`"];
Needs["DifferentialEquations`NDSolveUtilities`"];
Needs["FunctionApproximations`"];

Clear["Global`*"];

e = 1.60217662*10^-19;
c = 299792458;
m = 9.10938356*10^-31;
ϵ0 = 8.854187817*10^-12;
A = 10^17;
h = 6.62607004*10^-34; hbar = 1.0545718*10^-34;
k = 1/ϵ0 (2 e/h^3) (4 Pi/3);
x0 = (3/(4 *Pi)*A*1.60439572*10^-24/(4*10^14))^(1/3)/100;

xf = 10^0;

s = NDSolve[{1/x D[x y[x], {x, 2}] == 
    k ((e y[x]/c)^2 + 2 m*e y[x])^1.5, y[x0] == 10^6, y[xf] == 10^-5},
   y[x], {x, x0, xf}]

xf should be infinity but NDSolve doesn't allow infinity so I put a relatively large number.

y[xf] is also needed to be zero at infinity. I put a relatively small number since I expected that y[x] is close to zero when x is approaching infinity.

Here is the error:

At x == 0.112455460174063`, step size is effectively zero; singularity or stiff system suspected. >>

added 113 characters in body
Source Link
xzczd
  • 68.4k
  • 9
  • 174
  • 489
Loading
Source Link
Alex Leung
  • 91
  • 1
  • 1
  • 2
Loading