function y = f541(x)
y = cos(x)/2 - x*sin(x);