blob: 581c12430fe3b223100d04b75eac488519c913ad (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
//Initilization Of Variables
a=0 //Lower Limit oF the Integral
b=8 //Upper Limit of the Integral
n=10 //Interval of the integral
//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=75*t^2')
Rr=Trap_Composite1(f,a,b,n)/(2*8) //lb
//Moment calculations
M=Trap_Composite1(f,a,b,n) //ft-lb
//Result
clc
printf('The reaction is %f lb and Moment is %f lb-ft',Rr,M)
//Decimal accuracy causes discrepancy in answers
|
