ni får ursäkta det här men har bråttom bråttom bråttom

function y = funktion(x)
y=exp(-x./3)./(2-cos(x.*pi));
end

quad(@funktion,0,2.6)

x=linspace(0,2.6,100)’;
f=funktion(x);
fi=linspace(0,2*pi,30);
X=x*ones(size(fi));
Y=f*cos(fi);
Z=f*sin(fi);
mesh(X,Y,Z);

function y = uppgift1(h)
b=2.6;
a=0;
n=(b-a)/h;
x=linspace(0,2.6,(n+1));
f=(pi.*(funktion(x).^2));
I=h*(sum(f)-(f(1)+f(n+1))/2);
n=2*n;
h=h/2;
a=linspace(0,2.6,(n+1));
b=(pi*(funktion(a).^2));
U=h.*(sum(b)-(b(1)+b(n+1))/2);
S=U+((I-U)./3)
A=I-S
I=S
end

function y = funktion(x)

y = 3*x.*exp(-x.^3);

end

function y = quadfunk(h)
b=4;
a=0;
n=(b-a)/h;
x=linspace(0,2.6,(n+1));
A=1;
while abs(A) > 1.E-6
f=funktion(x);
I=h*(sum(f)-(f(1)+f(n+1))/2);
n=2*n;
h=h/2;
a=linspace(0,2.6,(n+1));
b=funktion(a);
U=h.*(sum(b)-(b(1)+b(n+1))/2);
S=U+((I-U)./3)
A=I-S
I=S
end
end