% x is a vector where x(i) = 10^(-i) format long e x = zeros(1,8); for i = 1:8 x(i) = 1/2+10^(-i); end a = zeros(1,8); for j = 1:8 a(i) = fsin(x(i))-sin(x(i)); end b = zeros(1,8); for j = 1:8 b(i) = fcos(x(i))-cos(x(i)); end a b