BJD_Diplom4.mw

restart; 1 

`:=`(d, 0.1e-1) 

0.1e-1 (1)
 

 

`:=`(P, 14) 

14 (2)
 

 

`:=`(t0, .7) 

.7 (3)
 

 

`:=`(l, 2.0) 

2.0 (4)
 

 

`:=`(l:, 60) 

60 (5)
 

 

`:=`(l2, 2) 

2 (6)
 

 

(7)
 

`:=`(U;, .22) 

.22 (8)
 

 

 

`:=`(psi, 1.3) 

1.3 (9)
 

 

`:=`(b, 0.2e-1) 

0.2e-1 (10)
 

 

`:=`(P87<, 40) 

40 (11)
 



 

`:=`(I7, `+`(`*`(`/`(1, 350), `*`(U;, `*`(`+`(`*`(35, `*`(l:)), l2)))))) 

1.321257143 (12)
 

 

if `<=`(U;, 380) then `:=`(R7, 4) else `:=`(R7, `+`(`/`(`*`(125), `*`(I7)))) end if 

4 (13)
 

evalf(R7) 

4. (14)
 

 

 

`:=`(Pp, `*`(P87<, `*`(psi))) 

52.0 (15)
 

 

`:=`(R>2A, `*`(Pp, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(l))))), `*`(`+`(ln(`+`(`/`(`*`(2, `*`(l)), `*`(d)))), `*`(`/`(1, 2), `*`(ln(`/`(`*`(`+`(`*`(4, `+`(t0, `*`(`/`(1, 2), `*`(l)))), l)), `*`(`+`(`*`(4, `... 

`+`(`/`(`*`(81.82892180), `*`(Pi))) (16)
 

evalf(R>2A) 

26.04695478 (17)
 

 

`:=`(n?@81, round(`/`(`*`(R>2A), `*`(R7)))) 

7 (18)
 

evalf(n?@81); 1 

7. (19)
 

 

`:=`(a, `/`(`*`(P), `*`(n?@81))) 

2 (20)
 

evalf(a) 

2. (21)
 



 

 

evalf(`/`(`*`(a), `*`(l))) 

1.000000000 (22)
 

`:=`(n87, .65) 

.65 (23)
 



 

`:=`(n7, round(`/`(`*`(R>2A), `*`(n87, `*`(R7))))) 

10 (24)
 

evalf(n7); 1 

10. (25)
 


 

if `>`(n7, 20) then `:=`(l1, P) else `:=`(l1, `*`(a, `*`(`+`(n7, `-`(1))))) end if 

18 (26)
 

evalf(l1) 

18. (27)
 

 

 

`:=`(R?>;, `*`(Pp, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(l1))))), `*`(ln(`+`(`/`(`*`(2, `*`(`^`(l1, 2))), `*`(b, `*`(t0))))))))) 

`+`(`/`(`*`(15.51707250), `*`(Pi))) (28)
 

evalf(R?>;); 1 

4.939237580 (29)
 



 

`:=`(a1, `/`(`*`(P), `*`(n7))) 

`/`(7, 5) (30)
 

evalf(a1); 1 

1.400000000 (31)
 

evalf(round(`/`(`*`(a1), `*`(l)))) 

1. (32)
 

 

`:=`(n8?, .62) 

.62 (33)
 




 

`:=`(R?>;1, `/`(`*`(R?>;), `*`(n8?))) 

`+`(`/`(`*`(25.02753629), `*`(Pi))) (34)
 

evalf(R?>;1) 

7.966512226 (35)
 

 

`:=`(R?>;!, `/`(`*`(R?>;1, `*`(R7)), `*`(`+`(R?>;1, `-`(R7))))) 

`+`(`/`(`*`(100.1101452), `*`(Pi, `*`(`+`(`/`(`*`(25.02753629), `*`(Pi)), `-`(4)))))) (36)
 

evalf(R?>;!) 

8.033770502 (37)
 

 

`:=`(n72, `/`(`*`(R>2A), `*`(R?>;!, `*`(n87)))) 

`+`(`/`(`*`(31.47266222), `*`(Pi)), `-`(5.030085560)) (38)
 

evalf(R>2A) 

26.04695478 (39)
 

evalf(R?>;!) 

8.033770502 (40)
 

n87 

.65 (41)
 

 

trunc(evalf(n72)) 

4 (42)