For questions regarding the underlying implementation of built-in functions.

learn more… | top users | synonyms

4
votes
1answer
402 views

How does Mathematica calculate the nth prime?

When I enter Prime[2000000000000], the two-trillionth prime, Mathematica gives 61427839512211 for the answer after several ...
4
votes
1answer
152 views

What are the hidden specifications for FindRoot

The Help page of FindRoot says: "by default, FindRoot uses Newton's method (Newton-Raphson) to solve a nonlinear system" (or a nonlinear equation I ...
4
votes
0answers
118 views

Given two samples, how does LocationTest select its AutomaticTest?

How does LocationTest select its "AutomaticTest" when testing two samples? Background: A previous post poses a similar question ...
4
votes
0answers
87 views

JewishNewYear date limitation

The documentation for the function "JewishNewYear" in the "calendar" package declares that the year argument must be between 1900 and 2100. No explanation is given for this limitation. The function ...
3
votes
2answers
222 views

What does MinimalPolynomial do?

suppose $x^5-x-1=0$, $y^7-y-1=0$ and $z=x+y$. I want to find a minimal polynomial expression of $z$, such that $p(z)=0$, and the code can be written like this: ...
3
votes
3answers
141 views

Pick bug? level 0 matching

A previous discussion of this issue (Is there a bug in Pick?) accepted that matching at level 0 is not a bug in Pick. Was that correct? I am not yet persuaded. ...
3
votes
2answers
208 views

How does RootOfUnityQ work?

How does Mathematica's RootOfUnityQ function work? That is, how does Mathematica know if a number is a root of unity? Example: Let $x = \frac{1-i \sqrt{2+\sqrt{5}}}{1+i \sqrt{2+\sqrt{5}}}$. Then ...
3
votes
1answer
201 views

What method does `NIntegrate` uses by default?

There is a variety of algorithm for performing numerical integration (See wiki). What method does NIntegrate uses by default? I looked on the documentation page and I saw that the function ...
3
votes
1answer
262 views

Daubechies scaling function

In Mathematica in wavelet package with the following codes Phi[x_, j_, k_] := 2^(j/2) WaveletPhi[DaubechiesWavelet[3], 2^(j) x - k] Mathematica computes ...
3
votes
1answer
144 views

Root finding: zeroes of Mathieu function

I am finding the roots of the Mathieu sine function, and find Mathematica and Maple do not agree on the solutions. For example, consider the solutions of ...
3
votes
0answers
63 views

Custom Methods for NMinimize/NonLinearModelFit

Is there a way to write custom methods for NonLinearModelFit or NMinimize? To give background, I currently have a batch of code that uses NonLinearModelFit with NMinimize as the method to do ...
3
votes
0answers
146 views

Distance Time Warp (DTW) function implementation in mathematica

The DTW (Distance time warp) is given here http://en.wikipedia.org/wiki/Dynamic_time_warping There is an implementation of DTW in R as well as MATLAB. However i am unable to find something for DTW in ...
3
votes
0answers
288 views

How to determine BLAS/LAPACK implementation used internally for numerical matrix operations?

Is there a command which reveals which implementation of BLAS and LAPACK are used in Mathematica's matrix operations such as Eigensystem? I asked a related question ...
2
votes
1answer
314 views

How does FindShortestTour work?

How are FindShortestTour and the algorithms therewithin implemented internally? Also: Can I access some form of incremental, in-progress result during evaluation? ...
2
votes
1answer
112 views

What is the formula Mathematica uses for ZetaZero?

What is the formula/algorithm Mathematica uses for the ZetaZero command?
2
votes
1answer
153 views

Implementation of Decompose

I'm curious as to how Decompose works so I decided to use Trace with the option ...
2
votes
2answers
111 views

Derivative of Vectoral NDSolve Interpolated Function / Rule Solutions,/., and Vector NDSolve [duplicate]

I need to find the instantaneous phase and instantaneous frequency for a rotating system. I am having difficulty making the math work with interpolating functions. First, here is a problem without ...
2
votes
1answer
102 views

Anyone knows the algorithm used by NullSpace function?

NullSpace function gives a list of vectors that forms a basis for the null space of the input matrix. When the rank of the input argument matrix $M_{m\times n}$ is ...
2
votes
1answer
339 views

Problem with Covariance Matrix Output in NonlinearModelFit

I am running NonlinearModelFit based off of some simulated data and trying to fit to a function with more than one parameter. Eventually, I would like to fit to 5 ...
2
votes
0answers
36 views

The idea behind Stiffness switching method with NDsolve [closed]

Does the Stiffness switching method with NDsolve switch just between multiple variants of 4th order Runge Kutta method or it uses also other methods?
2
votes
0answers
39 views

How does Mathematica automatically choose which classifier method to use?

I'm interested in what heuristics are used in Classify's choice of method.
2
votes
0answers
71 views

PlotRange->Automatic the exact function used to calculate outliers

I know that PlotRange->Automatic does the following: "the distribution of coordinate values is found, and any points sufficiently far out in the distribution ...
2
votes
0answers
51 views

Implementation of Series function

I'm trying to evaluate the partial derivatives of some function F[x,y] at some point, i.e. ...
2
votes
0answers
57 views

Can BreadthFirstScan be modified?

I would like to obtain a "slick" implementation of A*-search algorithm on Graph objects, using the built-in BreadthFirstScan (if that is possible at all). I noticed that BreadthFirstScan is static in ...
2
votes
0answers
82 views

What are the limits of the Prime-functions? [duplicate]

Will Prime and PrimePi test up to 3 * 10^12? Where can I find the limits for these funtions? ...
2
votes
0answers
140 views

PowersRepresentations Algorithm

I'm trying to understand the mathematics behind counting the number of representations of a positive integer by $n$ distinct $k$th powers, i.e. I would really like to know how to do the Mathematica ...
1
vote
2answers
93 views

Does FindFit use symbolic differentiation?

Does FindFit use symbolic differentiation when optimizing by least squares a non-linear function?
1
vote
2answers
273 views

How to generate a unique file name similar as that created by CreateTemporary

CreateTemporary can generate a unique file name like ...
1
vote
1answer
137 views

Open an MX file used by the kernel [closed]

I am working on a graph, and I am using Mathematica function FindGraphCommunities to analyze the properties of this graph, but I have no idea about the algorithm ...
1
vote
1answer
442 views

Method used by FindFit

Using FindFit, Mathematica selects an optimal algorithm when Automatic is selected. I'd like to know the name of the method ...
1
vote
0answers
28 views

Apparent issue with derivative of FractionalPart

This issue came to my attention from Math.SE: http://math.stackexchange.com/questions/1015325/is-the-derivative-of-x-on-0-1-always-equal-to-1/1015342#1015342 To summarize, it appears that ...
1
vote
0answers
137 views

How to write my own interpolation method

I was comparing results of spline interpolation in Mathematica and C++/gsl-library Results of interpolation are different if I use ...
1
vote
0answers
36 views

Internals of Prime function [duplicate]

Does anyone knows how exactly Prime works ? For example how Prime[1000000000] is calculated ? The only information I found was ...
1
vote
0answers
973 views

How to Output Chi-Squared Statistics when using NonLinearModelFit

I am using NonLinearModelFit for some curve fitting and I was wondering if NLM is able to output chi-squared/leastsquared statistics from the best-fit parameters and confidence intervals. From my ...
1
vote
0answers
83 views

What algorithm is Mathematica using for FirstPassageTimeDistribution[]?

What algorithm is Mathematica 9 using for FirstPassageTimeDistribution[]? Is it attempting to use PseudoInverse[] on the ...
0
votes
1answer
71 views

Mathematica Lattice Reduce Command

I'm going through a very old copy of "Mathematica: A System for Doing Mathematics by Computer" for self practice. I'm on chapter 3, and ran into the LatticeReduce ...
0
votes
0answers
29 views

Algorithm for converting truth tables to Boolean functions

The function BooleanFunction is able to convert a truth table into a Boolean function. Does anyone know what is the mathematical background and/or algorithm that ...
0
votes
0answers
51 views

How to make work Plus function? [duplicate]

I have a little problem. I have k = 1.56343*10^9; r = 9.10492*10^18. And I should receive a (2*r*r-k*k)/2*r*r, it equals to ...
0
votes
0answers
60 views

Plot function failure, and understanding how mathematica selects recursion points for correction

I think it is best to start out by directing everyone to this link. The discussion concerns why the Plot function fails to accurately generate a simple sine wave for given parameters. The answer given ...
0
votes
0answers
61 views

Birch clustering in Mathematica

I want to implement birch clustering in Mathematica. MAthematica only has 4 methods of clustering, is there any implementation for Birch clustering? Birch clustering is much faster than K means ...
0
votes
0answers
51 views

What Algorithm Does Mathematica use to Calculate Betweenness Centrality for Unweighted Graphs?

Recently asked this question http://math.stackexchange.com/questions/888369/betweenness-centrality-how-long-does-mathematica-take on the Mathematics SE. Does anyone know, as a follow up, what ...
0
votes
0answers
265 views

What algorithms does NSolve use?

I've been using NSolve to find roots and now I've been wondering what underlying algorithm mathematica uses. It's probabaly similiar to ...
0
votes
0answers
126 views

How does Mathematica solve eigenvalues?

I believe mathematica uses some techniques beyond our book. According to the text book it is not possible to solve eigenvalues of which the matrix has not completely eigenvectors, or of those matrix ...
-2
votes
1answer
110 views

Mathematica's steps towards the result of a definite integral

I am interested in this definite integral: $$ \eqalign{\int_0^\infty {x^2 \log(1-e^{-x})}\; dx\cr } $$ That is: ...
-3
votes
1answer
65 views

Miller-Rabin algorithm [closed]

I want to implement the Miller-Rabin algorithm in Mathematica to check if a number is prime with at least 99.99% probability. I used this: ...