blob: 691ca3f32afa3a1547c8c10fa3b48da18a57fb39 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
//Initilization Of Variables
a=0 //Lower Limit oF the Integral
b=0.3 //Upper Limit of the Integral
n=10 //Interval of the integral
g=9.8 //m/s^2
rho=1000 //kg/m^3
//Calculation
//Using Trapezoidal Rule for Intergration
function[I1]=Trap_Composite1(f,a,b,n)
h=(b-a)/n
t=linspace(a,b,n+1)
I1=(h/2)*((2*sum(f(t)))-f(t(1))-f(t(n+1)))
endfunction
deff('[y]=f(t)','y=(g*0.6*rho*1.2*t)-(0.6*g*rho*t^2)')
B=Trap_Composite1(f,a,b,n)/(2*b) //N
//Result
clc
printf('The value of B is %f N',B)
|
