Programming for Lorenz's System - What to Use?

  • Thread starter Thread starter javiergra24
  • Start date Start date
  • Tags Tags
    Model
javiergra24
Messages
19
Reaction score
0
Hi guys

I have to write a program for the lorenz's system of equations. What program/programming language should I use? Fortran, Matlab, Python, Octave, R? I can´t use a preprogrammed module, I have to write the program by myself (using maybe the finite differences method or Runge Kutta). What do yoy recommend me?

Thanks
 
Physics news on Phys.org
I stay away from Fortran, Octave and R (although these would all work). Matlab, python or C would all work well. To decide which one is appropriate, we would need to know more details about what you are doing / trying to accomplish; or at least what kind of programming background you have.
 
Hi
Unfortunately I can't use C, I'll try it with python + numpy (in my opinion it's more advanced than C). I have enough time (2 months) to do it.

I'll post the exact problem later. But the basic tasks are the solutions of the lorenzt's equation for fixes values of sigma and rho. Main task is the numerical algorithm for solving the equation.

Thank you very much.
 
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top