Lab 8.

wxmx 1 / 19

LABORATORIO 8
JORGE SERRANO
8-934-297
Figure 1:

Figure 2:

m1:1;m2:1;m3:1;
1
1
1
k1:100;k2:100;k3:100;k4:100;
100
100
100
100
c1:1;c2:1;c3:1;c4:1;
1
1
1
1
Fo:10;
10
w:1;
1
Matriz
M: matrix(
[m1,0,0],
[0,m2,0],
[0,0,m3]
);
Lab 8.wxmx 2 / 19

1 0 0
0 1 0
0 0 1
Matriz de amortiguamiento
C: matrix(
[c1+c2,−c2,0],
[−c2,c2+c3,−c3],
[0,−c3,c3+c4]
);
2 −1 0
−1 2 −1
0 −1 2
Matriz de rigidez
K: matrix(
[k1+k2,−k2,0],
[−k2,k2+k3,−k3],
[0,−k3,k3+k4]
);
200 − 100 0
− 100 200 − 100
0 − 100 200
Polinomio Caracteristico, Ecuacion caracteristica
PC:determinant(−λ·M+K);
2
( 200 −λ ) −10000 ( 200 −λ ) −10000 ( 200 −λ )
expand(%);
3 2
−λ +600 λ −100000 λ +4000000
load(assign);
C:/maxima−5.45.0/share/maxima/5.45.0/share/assign.mac
allroots(PC);
[λ =58.5786437626905 , λ =200.0 , λ =
341.4213562373096]
Lab 8.wxmx 3 / 19

assign(%);
[ 58.5786437626905 , 200.0 ,341.4213562373096 ]
Frecuencias Naturales, Frecuencia de resonancia ,
Frecuencia de oscilacion
ω:sqrt(%);
[ 7.653668647301796 , 14.14213562373095,
18.47759065022574]
Valores Propios
λ1:%o21[1];
58.5786437626905
λ2:%o21[2];
200.0
λ3:%o21[3];
341.4213562373096
Vectores propios (Modo de vibración)
V1:−λ1·M+K;
141.4213562373095 − 100 0
− 100 141.4213562373095 − 100
0 − 100 141.4213562373095
E1:echelon(%);
rat: replaced 141.4213562373095 by 331531257/2344280 = 141.421356
rat: replaced 141.4213562373095 by 331531257/2344280 = 141.421356
rat: replaced 141.4213562373095 by 331531257/2344280 = 141.421356
234428000
1 − 0
331531257

77720209515996000
0 1 −
54956487184000049
0 0 1
v1: matrix(
[1],
[−1/E1[1,2]],
[−(−1/E1[1,2])/E1[2,3]]
);
Lab 8.wxmx 4 / 19

1
331531257
234428000

54956487184000049
54956487184000000
float(%);
1.0
1.414213562373095
1.000000000000001
V2:−λ2·M+K;
−14
2.842170943040401 10 − 100 0
−14
− 100 2.842170943040401 10 − 100
−14
0 − 100 2.842170943040401 10
E2:echelon(%);
rat: replaced 2.842170943040401e−14 by 1/35184372088832 = 2.84217
rat: replaced 2.842170943040401e−14 by 1/35184372088832 = 2.84217
rat: replaced 2.842170943040401e−14 by 1/35184372088832 = 2.84217
1 − 3518437208883200 0
3518437208883200
0 1
12379400392853802748991242239999
0 0 1
v2:matrix(
[1],
[−1/E2[1,2]],
[−(−1/E2[1,2])/E2[2,3]]
);
1
1
3518437208883200

12379400392853802748991242239999
−
12379400392853802748991242240000
Lab 8.wxmx 5 / 19

float(%);
1.0
−16
2.842170943040401 10
− 1.0
V3:−λ3·M+K;
− 141.4213562373096 − 100 0
− 100 − 141.4213562373096 − 100
0 − 100 − 141.4213562373096
E3:echelon(%);
rat: replaced −141.4213562373096 by −233957168/1654327 = −141.421
rat: replaced −141.4213562373096 by −233957168/1654327 = −141.421
rat: replaced −141.4213562373096 by −233957168/1654327 = −141.421
41358175
1 0
58489292

2419010374162100
0 1
1710498639330639
0 0 1
v3:matrix(
[1],
[−1/E3[1,2]],
[−(−1/E3[1,2])/E3[2,3]]
);
1
58489292
−
41358175

1710498639330639
1710498639330625
float(%);
1.0
− 1.414213562373098
1.000000000000008
Lab 8.wxmx 6 / 19

Generar matriz ψ cambio de bases ortogonales

(Vectores propios)
fi: matrix(
[v1[1,1],v2[1,1],v3[1,1]],
[v1[2,1],v2[2,1],v3[2,1]],
[v1[3,1],v2[3,1],v3[3,1]]
);
1 1 1
331531257 1 58489292
−
234428000 3518437208883200 41358175

54956487184000049 12379400392853802748991242239999 1710498639330639

−
54956487184000000 12379400392853802748991242240000 1710498639330625
ψ:float(%);
1.0 1.0 1.0
−16
1.414213562373095 2.842170943040401 10 − 1.414213562373098
1.000000000000001 − 1.0 1.000000000000008
Matriz M*...Ms(diagonal)...Ms=ψT*M*ψ
Ms:transpose(ψ).M.ψ;
−16 −15
4.000000000000003 − 4.440892098500626 10 4.218847493575595 10
−16 −15
− 4.440892098500626 10 2.0 − 8.659739592076221 10
−15 −15
4.218847493575595 10 − 8.659739592076221 10 4.000000000000025
Cs:transpose(ψ).C.ψ;
−15 −15
2.343145750507621 − 1.77635683940025 10 4.440892098500626 10
−15 −14
− 1.554312234475219 10 4.0 − 1.77635683940025 10
−15 −14
3.885780586188048 10 − 1.798561299892754 10 13.65685424949246
Ks:transpose(ψ).K.ψ;
Lab 8.wxmx 7 / 19

−13 −13
234.3145750507622 − 1.4210854715202 10 4.547473508864641 10
−13 −12
− 1.49213974509621 10 400.0 − 1.762145984685048 10
−13 −12
3.694822225952521 10 − 1.790567694115452 10 1365.685424949246
Generar la matriz fuerza (Base ortogonal)
F: matrix(
[Fo·cos(w·t)],
[0],
[0]
);
10 cos ( t )
0
0
Fs:transpose(ψ).F;
10.0 cos ( t )

10.0 cos ( t )

10.0 cos ( t )
Generar las ecuaciones diferenciales que gobiernan el
sistema dinamico
ED1:Ms[1,1]·'diff(z1,t,2)+Cs[1,1]·'diff(z1,t)+Ks[1,1]·z1=Fs[1,1];
2
d
4.000000000000003 z1 +
2
dt
d
2.343145750507621 z1 +234.3145750507622 z1 =
dt
10.0 cos ( t )
ode2(ED1,z1,t)$
rat: replaced −10.0 by −10/1 = −10.0
rat: replaced 234.3145750507622 by 137324743/586070 = 234.3145750
rat: replaced 2.343145750507621 by 81108829/34615358 = 2.34314575
rat: replaced 4.000000000000003 by 4/1 = 4.0
subst([%k1=A1,%k2=B1],float(%));
Lab 8.wxmx 8 / 19

z1 = A1 sin ( 7.648062324867885 t ) +B1

cos ( 7.648062324867885 t ) /
81108829 t
276922864 −31
2.718281828459045 +3.6640851369145 10 (
27
1.205440987002163 10 sin ( t ) +1.184862822169846
29
10 cos ( t ) )
expand(%);
A1 sin ( 7.648062324867885 t )
z1 = +
81108829 t
276922864
2.718281828459045
B1 cos ( 7.648062324867885 t ) −4
+4.416838403902171 10
81108829 t
276922864
2.718281828459045
sin ( t ) +0.043414382559951 cos ( t )
assign(%);
A1 sin ( 7.648062324867885 t )
+
81108829 t
276922864
2.718281828459045
B1 cos ( 7.648062324867885 t ) −4
+4.416838403902171 10
81108829 t
276922864
2.718281828459045
sin ( t ) +0.043414382559951 cos ( t )
ED2:Ms[2,2]·'diff(z2,t,2)+Cs[2,2]·'diff(z2,t)+Ks[2,2]·z2=Fs[2,1];
2
d d
2.0 z2 +4.0 z2 +400.0 z2 = 10.0 cos ( t )
2 dt
dt
ode2(ED2,z2,t);
rat: replaced −10.0 by −10/1 = −10.0
rat: replaced 400.0 by 400/1 = 400.0
rat: replaced 4.0 by 4/1 = 4.0
rat: replaced 2.0 by 2/1 = 2.0
p p
Lab 8.wxmx 9 / 19

−t p p
z2 = %e ( %k1 sin ( 199 t ) +%k2 cos ( 199 t ) ) +
2 sin ( t ) + 199 cos ( t )
7921
subst([%k1=A2,%k2=B2],float(%));
z2 = A2 sin ( 14.10673597966588 t ) +B2
t
cos ( 14.10673597966588 t ) / 2.718281828459045 +
−4
1.262466860244918 10 ( 2.0 sin( t ) +199.0 cos ( t ) )
assign(%);
A2 sin ( 14.10673597966588 t ) +B2
t
cos ( 14.10673597966588 t ) / 2.718281828459045 +
−4
1.262466860244918 10 ( 2.0 sin( t ) +199.0 cos ( t ) )
ED3:Ms[3,3]·'diff(z3,t,2)+Cs[3,3]·'diff(z3,t)+Ks[3,3]·z3=Fs[3,1];
2
d
4.000000000000025 z3 +
2
dt
d
13.65685424949246 z3 +1365.685424949246 z3 =
dt
10.0 cos ( t )
ode2(ED3,z3,t);
rat: replaced −10.0 by −10/1 = −10.0
rat: replaced 1365.685424949246 by 516644259/378304 = 1365.685424
rat: replaced 13.65685424949246 by 92085669/6742817 = 13.65685424
rat: replaced 4.000000000000025 by 160842843834661/4021071095866
3702830219594297871885 t
−
2169067723473394760074
z3 = %e (%k1 sin ( 3 sqrt(
65120746454367479551748683545074107581758225055
)t/ )+%k2 cos ( 3 sqrt(
65120746454367479551748683545074107581758225055
)t/ )) +
21753670417110085422730104239484890653532992512000
sin ( t ) +
Lab 8.wxmx 10 / 19

2168995539161536060970748124122307877979215022441600
cos ( t ) /
295378669916247561928299886963248016117730230102900441
assign(%);
3702830219594297871885 t
−
2169067723473394760074
%e (%k1 sin ( 3 sqrt(
65120746454367479551748683545074107581758225055
)t/ )+%k2 cos ( 3 sqrt(
65120746454367479551748683545074107581758225055
)t/ )) +
21753670417110085422730104239484890653532992512000
sin ( t ) +
2168995539161536060970748124122307877979215022441600
cos ( t ) /
295378669916247561928299886963248016117730230102900441
z3:subst([%k1=A3,%k2=B3],float(%));
A3 sin ( 18.39856360357886 t ) +B3
cos ( 18.39856360357886 t ) /
3702830219594297871885 t
2169067723473394760074
2.718281828459045 +
−54 49
3.38548480932473 10 (2.175367041711009 10
51
sin ( t ) +2.168995539161536 10 cos( t ) )
expand(%);
A3 sin ( 18.39856360357886 t )
+
3702830219594297871885 t
2169067723473394760074
2.718281828459045
B3 cos ( 18.39856360357886 t )
+
3702830219594297871885 t
2169067723473394760074
2.718281828459045
−5
7.364672074418298 10 sin ( t ) +
0.007343101449324484 cos ( t )
Lab 8.wxmx 11 / 19

Z: matrix(
[z1],
[z2],
[z3]
);
A1 sin ( 7.648062324867885 t ) B1 cos ( 7.648062324867885 t )
+ + 4.4168384
81108829 t 81108829 t
276922864 276922864
2.718281828459045 2.718281828459045

A2 sin ( 14.10673597966588 t ) + B2 cos ( 14.10673597966588 t )

+ 1.2624
t
2.718281828459045

A3 sin ( 18.39856360357886 t ) + B3 cos ( 18.39856360357886 t ) −54

+ 3.38548480932473 10 2
3702830219594297871885 t
2169067723473394760074
2.718281828459045
X:ψ.Z;
A3 sin ( 18.39856360357886 t ) + B3 cos ( 18.398563603578
1.0 3702830219594297871885 t
2169067723473394760074
2.718281828459045

A3 sin ( 18.39856360357886 t ) + B3 cos ( 18.39856360357886 t )

+ 3.38
− 1.414213562373098 3702830219594297871885 t
2169067723473394760074
2.718281828459045

A3 sin ( 18.39856360357886 t ) + B3 cos ( 18.398563603578

1.000000000000008 3702830219594297871885 t
2169067723473394760074
2.718281828459045

x1:X[1,1];
1.0 ( A3 sin ( 18.39856360357886 t ) +B3
cos ( 18.39856360357886 t ) /
3702830219594297871885 t
2169067723473394760074
2.718281828459045 +
−54 49
3.38548480932473 10 (2.175367041711009 10
51
sin ( t ) +2.168995539161536 10 cos( t ) )) +1.0 ( A2
sin ( 14.10673597966588 t ) +B2
Lab 8.wxmx 12 / 19

t
cos ( 14.10673597966588 t ) / 2.718281828459045 +
−4
1.262466860244918 10 ( 2.0 sin( t ) +199.0 cos ( t ) ) )+1.0
A1 sin ( 7.648062324867885 t )
( +
81108829 t
276922864
2.718281828459045
B1 cos ( 7.648062324867885 t ) −4
+4.416838403902171 10
81108829 t
276922864
2.718281828459045
sin ( t ) +0.043414382559951 cos ( t ) )
x2:X[2,1]$;
x3:X[3,1];
1.000000000000008 ( A3
sin ( 18.39856360357886 t ) +B3
cos ( 18.39856360357886 t ) /
3702830219594297871885 t
2169067723473394760074
2.718281828459045 +
−54 49
3.38548480932473 10 (2.175367041711009 10
51
sin ( t ) +2.168995539161536 10 cos( t ) )) −1.0 ( A2
sin ( 14.10673597966588 t ) +B2
t
cos ( 14.10673597966588 t ) / 2.718281828459045 +
−4
1.262466860244918 10 ( 2.0 sin( t ) +199.0 cos ( t ) ) )+
A1 sin ( 7.648062324867885 t )
1.000000000000001 ( +
81108829 t
276922864
2.718281828459045
B1 cos ( 7.648062324867885 t ) −4
+4.416838403902171 10
81108829 t
276922864
2.718281828459045
sin ( t ) +0.043414382559951 cos ( t ) )
Aplicamos condiciones iniciales de posición
Lab 8.wxmx 13 / 19

x1o:0;x2o:0;x3o:0;
0
0
0
eq1:subst([t=0],x1)=x1o;
1.0 ( 1.0 B3+ 0.007343101449324484 ) +1.0
( 1.0 B2+ 0.02512309051887388 ) +1.0
( 1.0 B1+ 0.043414382559951 ) =0
expand(%);
1.0 B3+ 1.0 B2+ 1.0 B1+ 0.07588057452814936 =0
eq2:subst([t=0],x2)=x2o;
−1.414213562373098
( 1.0 B3+ 0.007343101449324484 ) +2.842170943040401
−16
10 ( 1.0 B2+ 0.02512309051887388 ) +
1.414213562373095 ( 1.0 B1+ 0.043414382559951 ) =0
expand(%);
−1.414213562373098 B3+ 2.842170943040401
−16
10 B2+ 1.414213562373095 B1+
0.05101249495882046 =0
eq3:subst([t=0],x3)=x3o;
1.000000000000008
( 1.0 B3+ 0.007343101449324484 ) −1.0
( 1.0 B2+ 0.02512309051887388 ) +1.000000000000001
( 1.0 B1+ 0.043414382559951 ) =0
expand(%);
1.000000000000008 B3−1.0 B2+
1.000000000000001 B1+ 0.02563439349040171 =0
linsolve([eq1,eq2,eq3],[B1,B2,B3]);
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 0.043414382559951 by 7606051/175196572 = 0.04341438
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 0.02512309051887388 by 199/7921 = 0.025123090518873
Lab 8.wxmx 14 / 19

rat: replaced 1.0 by 1/1 = 1.0

rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 0.007343101449324484 by 1152319/156925382 = 0.00734
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 1.414213562373095 by 22619537/15994428 = 1.41421356
rat: replaced 0.043414382559951 by 7606051/175196572 = 0.04341438
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 2.842170943040401e−16 by 1/3518437208883200 = 2.842
rat: replaced 0.02512309051887388 by 199/7921 = 0.025123090518873
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced −1.414213562373098 by −22619537/15994428 = −1.41421
rat: replaced 0.007343101449324484 by 1152319/156925382 = 0.00734
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 1.000000000000001 by 1/1 = 1.0
rat: replaced 0.043414382559951 by 7606051/175196572 = 0.04341438
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced −1.0 by −1/1 = −1.0
rat: replaced 0.02512309051887388 by 199/7921 = 0.025123090518873
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 1.000000000000008 by 121718908847852/1217189088478
rat: replaced 0.007343101449324484 by 1152319/156925382 = 0.00734
rat: replaced 1.0 by 1/1 = 1.0
7606051 199
[B1= − ,B2= − ,B3= −
175196572 7921
1152319
]
156925382
float(%);
[B1= −0.043414382559951, B2= −
0.02512309051887388 , B3= −0.007343101449324495 ]
assign(%);
[−0.043414382559951, −0.02512309051887388 , −
0.007343101449324495 ]
Condiciones iniciales de velocidad
v1o:0;v2o:0;v3o:0;
0
0
Lab 8.wxmx 15 / 19

0
subst([t=0],diff(x1,t))=v1o;
1.0 (−1.707106781186547 B3+
−5
18.39856360357886 A3 +7.364672074418298 10 )+1.0
(−1.0 B2+ 14.10673597966588 A2 +2.524933720489837
−4
10 )+1.0 (−0.2928932188134527 B1+
−4
7.648062324867885 A1 +4.416838403902171 10 )=0
eq4:''%;
1.0 (18.39856360357886 A3 +
0.01260910499982679 ) +1.0 (14.10673597966588 A2 +
0.02537558389092286 ) +1.0 (7.648062324867885 A1 +
0.01315746209117289 ) =0
expand(%);
18.39856360357886 A3 +14.10673597966588 A2 +
7.648062324867885 A1 +0.05114215098192254 =0
subst([t=0],diff(x2,t))=v2o;
−1.414213562373098 ( −1.707106781186547 B3+
−5
18.39856360357886 A3 +7.364672074418298 10 )+
−16
2.842170943040401 10 (−1.0 B2+
−4
14.10673597966588 A2 +2.524933720489837 10 )+
1.414213562373095 ( −0.2928932188134527 B1+
−4
7.648062324867885 A1 +4.416838403902171 10 )=0
eq5:''%;
−1.414213562373098 (18.39856360357886 A3 +
−16
0.01260910499982679 ) +2.842170943040401 10 (
14.10673597966588 A2 +0.02537558389092286 ) +
1.414213562373095 (7.648062324867885 A1 +
0.01315746209117289 ) =0
expand(%);
−26.01949817636528 A3 +4.009375510254894
−15
10 A2 +10.81599346570287 A1 +7.754940356050933
Lab 8.wxmx 16 / 19

−4
10 =0
subst([t=0],diff(x3,t))=v3o;
1.000000000000008 ( −1.707106781186547 B3+
−5
18.39856360357886 A3 +7.364672074418298 10 )−1.0
(−1.0 B2+ 14.10673597966588 A2 +2.524933720489837
−4
10 )+1.000000000000001 ( −0.2928932188134527 B1+
−4
7.648062324867885 A1 +4.416838403902171 10 )=0
eq6:''%;
1.000000000000008 (18.39856360357886 A3 +
0.01260910499982679 ) −1.0 (14.10673597966588 A2 +
0.02537558389092286 ) +1.000000000000001 (
7.648062324867885 A1 +0.01315746209117289 ) =0
expand(%);
18.39856360357901 A3 −14.10673597966588 A2 +
−4
7.648062324867892 A1 +3.909832000769334 10 =0
linsolve([eq4,eq5,eq6],[A1,A2,A3]);
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 0.01315746209117289 by 1586607/120586097 = 0.013157
rat: replaced 7.648062324867885 by 64353565/8414362 = 7.648062324
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 0.02537558389092286 by 201/7921 = 0.025375583890922
rat: replaced 14.10673597966588 by 84559415/5994258 = 14.10673597
rat: replaced 1.0 by 1/1 = 1.0
rat: replaced 0.01260910499982679 by 6589884/522629005 = 0.012609
rat: replaced 18.39856360357886 by 49782815/2705799 = 18.39856360
rat: replaced 1.414213562373095 by 22619537/15994428 = 1.41421356
rat: replaced 0.01315746209117289 by 1586607/120586097 = 0.013157
rat: replaced 7.648062324867885 by 64353565/8414362 = 7.648062324
rat: replaced 2.842170943040401e−16 by 1/3518437208883200 = 2.842
rat: replaced 0.02537558389092286 by 201/7921 = 0.025375583890922
rat: replaced 14.10673597966588 by 84559415/5994258 = 14.10673597
rat: replaced −1.414213562373098 by −22619537/15994428 = −1.41421
rat: replaced 0.01260910499982679 by 6589884/522629005 = 0.012609
rat: replaced 18.39856360357886 by 49782815/2705799 = 18.39856360
Lab 8.wxmx 17 / 19

rat: replaced 1.000000000000001 by 1/1 = 1.0

rat: replaced 0.01315746209117289 by 1586607/120586097 = 0.013157
rat: replaced 7.648062324867885 by 64353565/8414362 = 7.648062324
rat: replaced −1.0 by −1/1 = −1.0
rat: replaced 0.02537558389092286 by 201/7921 = 0.025375583890922
rat: replaced 14.10673597966588 by 84559415/5994258 = 14.10673597
rat: replaced 1.000000000000008 by 121718908847852/1217189088478
rat: replaced 0.01260910499982679 by 6589884/522629005 = 0.012609
rat: replaced 18.39856360357886 by 49782815/2705799 = 18.39856360
13350285649734 1204845858
[A1 =− ,A2 =− ,
7760145231385805 669795126215
17830901537316
A3 =− ]
26017943069549075
float(%);
[A1 =−0.001720365437973886 , A2 =−
−4
0.001798827448638753 , A3 =−6.853309460187482 10
]
assign(%);
[−0.001720365437973886 , −
−4
0.001798827448638753 , −6.853309460187482 10 ]
x1t:''x1;
−4
1.0 ( −6.853309460187482 10
sin ( 18.39856360357886 t ) −0.007343101449324495
cos ( 18.39856360357886 t ) /
3702830219594297871885 t
2169067723473394760074
2.718281828459045 +
−54 49
3.38548480932473 10 (2.175367041711009 10
51
sin ( t ) +2.168995539161536 10 cos( t ) )) +1.0 ( −
0.001798827448638753 sin( 14.10673597966588 t ) −
0.02512309051887388 cos( 14.10673597966588 t ) /
t −4
2.718281828459045 +1.262466860244918 10
( 2.0 sin( t ) +199.0 cos( t ) ) )+1.0 (−
Lab 8.wxmx 18 / 19

0.001720365437973886 sin ( 7.648062324867885 t )

−
81108829 t
276922864
2.718281828459045
0.043414382559951 cos ( 7.648062324867885 t )
+
81108829 t
276922864
2.718281828459045
−4
4.416838403902171 10 sin ( t ) +0.043414382559951
cos ( t ) )
x2t:''x2;
−4
−1.414213562373098 ( −6.853309460187482 10
sin ( 18.39856360357886 t ) −0.007343101449324495
cos ( 18.39856360357886 t ) /
3702830219594297871885 t
2169067723473394760074
2.718281828459045 +
−54 49
3.38548480932473 10 (2.175367041711009 10
51
sin ( t ) +2.168995539161536 10 cos( t ) )) +
−16
2.842170943040401 10 ( −0.001798827448638753
sin ( 14.10673597966588 t ) −0.02512309051887388
t
cos ( 14.10673597966588 t ) / 2.718281828459045 +
−4
1.262466860244918 10 ( 2.0 sin( t ) +199.0 cos ( t ) ) )+
1.414213562373095 ( −
0.001720365437973886 sin ( 7.648062324867885 t )
−
81108829 t
276922864
2.718281828459045
0.043414382559951 cos ( 7.648062324867885 t )
+
81108829 t
276922864
2.718281828459045
−4
4.416838403902171 10 sin ( t ) +0.043414382559951
cos ( t ) )
Lab 8.wxmx 19 / 19

x3t:''x3;
−4
1.000000000000008 ( −6.853309460187482 10
sin ( 18.39856360357886 t ) −0.007343101449324495
cos ( 18.39856360357886 t ) /
3702830219594297871885 t
2169067723473394760074
2.718281828459045 +
−54 49
3.38548480932473 10 (2.175367041711009 10
51
sin ( t ) +2.168995539161536 10 cos( t ) )) −1.0 ( −
0.001798827448638753 sin( 14.10673597966588 t ) −
0.02512309051887388 cos( 14.10673597966588 t ) /
t −4
2.718281828459045 +1.262466860244918 10
( 2.0 sin( t ) +199.0 cos( t ) ) )+1.000000000000001 ( −
0.001720365437973886 sin ( 7.648062324867885 t )
−
81108829 t
276922864
2.718281828459045
0.043414382559951 cos ( 7.648062324867885 t )
+
81108829 t
276922864
2.718281828459045
−4
4.416838403902171 10 sin ( t ) +0.043414382559951
cos ( t ) )
wxplot2d([x1t,x2t,x3t],[t,0,40])$