Example 7​Universal Harmony
Book:
RisingoftheJamesWebbSpaceTelescopeanditsFundamentalBlindnes
:Page150​​IsaacNewtonbelievedlikemanyscientistsofhistimeina
Universal(Cosmic)HarmonyinPhysicsandinArts.
Likein
Music
:​​
NewtonandtheMysteryoftheMajorSixth
andin
Paintings
.ThisEquationhasbeenpresentedinthe
book
atPage151,EquationB-9​​​​
- 
1
2
c
∂
¯
E
×
¯
Η

∂t
+
ε
0
¯
E
∇. 
¯
E
−
ε
0
¯
E
×∇×
¯
E
+
μ
0
¯
Η
∇. 
¯
Η
−
μ
0
¯
Η
×∇×
¯
Η

= 0
(B-9)​​​​ExampleofNewton's
Universal(Cosmic)Harmony
withinanarbitrarybeamoflight,withanarbitrary'ElectricandMagneticFieldIntensityDivision'g[x,y]inthe(x,y)plane,propagatinginthez-directionwiththespeedoflightc=
1
ϵ0μ0
.
In[]:=
Needs["DifferentialEquations`NDSolveProblems`"]
In[]:=
Needs["DifferentialEquations`NDSolveUtilities`"]
In[]:=
Needs["DifferentialEquations`InterpolatingFunctionAnatomy`"]
In[]:=
ϵ0=.
In[]:=
μ0=.
In[]:=
x=.
In[]:=
y=.
In[]:=
z=.
In[]:=
t=.
In[]:=
Get["VectorAnalysis`"]
General
:VectorAnalysis` is now obsolete. The legacy version being loaded may conflict with current functionality. See the Compatibility Guide for updating information.
In[]:=
InverseFunctions->True
Out[]=
InverseFunctionsTrue
In[]:=
Needs["DifferentialEquations`NDSolveProblems`"]
In[]:=
Needs["DifferentialEquations`NDSolveUtilities`"]
In[]:=
Needs["DifferentialEquations`InterpolatingFunctionAnatomy`"]
In[]:=
SetCoordinates[Cartesian[x,y,z]]
Out[]=
Cartesian[x,y,z]
In[]:=
{Coordinates[Cartesian],CoordinateRanges[Cartesian]}
Out[]=
{{x,y,z},{-∞<x<∞,-∞<y<∞,-∞<z<∞}}
In[]:=
g[z,t]=.
Unset
:Assignment on g for g[z,t] not found.
Out[]=
$Failed
In[]:=
f[x,y]=.
In[]:=
f[x,y]=.
Unset
:Assignment on f for f[x,y] not found.
Out[]=
$Failed
In[]:=
K1=1
Out[]=
1
In[]:=
ev={f[x,y]*g[t-(K1/z+1)*z*Sqrt[ϵ0]*Sqrt[μ0]],0,0}
Out[]=
f[x,y]gt-1+
1
z
z
ϵ0
μ0
,0,0
In[]:=
mv=(1/Sqrt[μ0])*Sqrt[ϵ0]*​​{0,f[x,y]*g[t-(K1/z+1)*z*Sqrt[ϵ0]*Sqrt[μ0]],0}
Out[]=
0,
ϵ0
f[x,y]gt-1+
1
z
z
ϵ0
μ0

μ0
,0
In[]:=
Book:
RisingoftheJamesWebbSpaceTelescopeanditsFundamentalBlindnes
:Page15,Equation(19)​​ExampleofaLASER-BEAMwithaGaussianIntensitydivison
-K2
2
x
+
2
y

,combinedwithanarbitrarydivisiong[x,y]inthe(x,y)plane,propagatinginthez-directionwiththespeedoflightc=
1
ϵ0μ0
.​​TheinputfortheElectricFieldIntensity{E(x,y,z,)=ev}andtheMagneticFieldIntensity{H(x,y,z,)=mv}hasbeensubstitutedintheFieldEquationfortheElectromagneticField(BookEquation19,page16).​​​​
- 
1
2
c
∂
¯
E
×
¯
Η

∂t
+
ε
0
¯
E
∇. 
¯
E
−
ε
0
¯
E
×∇×
¯
E
+
μ
0
¯
Η
∇. 
¯
Η
−
μ
0
¯
Η
×∇×
¯
Η

= 0
Equation(19)​​​​Whichcanbewrittenas:​​​​term1=-
1
2
c
∂
¯
E
×
¯
Η

∂t
​​​​term2=
ε
0
¯
E
(∇. 
¯
E
)​​​​term3=−
ε
0
¯
E
×(∇×
¯
E
)​​​​term4=
μ
0
¯
Η
∇. 
¯
Η
​​​​term5=−
μ
0
¯
Η
×∇×
¯
Η
​​​​Equation(19)equals:​​​​term1+term2+term3+term4+term5=0
In[]:=
Div[ev]
Out[]=
gt-1+
1
z
z
ϵ0
μ0

(1,0)
f
[x,y]