(%i1) x0:1;

(%i2) f(x):=x^2-2*x;

(%i3) f(x0-h);

(%i4) ratsimp(%);

(%i5) f(x0+h);

(%i6) ratsimp(%);

(%i7) x0:2;

(%i8) f(x)=x^4-8*x^3+20*x^2-16*x;

(%i9) f(x0+h);

(%i10) f(x):=x^4-8*x^3+20*x^2-16*x;

(%i11) f(x0+h);

(%i12) ratsimp(%);

(%i13) f(x0-h);

(%i14) ratsimp(%);

(%i15) x0:-1;

(%i16) f(x):=(x*(x+2))/(x+1)^4;

(%i17) f(x0+h);

(%i18) ratsimp(%);

(%i19) f(x0-h);

(%i20) ratsimp(%);

(%i21) x0:1;

(%i22) f(x):=2^x-2^(2-x);

(%i23) f(x0+h);

(%i24) ratsimp(%);

(%i25) f(x0-h);

(%i26) ratsimp(%);

(%i27) x0:-0.5;

(%i28) y0:0;

(%i29) y0-f(x0-h);

(%i30) f(x):=2*x^3+3*x^2+x;

(%i31) y0-f(x0-h);

(%i32) ratsimp(%);

(%i33) f(x0+h)-y0;

(%i34) ratsimp(%);

(%i35) x0:2;

(%i36) y0:0;

(%i37) f(x):=(x-2)/(x^2-4*x);

(%i38) y0-f(x0-h);

(%i39) ratsimp(%);

(%i40) f(x0+h)-y0;

(%i41) ratsimp(%);

(%i42)