Skip to main content

Unanswered Questions

646 questions with no upvoted or accepted answers
10 votes
0 answers
1k views

Modelling Hysteresis with a Differential Equation

I want to implement the bulk ferromagnetic hysteresis model (mostly the Jiles-Atherton Model), see http://drum.lib.umd.edu/bitstream/1903/6043/1/PhD_99-1.pdf page 44 equation (30). The needed ...
9 votes
0 answers
977 views

Solving the 2D Schrödinger equation with eigensystem, then verifying orthonormality of eigenfunctions with NIntegrate

I am solving the time-independent 2D Schrödinger equation for an interacting electron and hole in the case of anisotropic electron and hole masses, where the interaction is described a modified form ...
9 votes
0 answers
562 views

Spurious Error Messages from NDSolve when Using WhenEvent with Time Delays

Bug introduced between versions 10.1 and 10.4, and resolved in 11.3. Using either 11.2 or 10.4 on Windows 10 (64 bit), I am unable to reproduce the answer by March to question 99576. Specifically, ...
8 votes
0 answers
197 views

Bug in integral related to beta distribution

I've encountered a problem, which I'm pretty sure is a bug. I've contacted Wolfram's support and they were less than helpful. I'd like to bring the issue up here to get either a confirmation that it ...
7 votes
0 answers
436 views

Integrating rational functions of several variables over $\mathbb{H}^4$

Let $W$ be a rational function of $8$ variables $a,b,c,d,e,f,g,h$ from this file, e.g.: ...
6 votes
0 answers
186 views

ParametricNDSolveValue causes kernel to crash

Bug introduced in 12.1.1 and persisting through 13.2.0. In the course of answering question 228693, I found that ...
6 votes
0 answers
7k views

What's the most difficult multidimensional integral that Mathematica has solved?

I am interested in benchmarking numerical integration methods and am trying to develop a wider set of difficult multivariate examples. For my particular methods, I only want to look at non-negative-...
6 votes
0 answers
211 views

Details of NDSolve calling LSODA

Inspired by my question regarding the computation time of NDSolve using the LSODA backend I was wondering how NDSolve is actually calling LSODA (what arguments are sent to LSODA), i.e. what are the ...
6 votes
0 answers
511 views

What are good/best practices to take the Fourier transform of an InterpolatingFunction?

I have a function which I have obtained from numerical integration of a differential equation, and I would like to take its Fourier transform. What are good practices for doing this? To make things ...
6 votes
0 answers
187 views

Modify NDSolve`StateData (if possible)

I am trying to solve a PDE that needs to be scaled constantly (refer to this). @andre suggests I modify NDSolve`StateData. Now, the problem is, I'm not used to ...
5 votes
0 answers
147 views

"Glueing" interpolating functions

Suppose we have two interpolating functions with exactly touching domains, in our case {-1, 1} and {1, 2}. How to "glue&...
5 votes
0 answers
249 views

Hypergeometric Function Integration Using Mellin-Barnes Representation

I have the following integral: $$ I=\int_0^1 d \alpha d \beta d \gamma r(\gamma) s(\alpha, \beta) T(\alpha, \beta, \gamma) $$ where I define $$r(\gamma)=(\gamma(1-\gamma))^{ -1 / 2+\epsilon / 2}$$ and,...
5 votes
0 answers
224 views

Are there any dedicated built-in methods to NIntegrate over the unit sphere or the 3D rotational group?

As part of a molecular-physics calculation, I need to perform an integral to find the average of a certain function over all the possible orientations of the molecule, $$ \langle f(\mathbf v)\rangle = ...
5 votes
0 answers
182 views

Numerical verification of the estimate:

How to verify numerically with considerable accuracy in Mathematica the following : $$\int_2^x\dfrac{1}{z\Gamma(\sin^2[π\Gamma(z)/(2z)])}dz\sim\ln(\ln(x))$$ ? I need more suitable and better code ...
5 votes
0 answers
276 views

StiffSystem or Singularity - a system of second order ODEs in the problem of geodesics

I would be extremely grateful for any help regarding the following code I wrote and the errors it produces. In this code I am investigating the behaviour of a massive particle trapped in the vicinity ...

15 30 50 per page
1
2 3 4 5
44