h[7]=-0.0105974017850890; h[6]= 0.0328830116668852; h[5]= 0.0308413818355607; h[4]=-0.1870348117190931; h[3]=-0.0279837694168599; h[2]= 0.6308807679398597; h[1]= 0.7148465705529154; h[0]= 0.2303778133088964; isgn=1; for(k=0;k<8;k++) { gg[k]=isgn*h[7-k]; isgn=-isgn; } for(k=0;k<8;k++) { g[k]=gg[7-k]; hh[k]=h[7-k]; } return; } float fun(float x) { float pi=3.1415926; float yx=30*exp(-x/40)*sin(2*pi*x/40); return(yx); }