It has been long known that a year after Schrödinger published his equation, Madelung also published a hydrodynamics version of Schrödinger equation. Quantum diffusion is studied via dissipative Madelung hydrodynamics. Initially the wave packet spreads ballistically, than passes for an instant through normal diffusion and later tends asymptotically to a sub‐diffusive law. In this paper we will review two different approaches, including Madelung hydrodynamics and also Bohm potential. Madelung formulation leads to diffusion interpretation, which after a generalization yields to Ermakov (...) equation. Since Ermakov equation cannot be solved analytically, then we try to find out its solution with Mathematica package. It is our hope that these methods can be verified and compared with experimental data. But we admit that more researches are needed to fill all the missing details. (shrink)

Following a hypothesis by Marciak-Kozlowska, 2011, we consider one-dimensional Schumann wave transfer phenomena. Numerical solution of that equation was obtained by the help of Mathematica.

In a recent paper in Astrophysics and Space Science Vol. 364 no. 11 (2019), S. Ershkov & D. Leschenko presented a new solving procedure for Euler-Poisson equations for solving momentum equations of the CR3BP near libration points for uniformly rotating planets having inclined orbits in the solar system with respect to the orbit of the Earth. The system of equations of the CR3BP has been explored with regard to the existence of an analytic way of presentation of (...) the approximated solution in the vicinity of libration points. A new and elegant ansatz has been suggested in their publication whereby, in solving, the momentum equation is reduced to a system two coupled Riccati ODEs. In this paper, we presented a numerical solution of such coupled Riccati ODEs using Mathematica software package. (shrink)

A new class of polynomials investigates the numerical solution of the fractional pantograph delay ordinary differential equations. These polynomials are equipped with an auxiliary unknown parameter a, which is obtained using the collocation and least-squares methods. In this study, the numerical solution of the fractional pantograph delay differential equation is displayed in the truncated series form. The upper bound of the solution as well as the error analysis and the rate of convergence theorem are also investigated in (...) this study. In five examples, the numerical results of the present method have been compared with other methods. For the first time, a -polynomials are used in this study to numerically solve delay equations, and accurate approximations have been displayed. (shrink)

In this paper, we investigate the numerical solution of the coupled fractional massive Thirring equation with the aid of He’s fractional complex transform. This study plays a significant aspect in the field of quantum physics, weakly nonlinear thrilling waves, and nonlinear optics. The main advantage of FCT is that it converts the fractional differential equation into its traditional parts and is also capable to handle the fractional order, whereas the homotopy perturbation method is employed to tackle the nonlinear terms (...) in the coupled fractional massive Thirring equation. An example is illustrated to present the efficiency and validity of the two-scale theory. The solutions are obtained in the form of series with simple and easy computations which confirm that the present approach is good in agreement and is easy to implement for such type of complex systems in science and engineering. (shrink)

An important task in mathematical sciences is to make quantitative predictions, which is often done via the solution of differential equations. In this paper, we investigate why, to perform this task, scientists sometimes choose to use numerical methods instead of analytical solutions. Via several examples, we argue that the choice for numerical methods can be explained by the fact that, while making quantitative predictions seems at first glance to be facilitated by analytical solutions, this is (...) actually often much easier with numerical methods. Thus we challenge the widely presumed superiority of analytical solutions over numerical methods. (shrink)

Volterra's (1926) equations for competition and predator-prey interactions are modified by introduction of root terms. A critical comparison with the original equations shows that the dynamic properties of the systems remain essentially alike, while the modification allows for explicit solution of the differential equations. Detailed solutions and numerical examples are given.

In this paper, some novel analytical and numerical techniques are introduced for solving and analyzing nonlinear second-order ordinary differential equations that are associated to some strongly nonlinear oscillators such as a quadratically damped pendulum equation. Two different analytical approximations are obtained: for the first approximation, the ansatz method with the help of Chebyshev approximate polynomial is employed to derive an approximation in the form of trigonometric functions. For the second analytical approximation, a novel hybrid homotopy with Krylov–Bogoliubov–Mitropolsky method (...) is introduced for the first time for analyzing the evolution equation. For the numerical approximation, both the finite difference method and Galerkin method are presented for analyzing the strong nonlinear quadratically damped pendulum equation that arises in real life, such as nonlinear phenomena in plasma physics, engineering, and so on. Several examples are discussed and compared to the Runge–Kutta numerical approximation to investigate and examine the accuracy of the obtained approximations. Moreover, the accuracy of all obtained approximations is checked by estimating the maximum residual and distance errors. (shrink)

In this paper, we present an intensive investigation of the finite volume method compared to the finite difference methods. In order to show the main difference in the way of approaching the solution, we take the Burgers equation and the Buckley–Leverett equation as examples to simulate the previously mentioned methods. On the one hand, we simulate the results of the finite difference methods using the schemes of Lax–Friedrichs and Lax–Wendroff. On the other hand, we apply Godunov’s scheme to simulate the (...) results of the finite volume method. Moreover, we show how starting with a variational formulation of the problem, the finite element technique provides piecewise formulations of functions defined by a collection of grid data points, while the finite difference technique begins with a differential formulation of the problem and continues to discretize the derivatives. Finally, some graphical and numerical comparisons are provided to illustrate and corroborate the differences between these two main methods. (shrink)

This paper provides equational semantics for Dung's argumentation networks. The network nodes get numerical values in [0,1], and are supposed to satisfy certain equations. The solutions to these equations correspond to the ?extensions? of the network. This approach is very general and includes the Caminada labelling as a special case, as well as many other so-called network extensions, support systems, higher level attacks, Boolean networks, dependence on time, and much more. The equational approach has its conceptual (...) roots in the nineteenth century following the algebraic equational approach to logic by George Boole, Louis Couturat, and Ernst Schroeder. (shrink)

To solve fractional delay differential equation systems, the Laguerre Wavelets Method is presented and coupled with the steps method in this article. Caputo fractional derivative is used in the proposed technique. The results show that the current procedure is accurate and reliable. Different nonlinear systems have been solved, and the results have been compared to the exact solution and different methods. Furthermore, it is clear from the figures that the LWM error converges quickly when compared to other approaches. When compared (...) with the exact solution to other approaches, it is clear that LWM is more accurate and gets closer to the exact solution faster. Moreover, on the basis of the novelty and scientific importance, the present method can be extended to solve other nonlinear fractional-order delay differential equations. (shrink)

In this paper, DSEK model with fractional derivatives of the Atangana-Baleanu Caputo is proposed. This paper gives a brief overview of the ABC fractional derivative and its attributes. Fixed point theory has been used to establish the uniqueness and existence of solutions for the fractional DSEK model. According to this theory, we will define two operators based on Lipschitzian and prove that they are contraction mapping and relatively compact. Ulam-Hyers stability theorem is implemented to prove the fractional DSEK model’s (...) stability in Banach space. Also, fractional Euler’s numerical method is derived for initial value problems with ABC fractional derivative and implemented on fractional DSEK model. The symmetric properties contribute to determining the appropriate method for finding the correct solution to fractional differential equations. The numerical solutions generated using fractional Euler’s method have been plotted for different values of α where α ∈ 0,1 and different step sizes h. Result discussion will be given, describing the changes that occur due to the step size h. (shrink)

In this research, the He-Laplace algorithm is extended to generalized third order, time-fractional, Korteweg-de Vries models. In this algorithm, the Laplace transform is hybrid with homotopy perturbation and extended to highly nonlinear fractional KdVs, including potential and Burgers KdV models. Time-fractional derivatives are taken in Caputo sense throughout the manuscript. Convergence and error estimation are confirmed theoretically as well as numerically for the current model. Numerical convergence and error analysis is also performed by computing residual errors in the entire (...) fractional domain. Graphical illustrations show the effect of fractional parameter on the solution as 2D and 3D plots. Analysis reveals that the He-Laplace algorithm is an efficient approach for time-fractional models and can be used for other families of equations. (shrink)

We study the different solution strategy for solving epidemic model in different imprecise environment, that is, a Susceptible-Infected-Susceptible model in imprecise environment. The imprecise parameter is also taken as fuzzy and interval environment. Three different solution procedures for solving governing fuzzy differential equation, that is, fuzzy differential inclusion method, extension principle method, and fuzzy derivative approaches, are considered. The interval differential equation is also solved. The numerical results are discussed for all approaches in different imprecise environment.

An effective numerical scheme based on Euler wavelets is proposed for numerically solving a class of neutral delay differential equations. The technique explores the numerical solution via Euler wavelet truncated series generated by a set of functions and matrix inversion of some collocation points. Based on the operational matrix, the neutral delay differential equations are reduced to a system of algebraic equations, which is solved through a numerical algorithm. The effectiveness and efficiency of the (...) technique have been illustrated by several examples of neutral delay differential equations. The main advantages and key role of using the Euler wavelets in this work lie in the performance, accuracy, and computational cost of the proposed technique. (shrink)

The one-dimensional case of the homogeneous Hamilton–Jacobi and Bernoulli equations St $${\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 2}}\right.\kern-0em}\!\lower0.7ex\hbox{$2$}}$$ S x 2 =0, where S(x, t) is Hamilton's principal function of a free particle and also Bernoulli's momentum potential of a perfect liquid, is considered. Non-elementary solutions are looked for in terms of odd power series in t with x-dependent coefficients and even power series in x with t-dependent coefficients. In both cases, and depending upon initial conditions, unexpected regularities are observed (...) in the terms of these expansions and this suggests that S(x, t) should be written as a product of the elementary solution x2/(2t) and a function f=f(φ) where φ=φ(x, |t|) owing to the symmetry property which is that S(x, −t)=−S(x, t). Requiring that this Ansatz satisfies the said equation and choosing the simplest realization of φ(x, |t|)=φ0 |t/t0|κ (x/x 0)λ≥0 with κ, λ∈ ℝ results in a soluble ordinary differential equation, of first order in u=ln φ and quadratic in f. This ODE has two fixed points: f=1, obviously, and f=0, a new fixed point which is often called “trivial.” The phase plane (fu, f) consists of a family of parabolas, all of which pass through the two fixed points. Explicit solutions of the general case are given close to these fixed points. A one-parameter family of solution is found to emerge from the “trivial” fixed point with non-trivial initial values S(x, 0). Detailed analyses of these findings will be reported elsewhere, bearing in mind that Bernoulli's equation has to be supplemented by the continuity equation satisfied by the density of the liquid. (shrink)

Cancer chemotherapy has been the most common cancer treatment. However, it has side effects that kill both tumor cells and immune cells, which can ravage the patient’s immune system. Chemotherapy should be administered depending on the patient’s immunity as well as the level of cancer cells. Thus, we need to design an efficient treatment protocol. In this work, we study a feedback control problem of tumor-immune system to design an optimal chemotherapy strategy. For this, we first propose a mathematical model (...) of tumor-immune interactions and conduct stability analysis of two equilibria. Next, the feedback control is found by solving the Hamilton–Jacobi–Bellman equation. Here, we use an upwind finite-difference method for a numerical approximate solution of the HJB equation. Numerical simulations show that the feedback control can help determine the treatment protocol of chemotherapy for tumor and immune cells depending on the side effects. (shrink)

It is shown that the admissible solutions of the continuity and Bernoulli or Burgers' equations of a perfect one-dimensional liquid are conditioned by a relation established in 1949–1950 by Pauli, Morette, and Van Hove, apparently, overlooked so far, which, in our case, stipulates that the mass density is proportional to the second derivative of the velocity potential. Positivity of the density implies convexity of the potential, i.e., smooth solutions, no shock. Non-elementary and symmetric solutions of the (...) above equations are given in analytical and numerical form. Analytically, these solutions are derived from the original Ansatz proposed in Ref. 1 and from the ensuing operations which show that they represent a particular case of the general implicit solutions of Burgers' equation. Numerically and with the help of an ad hoc computer program, these solutions are simulated for a variety of initial conditions called “compatible” if they satisfy the Morette–Van Hove formula and “anti-compatible” if the sign of the initial velocity field is reversed. In the latter case, singular behaviour is observed. Part of the theoretical development presented here is rephrased in the context of the Hopf–Lax formula whose domain of applicability for the solution of the Cauchy problem of the homogeneous Hamilton–Jacobi equation has recently been enlarged. (shrink)

This paper discusses the stability and Hopf bifurcation analysis of the diffusive Kaldor–Kalecki model with a delay included in both gross product and capital stock functions. The reaction-diffusion domain is considered, and the Galerkin analytical method is used to derive the system of ordinary differential equations. The methodology used to determine the Hopf bifurcation points is discussed in detail. Furthermore, full diagrams of the Hopf bifurcation regions considered in the stability analysis are shown, and some numerical simulations of (...) the limit cycle are used to confirm the theoretical outcomes. The delay investment parameter and diffusion coefficient can have great impacts on the Hopf bifurcations and stability of the business cycle model. The investment parameters for the gross product and capital stock as well as the adjustment coefficient of the production market are also studied. These parameters can cause instability in, and the stabilization of, the business cycle model. In addition, we point out that, as the delay investment parameter increases, the Hopf bifurcation points for the diffusion coefficient values decrease considerably. When the delay investment parameter has a very small value, the solution of the business cycle model tends to become steady. (shrink)

Plant protection problems are simulated by a system of ordinary differential equations with given initial conditions. The sensitivity and resistance of pathogen subpopulations to fungicide mixtures, fungicide weathering, plant growth, etc. are taken into consideration. The system of equations is solved numerically for each set of initial conditions and parameters of the disease and fungicide applications. Optimization algorithms were investigated and a computer program was developed for optimization of these solutions. 14 typical cases of the disease were (...) simulated and optimized in order to determine optimal fungicide treatments. The optimized strategy for fungicide application differs considerably from the commonly used method and seems to be an important new principle in plant protection. The approach developed in this study may be useful for a wide spectrum of purposes in the simulation of leaf diseases. It may also help the biologist to decrease or pinpoint experimental work and analyze its results and is perspective for plant disease control. (shrink)

This paper proposes the modified generalization of the HSS to solve a large and sparse continuous Sylvester equation, improving the efficiency and robustness. The analysis shows that the MGHSS converges to the unique solution of AX + XB = C unconditionally. We also propose an inexact variant of the MGHSS and prove its convergence under certain conditions. Numerical experiments verify the efficiency of the proposed methods.

We use the Tolman–Oppenheimer–Volkoff equations for a Chaplygin type fluid to study, analytically and numerically, the global behavior of static solutions of spherically symmetric objects. Two possible regimes are especially investigated. The first one is the phantom regime in which the pressure module exceeds the energy density. In this case the equator is absent and all the solutions have the geometry of a truncated spheroid with the same kind of singularity. The second case is the normal regime (...) for which we determine all the solutions, excluding the de Sitter one, corresponding to a tri-dimensional spheroidal geometry. Beyond the equator, three possible cases are considered; the first case has a closed spheroid characterized by a Schwarzschild-kind singularity with an infinite blue-shift at the south pole, the second case configuration has a regular spheroid and the third case has configuration of a truncated spheroid having a scalar curvature singularity to a finite value of the radial distance. We also compare all the geometric configurations with ones obtained in the special case of Chaplygin gas. (shrink)

In this paper, we consider the iterative algorithm for a boundary value problem of n -order fractional differential equation with mixed integral and multipoint boundary conditions. Using an iterative technique, we derive an existence result of the uniqueness of the positive solution, then construct the iterative scheme to approximate the positive solution of the equation, and further establish some numerical results on the estimation of the convergence rate and the approximation error.

Filtration combustion is described by Laplacian growth without surface tension. These equations have elegant analytical solutions that replace the complex integro-differential motion equations by simple differential equations of pole motion in a complex plane. The main problem with such a solution is the existence of finite time singularities. To prevent such singularities, nonzero surface tension is usually used. However, nonzero surface tension does not exist in filtration combustion, and this destroys the analytical solutions. However, a (...) more elegant approach exists for solving the problem. First, we can introduce a small amount of pole noise to the system. Second, for regularisation of the problem, we throw out all new poles that can produce a finite time singularity. It can be strictly proved that the asymptotic solution for such a system is a single finger. Moreover, the qualitative consideration demonstrates that a finger with 1 2 of the channel width is statistically stable. Therefore, all properties of such a solution are exactly the same as those of the solution with nonzero surface tension under numerical noise. The solution of the ST problem without surface tension is similar to the solution for the equation of cellular flames in the case of the combustion of gas mixtures. (shrink)

The role of mathematics in scientific practice is too readily relegated to that of formulating equations that model or describe what is being investigated, and then finding solutions to those equations. I survey the role of mathematics in: 1. Exact solutions of differential equations, especially conformal mapping; and 2. Simulations of solutions to differential equations via numerical methods and via agent-based models; and 3. The use of experimental models to solve equations (...) (a) via physical analogies based on similarity of the form of the equations, such as Prandtl’s soap-film method, and (b) the method of physically similar systems. Two major themes emerge: First, the role of mathematics in science is not well described by deduction from axioms, although it generally involves deductive reasoning. Creative leaps, the integration of experimental or observational evidence, synthesis of ideas from different areas of mathematics, and insight regarding analogous forms are required to find solutions to equations. Second, methods that involve mappings or transformations are in use in disparate contexts, from the purely mathematical context of conformal mapping where it is mathematical objects that are mapped, to the use of concrete physical experimental models, where one concrete thing is shown to correspond to another. (shrink)

Much of scientific inference involves fitting numerical data with a curve, or functional relation. The received view is that the fittest curve is the curve which best balances the conflicting demands of simplicity and accuracy, where simplicity is measured by the number ofparameters in the curve. The problem with this view is that there is no commonly accepted justification for desiring simplicity. This paper presents a measure of the stability of equations. It is argued that the fittest curve (...) is the curve which best balances stability and accuracy. The received view is defended with a proof that simplicity corresponds to stability, for linear regression equations. 1This paper is based on part of my doctoral dissertation. My thanks go to my thesis supervisor Professor Alasdair Urquhart for his encouragement, constructive criticism, and for directing me to several relevant articles: to my advisor Professor Ian Hacking for reminding me to concentrate on results that might have some application in the real world; and to my friend Wendy Brandts for sharing her ideas on a closely related problem. My thanks also to an anonymous referee of The British Journal of the Philosophy of Science for several helpful comments, to my friends and family for unfailing support, and to the Social Sciences and Humanities Research Council (awards 452-86-5885 and 453-87-0513) and the University of Toronto for financial assistance. (shrink)

Vesely (experimental physics, U. of Vienna, Austria) provides the basic numerical and computational techniques, followed by an explanation of specific problems of computational physics. Appendices address properties of computing machines and an outline of the technique of Fast Fourier Transformation. The first edition, published by Plenum Press, Ne.

In this paper, we show that reproducibility is a severe problem that concerns simulation models. The reproducibility problem challenges the concept of numerical solution and hence the conception of what a simulation actually does. We provide an expanded picture of simulation that makes visible those steps of simulation modeling that are numerically relevant, but often escape notice in accounts of simulation. Examining these steps and analyzing a number of pertinent examples, we argue that numerical solutions are importantly (...) different from usual mathematical solutions. They are do not merely approximate the latter, but introduce new problems, including issues of artificiality, stability, and well-posedness. Consequently, simulation modelling can attain reproducibility only to a certain degree because it is working with numerical solutions. (shrink)

It has been known for long time that most of the existing cosmology models have singularity problem. Cosmological singularity has been a consequence of excessive symmetry of flow, such as “Hubble’s law”. More realistic one is suggested, based on Newtonian cosmology model but here we include the vertical-rotational effect of the whole Universe. We review a Riccati-type equation obtained by Nurgaliev, and solve the equation numerically with Mathematica. It is our hope that the new proposed method can be verified with (...) observation data. (shrink)

The existence of singularities alerts that one of the highest priorities of a centennial perspective on general relativity should be a careful re-thinking of the validity domain of Einstein’s field equations. We address the problem of constructing distinguishable extensions of the smooth spacetime manifold model, which can incorporate singularities, while retaining the form of the field equations. The sheaf-theoretic formulation of this problem is tantamount to extending the algebra sheaf of smooth functions to a distribution-like algebra sheaf in (...) which the former may be embedded, satisfying the pertinent cohomological conditions required for the coordinatization of all of the tensorial physical quantities, such that the form of the field equations is preserved. We present in detail the construction of these distribution-like algebra sheaves in terms of residue classes of sequences of smooth functions modulo the information of singular loci encoded in suitable ideals. Finally, we consider the application of these distribution-like solution sheaves in geometrodynamics by modeling topologically-circular boundaries of singular loci in three-dimensional space in terms of topological links. It turns out that the Borromean link represents higher order wormhole solutions. (shrink)

We aim to construct exact and explicit solutions to a generalized Bogoyavlensky-Konopelchenko equation through the Maple computer algebra system. The considered nonlinear equation is transformed into a Hirota bilinear form, and symbolic computations are made for solving both the nonlinear equation and the corresponding bilinear equation. A few classes of exact and explicit solutions are generated from different ansätze on solution forms, including traveling wave solutions, two-wave solutions, and polynomial solutions.

In this article, the concept of system of symbolic 2-plithogenic linear equations and its solutions are introduced and studied. The Cramer's rule was applied to solve the system of symbolic 2-plithogenic linear equations. Also, provided enough examples for each case to enhance understanding.

On the basis of the exact solutions to the non-stationary spherically symmetric Einstein and Maxwell equations for dust matter and radial electromagnetic field, a model of a wormhole with the pulsating in time inner world and two static throats has been developed. It has been shown that such a wormhole with an arbitrary radius of the Gaussian curvature can connect both two different asymptotically flat space-times and two regions of the selfsame space-time. The problem of the fulfilment of (...) the energy conditions in this wormhole has been investigated, as well as the problem of its traversability investigation has been set. (shrink)

In 1904, Norwegian physicist Vilhelm Bjerknes published what would become a landmark paper in the history of meteorology. In that paper, he proposed that daily weather forecasts could be made by calculating later states of the atmosphere from an earlier state using the laws of hydrodynamics and thermodynamics (Bjerknes 1904). He outlined a set of differential equations to be solved and advocated the development of graphical and numerical solution methods, since analytic solution was out of the question. Using (...) these theory-based equations to produce daily forecasts, however, turned out to be more difficult than anticipated. Graphical solution techniques had limited success, and a first attempt to use .. (shrink)

The progress of Mathematics during the nineteenth century was characterised both by an enormous acquisition of new knowledge and by the attempts to introduce rigour in reasoning patterns and mathematical writing. Cauchy’s presentation of Mathematical Analysis was not immediately accepted, and many writers, though aware of that new style, did not use it in their own mathematical production. This paper is devoted to an episode of this sort that took place in Spain during the first half of the century: It (...) deals with the presentation of a method for numerically solving algebraic equations by José Mariano Vallejo, a late Spanish follower of the Enlightenment ideas, politician, writer, and mathematician who published it in the fourth (1840) edition of his book Compendio de Matemáticas Puras y Mistas, claiming to have discovered it on his own. Vallejo’s main achievement was to write down the whole procedure in a very careful way taking into account the different types of roots, although he paid little attention to questions such as convergence checks and the fulfilment of the hypotheses of Rolle’s Theorem. For sure this lack of mathematical care prevented Vallejo to occupy a place among the forerunners of Computational Algebra. (shrink)

The statistical properties of a single quantum object and an ensemble of independent such objects are considered in detail for two-level systems. Computer simulations of dynamic zero-point quantum fluctuations for a single quantum object are reported and compared with analytic solutions for the ensemble case.

The nonlinear differential equation governing the periodic motion of the one-dimensional, undamped, and unforced cubic-quintic Duffing oscillator is solved exactly by obtaining the period and the solution. The period is given in terms of the complete elliptic integral of the first kind and the solution involves Jacobian elliptic functions. We solve the cubic-quintic Duffing equation under arbitrary initial conditions. Physical applications are provided. The solution to the mixed parity Duffing oscillator is also formally derived. We illustrate the obtained results with (...) concrete examples. We give high accurate trigonometric approximations to the Jacobian function cn. (shrink)

1. One theory or many? In 2004 a very interesting and readable article by Lenore Blum, entitled “Computing over the reals: Where Turing meets Newton,” appeared in the Notices of the American Mathematical Society. It explained a basic model of computation over the reals due to Blum, Michael Shub and Steve Smale (1989), subsequently exposited at length in their influential book, Complexity and Real Computation (1997), coauthored with Felipe Cucker. The ‘Turing’ in the title of Blum’s article refers of course (...) to Alan Turing, famous for his explication of the notion of mechanical computation on discrete data such as the integers. The ‘Newton’ there has to do to with the well known numerical method due to Isaac Newton for approximating the zeros of continuous functions under suitable conditions that is taken to be a paradigm of scientific computing. Finally, the meaning of “Turing meets Newton” in the title of Blum’s article has another, more particular aspect: in connection with the problem of controlling errors in the numerical solution of linear equations and inversion of matrices,Turing (1948) defined a notion of condition for the relation of changes in the output of a computation due to small changes in the input that is analogous to Newton’s definition of the derivative. The thrust of Blum’s 2004 article was that the BSS model of computation on the reals is the appropriate foundation for scientific computing in general. By way of response, two years later another very interesting and equally readable article appeared in the Notices, this time by Mark Braverman and Stephen Cook (2006) entitled “Computing over the reals: Foundations for scientific computing,” in which the authors argued that the requisite foundation is provided by a quite different “bit computation” model, that is in fact prima facie incompatible with the BSS model. The bit computation model goes back to ideas due to Stefan Banach and Stanislaw Mazur in the latter part of the 1930s, but the first publication was not made until Mazur (1963).. (shrink)

The aesthetic field equations do not resemble the wave equation, nor was the motivation behind them the wave equation. Nevertheless, we show that there exists a solution to the field equations that satisfies the wave equation. Integrability is also satisfied by this solution. Previously we showed that the Aesthetic Field Equations have particle solutions. Now we see that the equations also have sinusoidal solutions.

The extended system method of Nosé and Hoover for the control of temperature of a classical ensemble if applied to the de Broglie-Bohm-Vigier formulation of quantum mechanics. This allows for the simulation of the motion of a quantum particle at a constant preset temperature. A specific algorithm for numerical solution of the resulting equations of motion, based on the application of the methods of molecular dynamics simulation, is provided.

Despite the many successes of the relativistic quantum theory developed by Horwitz et al., certain difficulties persist in the associated covariant classical mechanics. In this paper, we explore these difficulties through an examination of the classical. Coulomb problem in the framework of off-shell electrodynamics. As the local gauge theory of a covariant quantum mechanics with evolution paratmeter τ, off-shell electrodynamics constitutes a dynamical theory of ppacetime events, interacting through five τ-dependent pre-Maxwell potentials. We present a straightforward solution of the classical (...) equations of motion, for a test event traversing the field induced by a “fixed” event (an event moving uniformly along the time axis at a fixed point in space). This solution is seen to be unsatisfactory, and reveals the essential difficulties in the formalism at the classical levels. We then offer a new model of the particle current—as a certain distribution of the event currents on the worldline—which eliminates these difficulties and permits comparison of classisical off-shell electrodynamics with the standard Maxwell theory. In this model, the “fixed” event induces a Yukawa-type potential, permitting a semiclassical identification of the pre-Maxwell time scale λ with the inverse mass of the intervening photon. Numerical solutions to the equations of motion are compared with the standard Maxwell solutions, and are seen to coincide when λ≳10−6 seconds, providing an initial estimate of this parameter. It is also demonstrated that the proposed model provides a natural interpretation for the photon mass cut-off required for the renormalizability of the off-shell quantum electrodynamics. (shrink)

We study the existence of weak solutions for the coupled system of functional integral equations of Urysohn-Stieltjes type in the reflexive Banach spaceE. As an application, the coupled system of Hammerstien-Stieltjes functional integral equations is also studied.

Es soll ein Beitrag zur epistemischen Charakterisierung von Computersimulationen als jenseits von Experiment und Theorie geleistet werden. Es wird argumentiert, dass die in der Simulationstechnik eingesetzten Verfahren nicht numerische Lösungen liefern, sondern deren Dynamik mittels generativer Mechanismen imitieren. Die Computersimulationen in der Klimatologie werden als systematisches wie historisches Fallbeispiel behandelt. Erst "Simulationsexperimente" gestatten es, mittels Modellen eine Dynamik zu imitieren, ohne deren Grundgleichungen zu "lösen". /// Computer simulations will be characterized in epistemic respect as a method between experiment and theory. (...) It will be argued that simulations do not provide numerical solutions, rather they use generative mechanisms to imitate a dynamics. Climate science will be considered as a case both systematically and historically. Only simulation experiments allow to build models that imitate a dynamics without solving the relevant equations. (shrink)