%Störningsräkning på ekv exp(-ax)-x^2=0, a=1.0 +- 0.05
for a=[0.95,1,1.05]
x=1;h=1;
while abs(h)>1e-4
    f=exp(-a*x)-x^2;
    fprim=-a*exp(-a*x)-2*x;
    h=f/fprim;
    x=x-h;
end
rot=x;
resultat=[a,rot]
end
