DSolve; Lorentz Force equations of motion

Click For Summary
SUMMARY

This discussion focuses on using DSolve in Mathematica 6.0.1.0 to derive equations of motion for an electron in crossed electric and magnetic fields. The user successfully derived the second-order differential equations governing motion but encountered difficulties in obtaining the correct solutions for x[t] and y[t] using DSolve. Key issues included handling constants of integration and generating parametric plots from the equations of motion. Restarting the kernel resolved some errors, but challenges with variable substitution and plotting persisted.

PREREQUISITES
  • Understanding of second-order differential equations
  • Familiarity with the Lorentz Force equation
  • Basic knowledge of Mathematica syntax and functions
  • Experience with parametric plotting in Mathematica
NEXT STEPS
  • Learn how to effectively use DSolve in Mathematica for differential equations
  • Explore the use of constants of integration in Mathematica
  • Research techniques for generating parametric plots in Mathematica
  • Investigate the GeneratedParameters function in Mathematica for unique variable representation
USEFUL FOR

Students and researchers in physics, particularly those studying electromagnetism and differential equations, as well as Mathematica users seeking to solve complex equations of motion.

KleZMeR
Messages
125
Reaction score
0
Hi,

I am trying to use DSolve in Mathematica 6.0.1.0, in order to find equations of motion (x[t], y[t]) from second-order differential equations. I have looked through much documentation, and attempted numerous codes to figure this out, but get many errors or the same answer every time. This is my first Mathematica course. I derived the second order equations by hand, which is found in many books.
I do not see how the equations of motion I am to show, x[t], y[t], are produced by DSolve.
Below is stated problem and attempt.

Electron injected with initial velocity Vo into crossed field; E=Ey, B=Bz,
and: w=omega, Vd=(Ey/Bz)
Also given: w=1, Vox=Voy=2, Vd=1, (I do not know why these are given, unless they are initial conditions, and even so, they still do not change my outcome)

a) Use F = q (E + V x B) and show by hand that motion is governed by these equations (I did this):

x''[t] + (w^2)x = (w^2)*Vd*t - w*Voy

y''[t] + (w^2)y = -(w^2)*Vd*t + w*Vox



b) Use DSolve to show that the solutions to these equations are:

x[t] = (1/w)*(Vox - Vd)*sin(wt) + (Voy/w)*(cos(wt) - 1) + Vd*t

y[t] = (Voy/w)*sin(wt) - (1/w)*(Vox - Vd)*(cos(wt) - 1)




So I am having trouble with part b)... I am sure that I am missing much in my attempt at solving it, but all other inputs either give me an error, or do not change the output, any help or suggestions would be appreciated, thanks.
 

Attachments

Physics news on Phys.org
DSolve[x''[t] + Subscript[\[Omega], c]^2 x[t] ==
Subscript[\[Omega], c]^2 Subscript[v, d]*t -
Subscript[\[Omega], c]*Subscript[v, 0], x[t], t]i did it and it works

Try restarting the kernel (Evaluation - Quit Kernal - Local (or other)) and just running that one line.
 
Awsome! Thank you so much K.J. !

I would also like to replace the two constants of integration(C[1], C[2]) with two chosen variables

If I include this in DSolve:

{C[1] -> Vox, C[2] -> Voy},

it returns that I can not use these for variables.

I have also found a GeneratedParameters-> function, but it does not work either,.. and if I'm correct, the GeneratedParameters-> Module{C[1], C[2]..&} function is only to ensure that C[] values are all unique, and does not change their representation.

And last, these two equations of motion, with initial conditions, don't seem to give a parametric plot, just an empty axis. I am attaching this attempt at the parametric plot. This form does plot a different set of equations, so perhaps it is the equation representation,

Again, any help would be appreciated
Thanks
 

Attachments

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 20 ·
Replies
20
Views
2K
Replies
3
Views
2K
  • · Replies 18 ·
Replies
18
Views
3K
  • · Replies 11 ·
Replies
11
Views
910
  • · Replies 10 ·
Replies
10
Views
3K
Replies
0
Views
1K