# -*- coding: utf-8 -*-
"""
Created on Sat Jan 14 18:13:19 2017

@author: Mariusz
"""
from numpy import *

def chebyshev(a, b, N):
    # return Chebyshev's interpolation points on the interval [a,b]
    I = arange(1, N+1, 1)
    X = (b + a)/2 + (b - a)/2*cos((2*I - 1)*pi/(2*N))
    return X