O‘zbekiston respublikasi oliy va o‘rta maxsus ta’lim vazirligi m. M. Mirsaidov, T. M. Sobirjonov nazariy mexanika

> 
 sol:=fsolve(CHastot tenglamasi,omg=0..infinity); 
sol[2]/sol[1];#CHastotlarning oʻzaro bog’liqligi
 
:= 
sol
,
.9776197747 2.932988298
 
3.000131926
Tebranishning differentsial tenglamasi 
> 
eq1:=a11*diff(x(t),t$2)+a12*diff(s(t),t$2)=-c*x(t):
> 
eq2:=a12*diff(x(t),t$2)+a22*diff(s(t),t$2)=-c*s(t): 
Boshlang’ich shartlari berilgan sistemani sonli yechish. Yuklarning harakati
ularning tinch holatidan boshlanadi. 
> 
r:=dsolve({eq1,eq2,s(0)=0.1,x(0)=0.1,D(s)(0)=0,D(x)(0)=0}, 
{s(t),x(t)},type=numeric, 
output=listprocedure):
 
> 
with(plots):
 
> 
T0:=6.4:#Vaqt oraliqlarini tebranish davriga karrali qilib tanlaymiz.
 
>
s1:=odeplot(r,[t,x(t)],0..T0,numpoints=2.21,labels=[t,x],color=black,thickness
=2):
 

---

313 
> 
s2:=odeplot(r,[t,s(t)],0..T0,numpoints=2.21,labels=[t,x],color=red):
 
> 
display(s1,s2);
 
Warning, the name changecoords has been redefined 
3.23-rasm 
> 
rs := eval(s(t),r):
 
> 
rx := eval(x(t),r):
 
> 
n:=200: # animatsiya uchun nuqtalar (kadrlar) soni 
> 
k:=2000:# Rasm uchun masshtab koeffitsienti
> 
fd := fopen("C:\\DOS\\):# sfla.txt faylini ochish", WRITE fayl.
fprintf(fd, "ms=[ %3.0f",k*rs(0)):
 
> 
for i to n do
 
> 
tm:=T0/n*i:
 
>