Mark McClure
Reputation
25,044
200/100 score
 Sep 9 comment Möbius transformations revealed @YvesKlett Seems like a reasonable expectation. :) I do have an implementation based on simple graphics primitives that I'll post soon. Sep 9 comment Calculating a potential function using the finite element method @s.s.o I didn't notice the {0,0} issue - thanks! Sep 9 comment Calculating a potential function using the finite element method @s.s.o The value of MaxBoundaryCellMeasure is ignored in your code, which is why it alleviates the problem. Using ToElementMesh directly on reg allows one to at least reduce the value of MaxBoundaryCellMeasure somewhat before the problem starts. Sep 9 comment Calculating a potential function using the finite element method @Hugh Essentially, that's correct. Technically, though, it's interpolating, rather than extrapolating. :) Sep 9 comment Calculating a potential function using the finite element method The return value of NDSolveValue is an InterpolatingFunction, which describes how to compute values between points on the grid. That answers your question 4 and also indicates what's going on in question 5 - namely the interpolation order is too low to expect to do better. Unfortunately, I don't think it's so easy to increase that on an unstructured grid. Also, I'm not so sure how well MaxBoundaryCellMeasure works and your problem with the mes on the circle seems to be alleviated when you delete it. That option isn't even available to DiscretizeRegion. Sep 8 comment Should eigenvalues be ordered? So, what do you think of the output of this: Sort[{1, 2, 4, Sqrt[Pi]}]? I think it's reasonable in a symbolic system. Sep 8 comment How can Mathematica help me to find a real radical expression for roots of this polynomial?‎ See Casus Irreducibilis and/or this notebook. The roots can be expressed without the imaginary unit, if you are willing to accept trig functions - just hit your output with ComplexExpand. Sep 3 comment How to speed up my Project Euler code If you realize that the maximum solution is 10 Floor[Ceiling[Sqrt[1929394959697989990]]/10] and step down by 10 from there, you'll get the solution almost immediately. Sep 2 comment How to Nest a Function with three Arguments? I suspect he'd like it to work with arguments other than A, B, and dt. Perhaps, a more general pattern would be appropriate? Sep 2 comment How to find this limit correctly? I suppose that Assumptions in Limit might simply be applied to simplify the expression ahead of the computation. That would explain the difference in the results and be dismissed as designed. A genuine domain restriction would be more properly dealt with as a third argument. Sep 2 comment How to find this limit correctly? @DanielLichtblau That is what I thought. But, then, how do I explain the fact that Limit[n^2 Sin[2*n*Pi], n->Infinity, Assumptions->Element[n, Integers]]==0? Sep 1 comment How to find this limit correctly? @belisarius Well, yeah, but that's a discrete plot indicating the limit is $2\pi$ ( I assume, I'm on my iPhone). Sep 1 comment Importing High Precision (22 digits) data from a file @jason There is no import as string in my example. I used ImportString since it works exactly as Import but allows one to show a whole example without reference to an internal file. If you want that specific number, then that specific number will need one more zero in your input file. Keep in mind also that you should get the binary number with the specified precision closest to your input, so that trailing digits are often present. I can't verify your 8.24574098106 as I get 8.24574097909*^8, which seems quite right. Sep 1 comment Use Mathematica to calculate the area enclosed between two curves Is there an obfuscated Mathematica contest?? :) Aug 31 comment Specifying Range of RSolve @QuinnCulver I mean that $|x_{n}-3/4|<(|x_{n-1}-(3/4)|)/2$. In words, the distance between $x_n$ and $3/4$ is a little less than half the distance between $x_{n-1}$ and $3/4$. You can verify this numerically by looking at Ratios[N[Table[q[n],{n,0,9}]-3/4]] - you'll see a sequence converging to $1/2$. Analytically, you're iterating $f(x)=(-2x+3\sqrt{4x+1}-3)/2$ which has a fixed point at $x=3/4$ satisfying $f'(3/4)=1/2$. Aug 28 comment Exporting as GIF No problem - I'm glad you figured it out. Aug 19 comment Coordinate extraction and manipulating graphs as graphics objects Oh, it's just one tiny little change. I thought about editing your answer but you can do it, if you want. And, yes, I've posted the GraphConvert thing before but more recently learned of the more intuitive approach with Show. Aug 19 comment Coordinate extraction and manipulating graphs as graphics objects The GraphComputationGraphConvertToGraphics function is called by Show when applied to a Graph. Thus, Show[g] accomplishes the same thing. Aug 11 comment DSolve problem with system of linear ODEs I have to admit that I find it an odd syntactical difference that NDSolve works while DSolve does not. I don't think I'd call it a bug, though. Aug 11 comment DSolve problem with system of linear ODEs You need to define X: X[t_] = {x[t], y[t], z[t]}`.