function [BSPL,DBSPL]=bsplfunc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% loads data (up to degree 5):
%      (*) inner integral coefficients of b-splines for inner cells
%      (*) b-spline equations to compute integrals for boundary cells
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
global n;
cwd=pwd;cd('input');
if n==1,load inoutinteg1Dn1.mat;
elseif n==2,load inoutinteg1Dn2.mat;
elseif n==3,load inoutinteg1Dn3.mat;
elseif n==4,load inoutinteg1Dn4.mat;
elseif n==5,load inoutinteg1Dn5.mat;
end, cd(cwd);
