2011過去問 - suchikeisan @ ウィキ

[1]

 (1.1.1) 

 #[1]1. 

NumIntTpz:=proc(f,N,a,b)

local S, h, xi, i; 
h:=(b-a)/N: 
S:=f(a)/2; 
for i from 1 to N-1 do 
  xi:=a+i*h; 
  S:=S+f(xi); 
end do: 
S:=S+f(b)/2; 
return S*h; 

end proc:

 #[1]2. 

dX2:=[]:
for i from 1 to nops(N) do

dX2:=[op(dX2),abs(s-NumIntTpz(y,N[i],a,b))]; 

printf("%20.15f\n", dX2[i]);
end do:

 0.409541846000000 
 0.145883586500000 
 0.051773510300000 
 0.018339615630000 
 0.006490230280000 
 0.002295740170000 
 0.000811865000000 
 0.000287064005000 

 (1.1.2) 

[2]

 (1.2.1) 

[3]

(1.4.1) 

[4]直線による近似

(1.5.1) 

(1.5.2) 

[4]曲線による近似

 (1.7.1)