moler < function(){
Fd < matrix(0,nrow=16,ncol=2)
f < function(x) exp(x)/(sin(x)^3 + cos(x)^3)
fcs < function(x,h) Im(f(x + 1i*h))/h
fdf < function(x,h) (f(x+h)  f(xh))/(2*h)
x0 < pi/4
h < 10^(seq(1,16))
Fd[,1] < fcs(x0,h)
Fd[,2] < fdf(x0,h)
return(Fd)
}
