blob: aab8193adbbb08f7a12e2410a423332adaa22a5a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
//Eg-11.5
//pg-477
clear
clc
//defining an inline function for simplicity and integrating it from x = 0 to 0.9
deff('out = func(in)','out = 0.3 + 1.6*in + 0.027*in^2')
h = (0.9-0)/3; //Stepsize
b = 0.9;
a = 0;
//Integrating using the simpsons 1/3rd rule
I = (3*h/8)*(func(a) + 3*func(a+h) + 3*func(a+2*h) + func(b));
printf('The value of the integration is %f\n',I)
printf(' Since the function to be integrated is a quadratic, the value of the integral is exact')
|
