5 years, 16 days

## MaplePrimes Activity

### These are replies submitted by michaelvio

I want to calculus the Time Quanta of the electron in the 1 S orbital of the Hydrogen atom thus I must have the exact value of the wave function of the electron!

You can see details at http://michaelvio.orgfree.com/T&SQ.pdf

electronFlux1.mw

## Thank 's...

electroSum.mw

I try to write the electron 1 S Hydrogen atom equation This is my site: http://michaelvio.orgfree.com/

Please contact me on my email at michaelvio@yahoo.com or michaelvio@gmail.com to give you all the details.

a=1/137 c= light speed E= Rythgerg energy rb is RadiusBohr and R is a distant constant that results from the Theory that is on the site.

Excellent solution!electronFlux1.mw

I want to consult with you on the calculus that I made ant for Time Quanta of the electron, please!

In the last file, I put the calculus of electron and I want to calculate the Time Quanta of the electron that must be rb/c the time between two consecutive o wave functions of the Hydrogen electron 1 S orbital with the boundary condition the velocity is c/137 all the time at radius R wave function =0 and

in origin, the derivate is R with the norm condition int(P(r,o),r=0..R)=1

## Change of var file...

The attached file is elecSVL4Ps.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/elecSVL4Ps.mw .

I look forward to your's opinion

## Change of Var...

```eq1 := 2*m*(E + 8*Pi*epsilon/r)*f(r, t)/h^2 + R*diff(f(r, t), r \$ 2)/r - diff(f(r, t), t \$ 2)/(a^2*c^2) = 0;
iv1 := f(rb, 0) = 0, f(R, t) = 0, D[1](f)(0, 0) = R, D[2](f)(rb,0)=a*c;
sys := [eq1, iv1]:
Sol := pdsolve(sys, HINT=`*`);
pt := pdetest(Sol, sys);
```

These are the correct initial value, it works with the first 3 conditions but when I introduced the fourth iv:D[2](f)(rb,0)=a*c the compilation time exceed very much. To solve it I used Change of Var :r=rb*rho;

Please help me find a solution that must not contain the _c[1] it must result from the 2 conditions.

I also want to find the exact value of f(r,t) with the condition f(rb, 0) = 0; and diff(f(rb,t),t)=a*c for the value t=0 and if is possible the pulsation of the sinusoidal solution of f(r,t). [the solution is a combination of AiryAi ; AiryBi and sinusoidal sin(a*c*sqrt(-2*E*m - _c[1])*t/h)]. I didn't find the value of _c[1] for the 2 additional condition above.

The issue is the period of time between 2 consecutive zero of the f(r,t)=0

introducing the iv1 f(rb,0) didn't work The compiling time over 15000sec and still not complete!

```infolevel[pdsolve] := 5: # optional
eq1 := 2*m*(E + 8*Pi*epsilon/r)*f(r, t)/h^2 + R*diff(f(r, t), r \$ 2)/r - diff(f(r, t), t \$ 2)/(a^2*c^2) = 0;
iv1 := f(r, 0) = 0, f(R, t) = 0, D[1](f)(0, 0) = R, f(rb,0)=0;
sys := [eq1, iv1]:
Sol := pdsolve(sys, HINT=`*`);
pt := pdetest(Sol, sys);```

Where f(r,t) is the function of variable r and t  in spherical coordinate and m, E, h, R, rb, a, and c are constants.

I also want to find the exact value of f(r,t) with the condition f(rb, 0) = 0; and diff(f(rb,t),t)=a*c for the value t=0 and if is possible the pulsation of the sinusoidal solution of f(r,t). [the solution is a combination of AiryAi ; AiryBi and sinusoidal sin(a*c*sqrt(-2*E*m - _c[1])*t/h)]. I didn't find the value of _c[1] for the 2 additional condition above.

The issue is the period of time between 2 consecutive zero of the f(r,t)=0

The solution must not contain the _c[1] it must result from the 2 conditions.

## Excelent answer but there is a problem!...

Excellent answer but something is not ok is entirely my fault the equation is rH''(r)+H'(r)+(rk^2-r^2*b^2/R^2)H(r)=0 where k, b, and R are real constant positive number, with condition H(R)=0 and H'(1/R)=R to be solved into series of power. The SECOND CONDITION IS TO DERIVATE H'(1/R)=R that's why I'm interested in an approximate solution based on series, or any results as long as it satisfied the too condition H(R)=0 and H(1/R)=R of the real function H(r) and if is possible a plot for R=370, k=100 and b=35. Very nice I learn a lot about it! could you add an extra help, please!

## sorry I missed H(r)...

Sorry, the equation is rH''(r)+H'(r)+(rk^2-r^2*b^2/R^2)H(r)=0 where k, b, and R are real constant positive number, with condition H(R)=0 and H(1/R)=R to be solved into series of power. I know from the literature that xy''+y'+xy=0, can't be solved in terms of elementary function(see G.Nagy-ODE-November 29, 2017) that's why I'm interested in an approximate solution based on series, or any results as long as it satisfied the too condition H(R)=0 and H(1/R)=R of the real function H(r) and if is possible a plot for R=370, k=100 and b=35.

## ok asympt is a good command!...

ok asympt is a good command!

asympt(G, r):

but I'm also interested in a general solution where R= an arbitrary constant, not a value!

## Not quite so!...

I'm very interested in an approximate solution but unfortunately is not correct you put R=0 and convert  it into series

`g1:=convert(series(%, R=0), polynom);  `

I have my own workout

 >
 >
 >
 (1)
 >
 (2)
 >
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)
 >
 (8)
 >
 (9)
 >
 (10)
 >
 >
 (11)
 >
 >

but I didn't resolve in a recurrent mod with the expression of the general term.