Advanced Search

Question: Error when executing procedure

Question:Error when executing procedure

GambiaMan 15 Maple
ode differential-equation numeric odeplot
November 04 2015
1

Hello,

I'm writing a code and I seem to have an issue when trying to implement a procedure. Here is the code:

 

> 

with(plots):
> 

Z := 75; A := 189; k := 14.6; Rm := 8*R; r0 := 10^(-8)*R; c := 137.036; ms := 105.66/(.51100)
> 

fmtoau := 10^(-15)/(0.529177e-1*10^(-9)):
> 

R := 1.1*fmtoau*A^(1.0/(3.0)):
> 

f := proc (x) options operator, arrow; 1/(1+exp(k*(x-1))) end proc:
> 

n0 := 3*Z*k^2/(4*Pi*(Pi^2+k^2)*R^3):
> 

n := proc (r) options operator, arrow; 4*Pi*n0*f(r/R) end proc:
> 

int(r^2*n(r), r = 0 .. Rm);

73.00000007

 (1)
> 

plot(n/n0, 0 .. 2*R);

 
> 

v1 := unapply(int(x^2*f(x), x), x):
> 

Vfermi := proc (r) options operator, arrow; -4*Pi*n0*R^2*(R*(v1(r/R)-v10)/r+v2Rm-v2(r/R)) end proc:
> 

Vuniform := proc (r) options operator, arrow; piecewise(r < R, -Z*(3/2-(1/2)*r^2/R^2)/R, -Z/r) end proc:
> 

plot([Vuniform(r), Vfermi(r), -Z/r], r = r0 .. 2*R, V = 1.2*Vfermi(r0) .. 0, legend = ["uniform charge", "Fermi distribution", "point charge"]);

 
> 

``
> 

plotsol1s := proc (E, K, r0, S, col) local Eqns, ICs, fnl, gnl, r, soln; global ms, c; Eqns := diff(fnl(r), r) = gnl(r)*[E+2*ms*c^2-Vuniform(r)]/c-(K+1)*fnl(r)/r, diff(gnl(r), r) = -fnl(r)*[E-Vuniform(r)]/c-(1-K)*fnl(r)/r; ICs := fnl(r0) = 1, gnl(r0) = 0; soln := dsolve({Eqns, ICs}, numeric); plots:-odeplot(soln, [r, fnl(r)], r0 .. S, color = col) end proc:
> 

plotsol1s(-3*10^5, -1, 10^(-10), Rm, red)

Error, (in f) unable to store '[HFloat(0.0)]' when datatype=float[8]

  
> 

``
> 

``


Any help would be greatly appreciated.

 

Gambia Man

 
Please Wait...