Cześć potrzebuję pomocy a mianowicie pętli która będzie podstawiała argument "h" do wzoru

function [y] = wzorA1()
h= 15;
x= (1/sqrt(3))
y = (3*zad1(x) - 14*zad1(x+h) + 26 *zad1(x+2*h) - 24* zad1(x+3*h) + 11* zad1(x+4*h) - 2*zad1(x+5*h))/(h^4);
end
 

w miejscu "zad1" jest pochodnia "6x5 - 2x</sup>3 + 3x2 - 47x - 7" której dokładny wynik wynosi "415,6921938" . Wzór w kodzie jest na przybliżoną wartość czwartej pochodnej , problem w tym że trzeba dobrać taki argument "h" by błąd różnicy po między dokładną wartością a przybliżoną był mniejszy niż 10 do -13 " 10 </sup>(-13).