function y=fd(x)
  y=1+sin(x)
endfunction