function FunS=FpexRlxNumSin(y,al,del);  
% sumS fr FPEX relaxation wth IC DEP k=2n={2,4,6,8,10,..} 
% indefinte integral to X (x^2*al-1)*sin(del*x) dx del=m*pi/(b-a)
% . for matix dimensions D. Bakewell 17:45 Su/15/8/1999
format long;
%FunS=sin(del*x);
yal2m1=y.^(2*al-1);
FunS=yal2m1.*sin(y*del);