function r=myfun25
fun=@(x)1./(exp(x)+exp(-x));
r=quad(fun,0,1);
end
function r=myfun26
x=[2.04 3.7 3.96 6.27 8.66 9.57 10.07 11.25 11.28 11.97];
y=[41.898 99.07 111.498 242.084 435.068 520.894 642.71 705.492 711.092 794.855];
fun=polyfit(x,y,2);
x0=3:0.5:40;
r=polyval(fun,x0);
end
function r&#