1.

Derive the equations of motion for the 5-storey

shear frame structure.
(a) Provide the mass matrix m.
(b) Provide the stiffness matrix k.

MASS MATRIX (M)

>>M= [245 0 0 0 0;0 147 0 0 0;0 0 196 0 0;0 0 0 49 0;0 0 0 0 98]
M=
245.0000 0 0 0 0
0 147.0000 0 0 0
0 0 196.0000 0 0
0 0 0 49.0000 0
0 0 0 0 98.0000

STIFFNESS MATRIX (K)

>> K=[328.125 -93.75 0 0 0;-93.75 187.5 -93.75 0 0;0 -93.75 281.25 -
187.5 0;0 0 -187.5 468.75 -281.25;0 0 0 -281.25 281.25]

K=
328.125 -93.750 0 0 0
-93.750 187.500 -93.750 0 0
0 -93.750 281.250 -187.500 0
0 0 -187.500 468.750 -281.250
0 0 0 -281.250 281.250
2. Solve the general eigenvalue problem to determine
the mode shapes and Eigen frequencies of the frame
structure.
(a) Calculate the circular natural frequencies ωi and the
natural periods Ti of each mode.
(b) Determine the mode shapes ϕi and normalize each
mode such that the maximum value is 1.0.

>> KT=K/1

KT =
328.125 -93.750 0 0 0
-93.750 187.500 -93.750 0 0
0 -93.750 281.250 -187.500 0
0 0 -187.500 468.750 -281.250
0 0 0 -281.250 281.250

>> [Eigenvectors, Eigenvalues] = eig(K,M)

Eigenvectors =
-0.0094 -0.0400 0.0439 0.0215 0.0000
-0.0306 -0.0521 -0.0311 -0.0468 -0.0007
-0.0473 0.0046 -0.0276 0.0440 0.0119
 -0.0511 0.0289 0.0206 -0.0104 -0.1282
-0.0528 0.0409 0.0469 -0.0428 0.0413

Eigenvalues =
0.0926 0 0 0 0
0 0.8413 0 0 0
0 0 1.6103 0 0
0 0 0 2.1698 0
0 0 0 0 11.7719

>> omegas=sqrt(Eigenvalues)
omegas =
0.3043 0 0 0 0
0 0.9172 0 0 0
0 0 1.2690 0 0
0 0 0 1.4730 0
0 0 0 0 3.4310
>> T=2*pi./omegas
T=
20.6471 Inf Inf Inf Inf
Inf 6.8502 Inf Inf Inf
Inf Inf 4.9514 Inf Inf
Inf Inf Inf 4.2655 Inf
Inf Inf Inf Inf 1.8313
>> T1= 2*pi./omegas(1,1)
T1 = 20.6471
>> T2= 2*pi./omegas(2,2)
T2 = 6.8502
>> T3= 2*pi./omegas(3,3)
T3 = 4.9514
>> T4= 2*pi./omegas(4,4)
T4 = 4.2655
>> T5= 2*pi./omegas (5,5)
T5 = 1.8313
>> phi1 =Eigenvectors(:,1)./Eigenvectors(1,1)
phi1 =
1.0000
3.2580
5.0429
5.4472
5.6288
>> phi2 =Eigenvectors(:,2)./Eigenvectors(1,2)
phi2 =
1.0000
1.3014
-0.1140
-0.7214
-1.0206
>> phi3 =Eigenvectors(:,3)./Eigenvectors(1,3)
phi3 =
1.0000
-0.7083
-0.6282
0.4693
1.0693
>> phi4 =Eigenvectors(:,4)./Eigenvectors(1,4)
phi4 =
1.0000
-2.1705
2.0437
-0.4847
-1.9871
>> phi5 =Eigenvectors(:,5)./Eigenvectors(1,5)
phi5 =
1.0e+03 *

0.0010
-0.0273
0.4477
-4.8242
1.5553
>> phi(:,1) = phi1; phi(:,2) = phi2; phi(:,3) = phi3; phi(:,4) =
phi4;phi(:,5)=phi5;
>> phi
phi =
1.0e+03 *

0.0010 0.0010 0.0010 0.0010 0.0010

0.0033 0.0013 -0.0007 -0.0022 -0.0273
0.0050 -0.0001 -0.0006 0.0020 0.4477
0.0054 -0.0007 0.0005 -0.0005 -4.8242
0.0056 -0.0010 0.0011 -0.0020 1.5553

3. Determine the analytical solution u(t) for the free

vibration response based on the modal expansion
technique.
(a) Calculate the generalized mass matrix.
(b) Determine the vector of the modal contribution
factors and its temporal derivative at t(0) s, i.e. q(0)
and p(0). The initial conditions of the frame structure
are given as:
(c)Provide the expression for the free vibration
response of the system.
(d) Evaluate the analytical solution for the time interval
0 d’ t d’ r5T1s1 using the time increment t1 and
determine the maximum and minimum values for the
displacement of floor 5, i.e maxru ∆t1 5 ptqs and minru
∆t1 5 ptqs.
>> Mstar=phi'*M*phi
Mstar =
1.0e+09 *

0.0000 0.0000 -0.0000 0.0000 -0.0000

0.0000 0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 0.0000 0.0000
-0.0000 -0.0000 -0.0000 0.0000 1.4168

>> phin=phi./4824.2
phin =
0.0002 0.0002 0.0002 0.0002 0.0002
0.0007 0.0003 -0.0001 -0.0004 -0.0057
0.0010 -0.0000 -0.0001 0.0004 0.0928
0.0011 -0.0001 0.0001 -0.0001 -1.0000
0.0012 -0.0002 0.0002 -0.0004 0.3224

>> u0=[-5;-2;-2;4;6]
u0 =
-5
-2
-2
4
6

>> coeff=inv(Mstar)*phi'*M*u0
coeff =
0.0192
-3.6925
-0.0960
-1.2306
-0.0001
>> syms t w1 w2 w3 w4 w5 Y
>> Y= coeff.*[cos(w1*t); cos(w2*t); cos(w3*t); cos(w4*t); cos(w5*t)]

Y=
(691190818189393*cos(t*w1))/36028797018963968
-(8314774808583819*cos(t*w2))/2251799813685248
-(864455435986445*cos(t*w3))/9007199254740992
-(5541984676715875*cos(t*w4))/4503599627370496
-(5201733627203681*cos(t*w5))/36893488147419103232
>> Y= vpa(Y,5)
Y=
0.019184*cos(t*w1)
-3.6925*cos(t*w2)
-0.095974*cos(t*w3)
-1.2306*cos(t*w4)
-0.00014099*cos(t*w5)

>> v= vpa(phi*Y,5)
v=
0.019184*cos(t*w1) - 3.6925*cos(t*w2) - 0.095974*cos(t*w3) -
1.2306*cos(t*w4) - 0.00014099*cos(t*w5)
0.062503*cos(t*w1) - 4.8053*cos(t*w2) + 0.067978*cos(t*w3) +
2.6709*cos(t*w4) + 0.003844*cos(t*w5)
0.096745*cos(t*w1) + 0.421*cos(t*w2) + 0.060287*cos(t*w3) -
2.5149*cos(t*w4) - 0.063125*cos(t*w5)
0.1045*cos(t*w1) + 2.6639*cos(t*w2) - 0.045041*cos(t*w3) +
0.59648*cos(t*w4) + 0.68018*cos(t*w5)
0.10799*cos(t*w1) + 3.7687*cos(t*w2) - 0.10262*cos(t*w3) +
2.4452*cos(t*w4) - 0.21928*cos(t*w5)

>> v=[1; 1; 1; 1; 1]
v=
1
 1
1
1
1

>> Y1=(phi1'*M*v)/(phi1'*M*phi1)
Y1= 0.2230
>> Y2=(phi2'*M*v)/(phi2'*M*phi2)
Y2 = 0.4464
>> Y3=(phi3'*M*v)/(phi3'*M*phi3)
Y3 = 0.2805
>> Y4=(phi4'*M*v)/(phi4'*M*phi4)
Y4 = 0.0501
>> Y5=(phi5'*M*v)/(phi5'*M*phi5)
Y5 = 1.4052e-08
>> check = Y1*phi1 + Y2*phi2 +Y3*phi3 +Y4*phi4+Y5*phi5
check =
1.0000
1.0000
1.0000
1.0000
1.0000

MAX= 5.999691
MINI= 0.025121

4. Determine the numerical solution ~u(t) for the free vibration

response based on Newmark’s method (use the constant average
acceleration method).
(a) Compute the numerical solution with the time increment t1 in the
same time interval (0 d’ t d’ r5T1s) and determine the maximum and
minimum values for the displacement of floor 5, i.e maxr~ ∆t1 5 ptqs
and minr~ ∆t1 5 ptqs.
(b) Plot the solution for u5ptq and u~5ptq in one diagram.
(c) Discuss the differences in the results.
(d) Determine a suitable time increment t2 for the numerical solution
such that the relative error in the maximum displacement is less than
0.1%.

(a) MAX= 22.45915

MINI= 0.09167
(b)
8

0
0 50 100 150 200 250 300
-2

-4

-6

-8
(c) It is a Sine wave. If the difference in the displacement (u) of the
structure is significant, then it may pose issue in the structural
stability.
(d) ERROR- 0.11454
APPENDIX-

loop P^ U U. U.. U AB

0 0 6 0 -17.2194 6

1 278841.9 5.928247 -1.56727 -17.0135 5.928247

2 268797.6 5.714703 -3.09705 -16.4006 5.714703

3 252324.3 5.364476 -4.55275 -15.3955 5.364476

4 229815.9 4.885943 -5.89956 -14.0222 4.885943

5 201810.9 4.290549 -7.10527 -12.3134 4.290549

6 168979 3.592534 -8.14104 -10.3102 3.592534

7 132105.5 2.808595 -8.98209 -8.06038 2.808595

8 92072.28 1.957479 -9.60831 -5.61777 1.957479

9 49836.94 1.059546 -10.0047 -3.04079 1.059546

10 6409.606 0.13627 -10.1618 -0.39108 0.13627

11 -37171 -0.79027 -10.0759 2.267981 0.790265

12 -79862.6 -1.6979 -9.74899 4.872796 1.697899

13 -120644 -2.56492 -9.18889 7.361066 2.564923

14 -158540 -3.3706 -8.40902 9.673275 3.370599

15 -192644 -4.09566 -7.42802 11.75412 4.095658

16 -222140 -4.72276 -6.26936 13.55383 4.722758

17 -246324 -5.2369 -4.96075 15.02937 5.236901

18 -264615 -5.62579 -3.53349 16.14544 5.625788

19 -276578 -5.88012 -2.02171 16.87534 5.880119

20 -281926 -5.99381 -0.46158 17.20162 5.99381

21 -280530 -5.96414 1.109585 17.11648 5.964143

22 -272425 -5.79183 2.654214 16.62195 5.791827

23 -257804 -5.48098 4.135361 15.72986 5.480983

24 -237017 -5.03905 5.5176 14.46155 5.039047

25 -210561 -4.47659 6.767869 12.84735 4.476587

26 -179069 -3.80706 7.856267 10.92587 3.807058

27 -143294 -3.04647 8.756759 8.743063 3.046472

28 -104092 -2.21302 9.44781 6.351144 2.213021

29 -62400 -1.32664 9.91289 3.807321 1.32664

30 -19215.6 -0.40853 10.14088 1.172434 0.408528

31 24428.43 0.519354 10.12631 -1.49049 0.519354

32 67488.17 1.434815 9.869553 -4.11777 1.434815

33 108933.7 2.315958 9.376735 -6.64656 2.315958

34 147773.9 3.141709 8.659646 -9.01638 3.141709

35 183079.6 3.892317 7.735438 -11.1706 3.892317

36 214006.4 4.549829 6.626215 -13.0575 4.549829

37 239814.7 5.09852 5.358508 -14.6322 5.09852

38 259887.2 5.525265 3.962637 -15.8569 5.525265

39 273743.7 5.819858 2.471988 -16.7024 5.819858

40 281052.9 5.975253 0.922216 -17.1484 5.975253

41 281639.9 5.987734 -0.64961 -17.1842 5.987734

42 275490.7 5.857001 -2.20591 -16.809 5.857001

43 262752.4 5.586182 -3.70944 -16.0318 5.586182

44 243729.7 5.181753 -5.12425 -14.8711 5.181753

45 218877.5 4.653389 -6.4165 -13.3548 4.653389

46 188790.2 4.013726 -7.55528 -11.519 4.013726

47 154187.4 3.278063 -8.51336 -9.40771 3.278063

48 115896.9 2.463997 -9.26781 -7.07142 2.463997

49 74834.35 1.590997 -9.8006 -4.566 1.590997

50 31981.93 0.679944 -10.099 -1.95137 0.679944

51 -11635.4 -0.24737 -10.1558 0.709932 0.247372

52 -54974.5 -1.16877 -9.96975 3.354254 1.168771

53 -96998.7 -2.06222 -9.54523 5.918349 2.062216

54 -136703 -2.90634 -8.89241 8.34089 2.906337

55 -173137 -3.68094 -8.0269 10.56394 3.680945

56 -205431 -4.36751 -6.9694 12.53432 4.367513

57 -232811 -4.94962 -5.74521 14.2049 4.949619

58 -254623 -5.41334 -4.38361 15.53574 5.413342

59 -270344 -5.74759 -2.91716 16.495 5.74759

60 -279600 -5.94437 -1.38094 17.05973 5.944368

61 -282168 -5.99897 0.188304 17.21643 5.99897

62 -277988 -5.91009 1.753048 16.96136 5.91009

63 -267158 -5.67985 3.275864 16.3006 5.679854

64 -249939 -5.31377 4.720328 15.24997 5.313768

65 -226742 -4.82059 6.051893 13.8346 4.820589

66 -198122 -4.21211 7.238709 12.08833 4.212112

67 -164762 -3.50289 8.252392 10.05294 3.502891

68 -127463 -2.70989 9.068696 7.777103 2.709888

69 -87114.3 -1.85207 9.668097 5.315255 1.852071

70 -44682.3 -0.94996 10.03626 2.726278 0.949956

71 -1181.59 -0.02512 10.16437 0.072094 0.025121

72 42347.37 0.900315 10.04938 -2.58381 0.900315

73 84863.47 1.804218 9.69403 -5.17792 1.804218

74 125349.8 2.664968 9.106819 -7.64819 2.664968

75 162838.1 3.461978 8.301793 -9.93552 3.461978

76 196431.6 4.176185 7.298206 -11.9852 4.176185

77 225327 4.790507 6.120063 -13.7483 4.790507

78 248833 5.29025 4.795541 -15.1825 5.29025

79 266387.5 5.663463 3.35632 -16.2536 5.663463

80 277570.5 5.901218 1.836824 -16.9359 5.901218

81 282114.8 5.997829 0.273395 -17.2132 5.997829

82 279911.4 5.950986 -1.29657 -17.0787 5.950986

83 271013.2 5.761808 -2.83553 -16.5358 5.761808

84 255633 5.43482 -4.30667 -15.5974 5.43482

85 234138.6 4.977844 -5.6748 -14.2859 4.977844

86 207044.1 4.401808 -6.9072 -12.6327 4.401808

87 174997.6 3.720491 -7.9744 -10.6774 3.720491

88 138765.5 2.950189 -8.85087 -8.46674 2.950189

89 99214.47 2.109324 -9.51564 -6.05354 2.109324

90 57290.44 1.218009 -9.95282 -3.49556 1.218009

91 13996.16 0.297562 -10.152 -0.85397 0.297562

92 -29632.9 -0.63 -10.1083 1.808042 0.630002

93 -72553.2 -1.5425 -9.82283 4.426812 1.542498

94 -113738 -2.4181 -9.30244 6.939703 2.418101

95 -152203 -3.23587 -8.55956 9.286611 3.235868

96 -187027 -3.97624 -7.61195 11.4114 3.97624

97 -217378 -4.62151 -6.48229 13.26326 4.62151

98 -242530 -5.15624 -5.19758 14.79789 5.156243

99 -261881 -5.56765 -3.78855 15.97859 5.567651

100 -274968 -5.84589 -2.28891 16.77711 5.845892

101 -281479 -5.98431 -0.73453 17.17437 5.984313

102 -281257 -5.9796 0.837422 17.16085 5.979602

103 -274309 -5.83187 2.389344 16.73688 5.831873

104 -260799 -5.54466 3.884119 15.9126 5.544658

105 -241052 -5.12483 5.285995 14.70773 5.124827

106 -215539 -4.58242 6.561441 13.15108 4.582422

107 -184872 -3.93042 7.679952 11.27989 3.930415

108 -149782 -3.1844 8.614776 9.138909 3.184402

109 -111110 -2.36222 9.343553 6.779344 2.362225

110 -69780.4 -1.48355 9.848854 4.257633 1.483548

111 -26781.9 -0.56939 10.11859 1.634088 0.569389

112 16857.24 0.358389 10.14632 -1.02854 0.358389

113 60093.16 1.277595 9.931361 -3.66657 1.277595

114 101891.8 2.166244 9.478871 -6.2169 2.166244

115 141253.4 3.003082 8.799668 -8.61854 3.003082

116 177236.5 3.768092 7.909996 -10.814 3.768092

117 208980.5 4.442978 6.831135 -12.7509 4.442978

118 235726.2 5.011597 5.588887 -14.3828 5.011597

119 256833.8 5.460351 4.212967 -15.6706 5.460351

120 271798.6 5.778505 2.736281 -16.5837 5.778505

121 280262.5 5.958449 1.19415 -17.1001 5.958449

122 282023.1 5.995882 -0.37654 -17.2076 5.995882

123 277038.4 5.889906 -1.93823 -16.9034 5.889906

124 265427.6 5.643056 -3.45356 -16.195 5.643056

125 247468.3 5.261237 -4.88629 -15.0992 5.261237

126 223590.1 4.753581 -6.20214 -13.6423 4.753581

127 194364.2 4.13223 -7.36966 -11.8591 4.13223

128 160489.4 3.412045 -8.36091 -9.79222 3.412045

129 122776.2 2.610252 -9.15219 -7.49116 2.610252

130 82126.38 1.746027 -9.72456 -5.01092 1.746027

131 39512.3 0.840041 -10.0643 -2.41083 0.840041

132 -4046.83 -0.08604 -10.1634 0.246916 0.086037

133 -47509.2 -1.01006 -10.0194 2.898759 1.010057

134 -89835.2 -1.90992 -9.63574 5.481271 1.909918

135 -130013 -2.7641 -9.02162 7.932682 2.764099

136 -167080 -3.55217 -8.19172 10.19436 3.552169

137 -200152 -4.25528 -7.16589 12.21221 4.255278

138 -228436 -4.85661 -5.96867 13.93798 4.856611

139 -251257 -5.34178 -4.62869 15.33038 5.341784

140 -268068 -5.69919 -3.178 16.35611 5.699194

141 -278468 -5.92029 -1.65131 16.99063 5.920292

142 -282207 -5.99979 -0.08511 17.21878 5.99979

143 -279196 -5.93579 1.483116 17.0351 5.935786

144 -269508 -5.72981 3.015871 16.44397 5.729811

145 -253374 -5.38679 4.476494 15.45954 5.386792

146 -231179 -4.91493 5.830049 14.10535 4.914932

147 -203456 -4.32552 7.044162 12.4138 4.325519

148 -170866 -3.63265 8.089794 10.42533 3.632648

149 -134189 -2.85289 8.941937 8.187512 2.852893

150 -94302.9 -2.0049 9.580208 5.753867 2.004903

151 -52161.2 -1.10896 9.989342 3.182603 1.10896

152 -8771.93 -0.18649 10.15955 0.535217 0.186493

153 34827.15 0.740434 10.08677 -2.12497 0.740434

154 77593.25 1.649652 9.772734 -4.73433 1.649652

155 118503.5 2.519413 9.224956 -7.23046 2.519413

156 156579.4 3.328916 8.456537 -9.55365 3.328916

157 190910.2 4.058799 7.485857 -11.6483 4.058799

158 220675 4.691604 6.336132 -13.4644 4.691604

159 245161.6 5.212197 5.034861 -14.9585 5.212197

160 263784.6 5.608125 3.613167 -16.0947 5.608125

161 276098.4 5.869919 2.105054 -16.8461 5.869919

162 281808.5 5.991319 0.546593 -17.1945 5.991319

163 280778.4 5.969419 -1.02494 -17.1316 5.969419

164 273032.7 5.804743 -2.57196 -16.659 5.804743

165 258756.7 5.501231 -4.05747 -15.788 5.501231

166 238291.8 5.066142 -5.44593 -14.5393 5.066142

167 212127.5 4.509883 -6.70413 -12.9429 4.509883

168 180889.5 3.845756 -7.80199 -11.0369 3.845756

169 145325.1 3.089648 -8.71324 -8.86697 3.089648

170 106284.9 2.259642 -9.41609 -6.48494 2.259642

171 64702.5 1.375591 -9.89372 -3.94781 1.375591

172 21572.59 0.458639 -10.1347 -1.31625 0.458639

173 -22073.3 -0.46928 -10.1333 1.346795 0.469283

174 -65191.2 -1.38598 -9.88956 3.977625 1.385981

175 -106750 -2.26953 -9.40926 6.513318 2.269529

176 -145755 -3.0988 -8.70391 8.893227 3.098796

177 -181275 -3.85395 -7.79038 11.06043 3.853945

178 -212458 -4.51692 -6.69052 12.96309 4.516917

179 -238560 -5.07185 -5.43064 14.55571 5.071855

180 -258957 -5.50548 -4.04088 15.80018 5.505485

181 -273159 -5.80744 -2.55446 16.66675 5.807436

182 -280829 -5.97049 -1.00695 17.13468 5.970486

183 -281781 -5.99073 0.564653 17.1928 5.990735

184 -275994 -5.8677 2.122746 16.8397 5.867699

185 -263606 -5.60432 3.630067 16.08383 5.604321

186 -244912 -5.2069 5.050565 14.94327 5.2069

187 -220362 -4.68494 6.350265 13.4453 4.684941

188 -190540 -4.05093 7.498081 11.62575 4.050929

189 -156161 -3.32003 8.466559 9.528143 3.320028

190 -118048 -2.50972 9.232536 7.202639 2.50972

191 -77110.3 -1.63938 9.777691 4.704864 1.639384

192 -34328.7 -0.72984 10.08899 2.09456 0.729838

193 9273.864 0.197165 10.15897 -0.56584 0.197165

194 52654.66 1.119451 9.985983 -3.21271 1.119451

195 94776.08 2.014963 9.574149 -5.78274 2.014963

196 134630.7 2.862281 8.933322 -8.21445 2.862281

197 171265.2 3.64114 8.078831 -10.4497 3.64114

198 203803.4 4.332911 7.031111 -12.435 4.332911

199 231467.1 4.921048 5.815223 -14.1229 4.921048

200 253594.7 5.391485 4.460248 -15.473 5.391485

201 269656.8 5.732969 2.998594 -16.453 5.732969

202 279269.3 5.937334 1.46522 -17.0395 5.937334

203 282202.3 5.999691 -0.1032 -17.2185 5.999691

204 278385.7 5.918548 -1.66915 -16.9856 5.918548

205 267910.7 5.695847 -3.19518 -16.3465 5.695847

206 251027.9 5.336914 -4.64478 -15.3164 5.336914

207 228141 4.850334 -5.9833 -13.92 4.850334

208 199797.5 4.247744 -7.1787 -12.1906 4.247744

209 166675.3 3.543558 -8.20241 -10.1697 3.543558

210 129566.6 2.754618 -9.02994 -7.90547 2.754618

211 89358.98 1.899794 -9.64149 -5.45221 1.899794

212 47014.07 0.999531 -10.0224 -2.86855 0.999531

213 3544.687 0.075361 -10.1637 -0.21628 0.075361

214 -40009.5 -0.85061 -10.0618 2.441168 0.850611

215 -82606.7 -1.75624 -9.71928 5.040226 1.756239

216 -123228 -2.61986 -9.1443 7.518734 2.619861

217 -160902 -3.42082 -8.35061 9.81741 3.420822

218 -194728 -4.13996 -7.35719 11.88128 4.139965

219 -223896 -4.76009 -6.1878 13.66097 4.760088

220 -247709 -5.26636 -4.87042 15.11392 5.266361

221 -265598 -5.64667 -3.43654 16.20538 5.646675

222 -277134 -5.89193 -1.92047 16.90924 5.891932

223 -282041 -5.99627 -0.35847 17.20868 5.996268

224 -280203 -5.95719 1.21211 17.09652 5.957186

225 -271663 -5.77562 2.753696 16.57544 5.775621

226 -256625 -5.45592 4.22942 15.65792 5.455917

227 -235450 -5.00572 5.603986 14.3659 5.005719

228 -208643 -4.4358 6.844517 12.73028 4.435795

229 -176845 -3.75978 7.921342 10.79018 3.759777

230 -140818 -2.99383 8.808707 8.591997 2.993834

231 -101423 -2.15628 9.485387 6.188316 2.156284

232 -59602.4 -1.26716 9.935197 3.636624 1.267162

233 -16355.9 -0.34773 10.14738 0.997952 0.347731

234 27281.74 0.580017 10.11686 -1.66459 0.580017

235 70266.89 1.493891 9.844366 -4.28732 1.493891

236 111571.4 2.372036 9.336418 -6.8075 2.372036

237 150207.4 3.193446 8.605163 -9.16486 3.193446

238 185250.7 3.938476 7.668092 -11.303 3.938476

239 215863.3 4.589307 6.547617 -13.1708 4.589307

240 241312.9 5.130372 5.270538 -14.7236 5.130372

241 260990.8 5.548729 3.867399 -15.9243 5.548729

242 274426.4 5.834373 2.37176 -16.7441 5.834373

243 281298.4 5.980473 0.819395 -17.1633 5.980473

244 281442.3 5.983532 -0.75257 -17.1721 5.983532

Fig- U vs Loop

0
0 50 100 150 200 250 300
-2

-4

-6

-8
Fig- U Error
8

0
0 50 100 150 200 250 300
-2

-4

-6

-8