tic;
v=fdm('primerG1.m',2^-13);
toc;
pause(10)
tic;
v=fdms('primerG1.m',2^-13);
toc;
pause(10)
tic; v=fdm1('primerG1', 2^(-13)); toc;

