# Easy Octave Algorithm 1

This page contains an Octave version of Algorithm 1, on page 55 of the B-series book. This function is intended to run equally well on Matlab, SciLab, as well as Octave.

Its purpose is to evaluate $a_1$, $a_2$, $\dots$, up to a maximum order of pmax which satisfy the equation $$ a_1+a_2x+a_3x^2 + \cdots = (1-x)^{-a_1} (1-x^2)^{-a_2} \cdots $$ in the sense that the coefficients of each power of $x$ on the two sides must agree.

#### Octave function definition

function Ntrees = alg1(pmax) Ntrees = ones(1, pmax) for i = 2: pmax - 1 for j = pmax: -1: 1 a = Ntrees(i) for k = 1: (j - 1)/i Ntrees(j) = Ntrees(j) + a*Ntrees(j - i*k) a = a*(Ntrees(i) + k)/(k + 1) end end end end

$\to$ B-series book

$\to$ Homepage

- Log in to post comments