Chebyshev interval inclusion function

This code is used to interval uncertainty analysis based on Chebyshev interval inclusion function proposed by Jinglai Wu.

% function [y_lb,y_ub]=CI_reg(fun_name,a,b,k,K,Expansion)
% Input
% fun_name the called function name
% a the lower bound vector of interval input
% b the upper bound vector of interval input
% k the order of CI expansion
% K the scanning (validation) points of each interval variable
% Expansion expansion type of Chebyshev polynomials-'full' or 'partial'

% Output
% y_lb lower bound of response
% y_ub upper bound of response

% Example
%[y_lb1,y_ub1]=CI_reg(@double_pendulum,[0.99 1.98]',[1.01 2.02]',4,10,'full');