function F = dsinfun(x,r)
% m-file for general function
F = r*pi*cos(pi*x);