Tag for the usage of RSolve, ReccurenceTable, and implementation of finite difference method(FDM) in Mathematica

learn more… | top users | synonyms (2)

20
votes
1answer
796 views

Has this implementation of FDM touched the speed limit of Mathematica?

Still, I'll use the implementation of the 1D FDTD method (you can simply understand it as a kind of explicit finite difference scheme for the Maxwell's equation) as the example. Just for completeness, ...
15
votes
1answer
281 views

Mathematica 9 cannot solve a difference equation that Mathematica 8 can

Consider RSolve[{a[n + 1] == 2 a[n] + n (1 + 2^n), a[1] == 1}, a[n], n] Mathematica 8 can solve this quickly, but Mathematica 9, even after being allowed to run ...
10
votes
1answer
509 views

Is it possible to find a limit of a sequence given by its recurrence relation?

I need to calculate a limit of a sequence given by its recurrence relation. I tried the following: ...
8
votes
2answers
181 views

Evolution of a trait in an asexual population (performance tuning)

I just coded a simple simulation module that looks at the evolution of a continuous trait in a haploid asexually reproducing population under density dependent competition in discrete time (i.e. ...
7
votes
1answer
512 views

Poisson PDE over a irregular region with FDM

The problem modeled by the Poisson PDE is related to the torsion of prismatic beams and I use the Finite Differences Method (FDM). I've managed to solve the equation over a rectangle region with ...
7
votes
1answer
215 views

How to find a recurrence relation for a sequence?

I have a sequence given by an explicit formula for n-th term: ...
6
votes
6answers
404 views

Recursion question

The first two terms of a sequence are $a_1 = 1$ and $a_2 = \frac {1}{\sqrt3}$. For $n\ge1$, $$a_{n + 2} = \frac {a_n + a_{n + 1}}{1 - a_na_{n + 1}}.$$ What is $|a_{2009}|$? So I tried (thanks to ...
6
votes
3answers
502 views

How can I solve a forward-backward-looking recurrence system?

I would like to (numerically) solve a forward-backward looking recurrence equation in Mathematica. The system is something like: ...
6
votes
1answer
186 views

RSolve[{y[0] == y0, y[t + 1] == (1 - (t + 1)/100) y[t]}, y[t], t]

When I evaluate this code: RSolve[{y[0] == y0, y[t + 1] == (1 - (t + 1)/100) y[t]}, y[t], t] I get an RSolve::bvnul error. I don't understand why. Yet, if I ...
5
votes
2answers
143 views

RSolve gives a Context::ssle error

Bug introduced in 9.0 or earlier and fixed in 10.0 I want to solve the following recurrence equation RSolve[{(k + 2) c[k + 2] q^(k + 1) - c[k] == 0}, c[k], k] ...
5
votes
2answers
70 views

Solution to Simultaneous Arithmetic/Geometric Mean Recursion Relations

I'm trying to solve the simultaneous convergent sequences of geometric/arithmetic means where $a_{n+1}=\frac{1}{2}(a_n+b_n)$ and $b_{n+1}=\sqrt{a_nb_n}$ and initial values are $a_0=1+x$ and $b_0=1-x$. ...
5
votes
1answer
64 views

Two dimensional recurrence

My actual problem is a bit more complicated, but let us first assume that we have a recurrence relation of the following type $a_{n,l} = a_{n-1,l} + a_{n-1,l-1} + a_{n,l-1}$ The boundary is given, ...
4
votes
1answer
679 views

RSolve Lotka-Volterra: A Closed Form Solution to this Difference Equation?

I would like to solve a difference equation, a Lotka-Volterra difference equation. I tried a model from http://www.stolaf.edu/people/mckelvey/envision.dir/nonDE.lotka-volt.html Using their ...
4
votes
5answers
293 views

Can I generate a “piecewise” list from a list in a fast and elegant way?

Suppose we have a number list a: n = 1000000; a = RandomReal[{0, 1}, {n}]; We want to generate a list ...
4
votes
1answer
327 views

Mathematica 10 issue: working with Entities slows down recurrence equations [duplicate]

I was looking at the Mathematica code given at http://community.wolfram.com/groups/-/m/t/326240 where example code is given to simulate the spreading of a supposed Ebola outbreak in Nigeria: ...
4
votes
1answer
143 views

How can I implement this recurrence equation in Mathematica?

I asked a similar question yesterday but today I realized that I was implementing the wrong function. I have tried changing the implementation but I can't get it to work. Problem: There is a river ...
3
votes
3answers
150 views

Tangled up in sequences and recurrence relations

I took a trip through sequences and recurrence relations and I got tangled up since the documentation didn't help much. I'm sure you can. Given a sequence ...
3
votes
1answer
102 views

How can I find the coefficients of the next recursive symbolic function?

I have the following function $$ f(x+3k) = -mf(x+2k)+f(x) $$. If we get $x=x+k$, then $$ f(x+4k) = -mf(x+3k)+f(x+k) = m^2f(x+2k)+f(x+k)-mf(x) $$. Also we obtain the coefficients of the final ...
3
votes
2answers
172 views

Rook walk and RecurrenceTable in Mathematica

This appears in a combinatorics book: $a(m,n)=2 a(m,n-1)+2 a(m-1,n)-3 a(m-1,n-1)$ It is a recurrence equation for the number of rook walks from $(0,0)$ to $(m,n)$. The initial conditions are: $ ...
3
votes
1answer
76 views

Kernel crash when evaluating NDSolveValue

I had a code that was working fine and I tried to add an extra few lines to find the solution of the same equation with different conditions. The original NDSolveValue solutions work fine but when it ...
3
votes
2answers
82 views

Show Factorial instead of Gamma in the result of RSolve

Now I wanted to solve for a recurring function with RSolve. Here's how I tried: ...
3
votes
1answer
91 views

Finite Difference - heat generation in a square (square within a square)

I am trying to develop a code that will allow me to solve the following problem. Lets say that I have a square block, inside that square block I have a smaller square box thats at a certain ...
3
votes
1answer
73 views

Random Variable in Recurrence Function

Following the previously published question, I'm looking for the solution of RecurrenceTable with explicit random variable. For example, something like ...
3
votes
2answers
603 views

RecurrenceTable with a recurrence function of two variables

I have to calculate the cool down process of a regenerative heat exchanger. I solved the problem in Excel before, but now I want to do it with Mathematica. I want to create a table with the variables ...
3
votes
2answers
396 views

Nonlinear recurrence relation

I entered the following recurrence relation into RSolve and it just gave the question back to me. $$a_{n+1} = \frac{(1+a_n +(a_{n−1})^3)} 3$$ Is there another way to get a formula for the nth term?
3
votes
2answers
267 views

RSolve with assumptions

I want to solve the following non-linear first order difference equation: $$ k_{t+1} = (1-a) k_t^{a} $$ with $0 < a < 1$, and $k_0 = \text{const.}$ . Of course, the solution is easily found, ...
3
votes
0answers
80 views

How to combine DifferenceRoot objects for odd and even-indexed terms

I'm trying to analyze a certain infinite sequence $S$, indexed by positive integers starting from $1$. It can be split into two subsequences: $S^{odd}$ by removing all even-indexed elements and ...
3
votes
0answers
129 views

An recurrence equation that can be solved in version 7 but not in version 8

Backslide introduced in 8.0, fixed at least in 9.0.1. A friend of mine showed me this code sample: ...
2
votes
2answers
208 views

How can I express this recurrence function in Mathematica?

\begin{equation} r_{i+1}=r_{i}+V_{i}{\Delta}t \end{equation} Where \begin{equation} V_{i}=(0.01, 0.40) \end{equation} and ${\Delta}t$ is the difference in time. My implementation (this is for ...
2
votes
3answers
126 views

Computing a sequence by recurrence

I am trying to create a list where each element of the list is calculated using the sum of all previous elements of the list and also uses the $N$th element of another list. The initial value of the ...
2
votes
1answer
116 views

How to solve coupled complex-valued equations faster and get simplified solution

I have six coupled equations. How to solve these equations symbolically faster in Mathematica. In these equations, variables p,q,x,y are ...
2
votes
3answers
97 views

Permutations and recurrence equation

First problem: I need to show all Permutations with length 4 of this elements {x, y, z, w, t} but with this condition: the element z mustn't be after element t. I ...
2
votes
2answers
138 views

Why does Mathematica not evaluate the difference recurrence using `RSolve` or `RecurrenceTable`?

Mathematica does not evaluate the following recurrence of two variables using RSolve and only echoes the input: ...
2
votes
1answer
96 views

Speeding up a differential equation system with nested modulo periodic events

I have the following differential equation system with periodic Events : ...
2
votes
1answer
83 views

RSolve and a recurrence relation

In the book An Introduction to the Analysis of Algorithms there is a recurrence relation for the number of compares used in Quicksort algorithm: $$ C_N=N+1+\sum_{0≤k≤N-1}\frac{1}{N}(C_k+C_{N−k−1}) $$ ...
2
votes
1answer
95 views

Getting the first $n$ terms of a recursion relation

Given a recursion relation, I would like to generate the first $n$ terms (automatically) without having to go term by term. Background Each term is of the form $a_{n}(t)$, that is, $a_{n}$ is a ...
2
votes
2answers
141 views

RSolve with conditional equation

I have the following conditional equation: $f(n) = 1+2 f(n/2)$ when n is even $f(n) = 1+f((n-1)/2)+f((n+1)/2)$ when n is odd I want to represent this as a conditional equation to be used in ...
2
votes
1answer
155 views

How can I get DSolve to solve my ODE?

I evaluated DSolve[{x[t] x''[t] + x'[t]^2 - g x[t] == 0, x[0] == 0, x'[0] == 0}, x, t] and the output was {}. But this ODE ...
2
votes
1answer
130 views

RSolve: Errors while solving Partial Difference Equation

I tried solving following partial difference equation in Mathematica by typing ...
2
votes
0answers
45 views

Why doesn't RSolve work here? [closed]

RSolve[ {θ[n + 1] == θ[n] + 0.25*ϕ[n], ϕ[n + 1] == ϕ[n] - 2.5*Sin[θ[n]], θ[0] == π/2, ϕ[0] == 0}, {θ[n], ϕ[n]}, n] What's wrong with my code? I know ...
2
votes
0answers
126 views

RSolve and incomplete gamma function

I am interested in how Mathematica uses the incomplete gamma function to solve difference equations. For example, if we have this inhomogeneous, 2nd order equation and use ...
1
vote
1answer
128 views

MAP a business cycles model

There is a Business Cycles Map on Phaser 3.0 and I've trying to run it on Mathematica the last four days without succes. The equations are $$ x_1 = (1.0 - a s)x_1 + a \left[s m + g \left(\frac{s ...
1
vote
1answer
48 views
1
vote
1answer
96 views

Computing a Sequence using Recurrence

I posted a similar question earlier today, but unfortunately I realized after several hours that I over-simplified the problem for this forum. There is a very good answer to the Computing a sequence ...
1
vote
1answer
91 views

How do I use RSolve to solve the problem below?

RSolve can't compute the solution. Appears that Mathematica is having a convergence problem. I don't know if the boundary conditions are properly defined... ...
1
vote
2answers
129 views

Solving matrix equation

How can I solve the forward equation S[t]=C+A*S[t+1]*(Inverse((I-B*S(t+1)))*A for t=1 to t=10 and ...
1
vote
1answer
100 views

Limit as n -> ∞ for RSolve with Multivariate Recurrence Relation

I'm trying to solve the following when n -> ∞: ...
1
vote
1answer
225 views

Partial differential equation, Finite difference Method

For my research project I am trying to solve following partial differential equation in Mathematica (V[r]*ψ[u, v])/4 + Derivative[1, 1][ψ][u, v] == 0 with ...
1
vote
1answer
70 views

Finite Difference equations - slow computations

The following code I have written runs very slow, in fact, so slow that it does not finish computation. ...
1
vote
2answers
59 views

RecurrenceTable - does it have a solution search option like FindRoot?

I have a second-order non-linear difference equation like $$0=f(x_{t+2},x_{t+1},x_t)$$ where $x_{t+2}$ is implicitly determined. I have tried the code like ...