clear all, clc
facit=-1.68043;
y=@(x) ((x-7).^(5/2)+2*sin(pi*sqrt(x)))./...
    (sqrt(x+4*log(x-2*pi))-1);

for h=[0.04 0.02 0.01 0.005]
    fp1=1/h*(y(7+h)-y(7));
    fp2=1/(2*h)*(-y(7+2*h)+4*y(7+h)-3*y(7));
    
    rad=[h fp1 fp1-facit fp2 fp2-facit];
    disp(rad)
end
