
function fsin = fsin(q)

fsin = sin(q) - (q - pi/4)*cos(q) - sin(pi/4);
