Home › Applications › Matlab

It looks like you're new here. If you want to get involved, click one of these buttons!

- 140.6K All Categories
- 104.3K Programming Languages
- 6.4K Assembler Developer
- 1.8K Basic
- 39.7K C and C++
- 4.2K C#
- 7.9K Delphi and Kylix
- 4 Haskell
- 9.6K Java
- 4.1K Pascal
- 1.3K Perl
- 1.9K PHP
- 505 Python
- 48 Ruby
- 4.3K VB.NET
- 1.6K VBA
- 20.8K Visual Basic
- 2.6K Game programming
- 309 Console programming
- 88 DirectX Game dev
- 1 Minecraft
- 108 Newbie Game Programmers
- 2 Oculus Rift
- 8.9K Applications
- 1.8K Computer Graphics
- 726 Computer Hardware
- 3.4K Database & SQL
- 520 Electronics development
- 1.6K Matlab
- 627 Sound & Music
- 254 XML Development
- 3.3K Classifieds
- 189 Co-operative Projects
- 175 For sale
- 189 FreeLance Software City
- 1.9K Jobs Available
- 598 Jobs Wanted
- 201 Wanted
- 2.9K Microsoft .NET
- 1.7K ASP.NET
- 1.1K .NET General
- 2.9K Miscellaneous
- 3 Join the Team
- 2 User Profiles
- 353 Comments on this site
- 58 Computer Emulators
- 1.8K General programming
- 137 New programming languages
- 596 Off topic board
- 161 Mobile & Wireless
- 35 Android
- 124 Palm Pilot
- 334 Multimedia
- 151 Demo programming
- 183 MP3 programming
- 0 Bash scripts
- 15 Cloud Computing
- 52 FreeBSD
- 1.7K LINUX programming
- 365 MS-DOS
- 0 Shell scripting
- 319 Windows CE & Pocket PC
- 4.1K Windows programming
- 879 Software Development
- 401 Algorithms
- 67 Object Orientation
- 84 Project Management
- 88 Quality & Testing
- 232 Security
- 7.5K WEB-Development
- 1.8K Active Server Pages
- 61 AJAX
- 2 Bootstrap Themes
- 55 CGI Development
- 19 ColdFusion
- 222 Flash development
- 1.4K HTML & WEB-Design
- 1.4K Internet Development
- 2.2K JavaScript
- 33 JQuery
- 284 WEB Servers
- 108 WEB-Services / SOAP

Welcome to the new platform of Programmer's Heaven! We apologize for the inconvenience caused, if you visited us from a broken link of the previous version. The main reason to move to a new platform is to provide more effective and collaborative experience to you all. Please feel free to experience the new platform and use its exciting features. Contact us for any issue that you need to get clarified. We are more than happy to help you.

cHenriques
Posts: **1**Member

in Matlab

hi everyone, I'm making a research work to my master's degree and I have a problem I will appreciate your help.

I want simulate classically (not quantum mechanics) a electron orbiting a proton (H atom) immerse in a random distribution electromagnetic field. This theory is known as SED http://www.bu.edu/simulation/publications/dcole/PD...

For now we will forget the random EM. Since I'm trying first simulate the electron orbit without influence of EM.

As you now classically the electron feel the Coulomb force and the reaction term. The reaction term make the electron spin down and shock with the nucleus, the proton.

The non-relativistic Lorentz/Dirac equation to the movement is:

m*z'' = (-e^2)*z/|z|^3 + REAC - e*(EM)

The EM are the E and B random equations but for now we will forget this term.

REAC is the reaction term =>

2/3 * e^2 /c^3 * z''' (so third order differential eq)

but REAC can be approximated to

2/3 * e^2 /c^3 * d/dt( (-e^2)/m * z / z^3 ) ---> first order equation

the problem is in 2D x-y, so the vector z have two components, the z is the radius = sqrt(x^2 + y^2)

I can plot the electron orbit only with the Coulomb term:

((-e^2)*z/|z|^3)

but I cannot plot correctly with the REAC term. I already try use the ode45, the ode23 and the ode23t

I also tried use the REAC in the first form (third order ode) and in the second form

I also tried use polar coordinates, I'm lost I need help the following code is the my successful implementation of ode without the REAC

function mov = semZPF2(t,y)

r=(sqrt(y(1)^2+y(2)^2));

e = 1.602177e-19;

m=9.10938291*10^-31;

mov(1)=y(3);

mov(2)=y(4);

mov(3)=(-e^2*y(1)/r^3)/m;

mov(4)=(-e^2*y(2)/r^3)/m;

return;

------------------------------------------------------

yi=[0 5.3*10^-11]; % initial position

vi=[23.0583421525753 0]; % initial velocity

[t,y]= ode45(@(t,y) semZPF2(t,y,), [0 1E-12] , [yi vi] );

plot (y(i,1), y(i,2))

Please I need your help to implement the REAC term

m*z'' = (-e^2)*z/|z|^3 + REAC

thanks,

Carlos Henriques

I want simulate classically (not quantum mechanics) a electron orbiting a proton (H atom) immerse in a random distribution electromagnetic field. This theory is known as SED http://www.bu.edu/simulation/publications/dcole/PD...

For now we will forget the random EM. Since I'm trying first simulate the electron orbit without influence of EM.

As you now classically the electron feel the Coulomb force and the reaction term. The reaction term make the electron spin down and shock with the nucleus, the proton.

The non-relativistic Lorentz/Dirac equation to the movement is:

m*z'' = (-e^2)*z/|z|^3 + REAC - e*(EM)

The EM are the E and B random equations but for now we will forget this term.

REAC is the reaction term =>

2/3 * e^2 /c^3 * z''' (so third order differential eq)

but REAC can be approximated to

2/3 * e^2 /c^3 * d/dt( (-e^2)/m * z / z^3 ) ---> first order equation

the problem is in 2D x-y, so the vector z have two components, the z is the radius = sqrt(x^2 + y^2)

I can plot the electron orbit only with the Coulomb term:

((-e^2)*z/|z|^3)

but I cannot plot correctly with the REAC term. I already try use the ode45, the ode23 and the ode23t

I also tried use the REAC in the first form (third order ode) and in the second form

I also tried use polar coordinates, I'm lost I need help the following code is the my successful implementation of ode without the REAC

function mov = semZPF2(t,y)

r=(sqrt(y(1)^2+y(2)^2));

e = 1.602177e-19;

m=9.10938291*10^-31;

mov(1)=y(3);

mov(2)=y(4);

mov(3)=(-e^2*y(1)/r^3)/m;

mov(4)=(-e^2*y(2)/r^3)/m;

return;

------------------------------------------------------

yi=[0 5.3*10^-11]; % initial position

vi=[23.0583421525753 0]; % initial velocity

[t,y]= ode45(@(t,y) semZPF2(t,y,), [0 1E-12] , [yi vi] );

plot (y(i,1), y(i,2))

Please I need your help to implement the REAC term

m*z'' = (-e^2)*z/|z|^3 + REAC

thanks,

Carlos Henriques

About & Contact / Terms of use / Privacy statement / Publisher: Lars Hagelin

Programmers Heaven articles / Programmers Heaven files / Programmers Heaven uploaded content / Programmers Heaven C Sharp ebook / Operated by CommunityHeaven LLC

© 1997-2013 Programmersheaven.com - All rights reserved.