Программа 36. Задача 50

>	restart; n:=4: # Число участков

>	m:=0,10,-18,15:# Моменты (кН)

>	L:=2,4,3,5:    # Длины участков(м)

>	dm:=0.09,0.18,0.25,0.32:#Диаметры (м)

>	G:=80e9:       # Модуль сдвига(Па)

>	M[0]:=-Ma: f[0]:=0: x[0]:=0:

>	for i to n do

>	Ji:=Pi*dm[i]^4/32:

>	Wi:=2*Ji/dm[i]:

>	M[i]:=M[i-1]+m[i]*1e3:

>	ffi:=M[i]*L[i]/G/Ji;

>	f[i]:=f[i-1]+ffi:

>	x[i]:=x[i-1]+L[i]:

>	tau[i]:=M[i]/Wi/1e6;# МПа

>

od:

>	Ma:=solve(f[n],Ma);#Момент на опоре (Нм)

>	plot([seq([x[i],f[i]],i=0..n)],title=`fi`);

>	plot([seq(op([[x[i],tau[i+1]],[x[i+1],tau[i+1]]]),i=0..n-1)],title=`tau`);

>	plot([seq(op([[x[i],M[i+1]],[x[i+1],M[i+1]]]),i=0..n-1)],title=`M`);

>	evalf(seq(tau[i],i=1..n));     # касательные напряжения (МПа)

>	evalf(seq(f[i]*1000,i=1..n-1));# углы*1000

>	evalf(seq(M[i],i=1..n));       # моменты (Нм)