blob: b38a5bb8370cba9949af64926641e27142651df3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
//Example No. 12_04
//Simpon's 1/3 rule
//Pg No.382
clear ;close ;clc ;
deff('F = f(x)','F = sqrt( sin(x) )');
x0 = 0 ;
xa = %pi/2 ;
//case(a) n = 4
n = 4 ;
h = ( xa-x0 )/n
I = 0
for i = 1:n/2
I = I + f(x0 + (2*i-2)*h) + 4*f(x0 + (2*i-1)*h) + f(x0 + 2*i*h) ;
end
I = h*I/3
disp(I,'Integral value for n = 4 is',h,'h = ')
//case(b) n = 6
n = 6
h = ( xa-x0 )/n
I = 0
for i = 1:n/2
I = I + f(x0 + (2*i-2)*h) + 4*f(x0 + (2*i-1)*h) + f(x0 + 2*i*h) ;
end
I = h*I/3
disp(I,'Integral value for n = 6 is',h,'h = ')
|
