X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=legendre_p_tilde.m;h=3bb30be664003fcab28ca19259d7b8f4821f045b;hp=2c9d2ab08916107fe11ca42ba863f32e4341117f;hb=b12c6c2a4bf4cef29b2e08b743c92889505c7ed9;hpb=e1b71b4ca7cfa08ac76744a17a3778d4ccfaa7e2 diff --git a/legendre_p_tilde.m b/legendre_p_tilde.m index 2c9d2ab..3bb30be 100644 --- a/legendre_p_tilde.m +++ b/legendre_p_tilde.m @@ -1,24 +1,24 @@ function P_tilde = legendre_p_tilde(n, a, b) - ## Return the `n`th Legendre polynomial scaled to the interval [a,b]. - ## - ## INPUT: - ## - ## * ``n`` - The index of the polynomial that we want. - ## - ## * ``a`` - The left endpoint of the interval. - ## - ## * ``b`` - The right endpoint of the interval. - ## - ## OUTPUT: - ## - ## * ``P_tilde`` - A polynomial function of one argument. - ## + % Return the `n`th Legendre polynomial scaled to the interval [a,b]. + % + % INPUT: + % + % * ``n`` - The index of the polynomial that we want. + % + % * ``a`` - The left endpoint of the interval. + % + % * ``b`` - The right endpoint of the interval. + % + % OUTPUT: + % + % * ``P_tilde`` - A polynomial function of one argument. + % if (n < 0) - ## Can't do anything here. Return nothing. + % Can't do anything here. Return nothing. P = NA; else - ## Compute the Legendre polynomial over [-1,1] and mangle it to fit - ## the interval [a,b]. + % Compute the Legendre polynomial over [-1,1] and mangle it to fit + % the interval [a,b]. P = legendre_p(n); P_tilde = @(x) P( (2/(b-a)).*x + 1 - (2*b)/(b-a) ); end