Function funspace::cheb_neumann

pub fn cheb_neumann<A: FloatNum>(n: usize) -> BaseR2r<A>

Function space with Neumann boundary conditions

 u'(-1)=0 and u'(1)=0

 \phi_k = T_k - k^{2} \/ (k+2)^2 T_{k+2}

Example

Transform array to function space.

use funspace::cheb_neumann;
use funspace::traits::BaseTransform;
use ndarray::Array1;
let ch = cheb_neumann::<f64>(10);
let mut y = Array1::<f64>::linspace(0., 9., 10);
let yhat: Array1<f64> = ch.forward(&mut y, 0);