Advertisement
Guest User

maple1

a guest
Oct 25th, 2016
60
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 14.02 KB | None | 0 0
  1. r1 := x*exp(x)-2 = 0;
  2. x exp(x) - 2 = 0
  3. plot(lhs(r1), x = 0 .. 1);
  4.  
  5. r2 := x = ln(2/x);
  6. /2\
  7. x = ln|-|
  8. \x/
  9. plot([abs(diff(rhs(r2), x)), 1], x = 0 .. 1);
  10.  
  11. eps := 10^(-8);
  12. 1
  13. ---------
  14. 100000000
  15. xs := 1.0;
  16. 1.0
  17. xn := 2/exp(xs);
  18. 0.7357588824
  19. i := 0;
  20. 0
  21. while abs(xn-xs) > eps do i := i+1; xs := xn; xn := 2/exp(xs) end do;
  22. 1
  23. 0.7357588824
  24. 0.9582834174
  25. 2
  26. 0.9582834174
  27. 0.7671014354
  28. 3
  29. 0.7671014354
  30. 0.9287141770
  31. 4
  32. 0.9287141770
  33. 0.7901227258
  34. 5
  35. 0.7901227258
  36. 0.9075782006
  37. 6
  38. 0.9075782006
  39. 0.8070004768
  40. 7
  41. 0.8070004768
  42. 0.8923888630
  43. 8
  44. 0.8923888630
  45. 0.8193518466
  46. 9
  47. 0.8193518466
  48. 0.8814344286
  49. 10
  50. 0.8814344286
  51. 0.8283767236
  52. 11
  53. 0.8283767236
  54. 0.8735153792
  55. 12
  56. 0.8735153792
  57. 0.8349627228
  58. 13
  59. 0.8349627228
  60. 0.8677813106
  61. 14
  62. 0.8677813106
  63. 0.8397642092
  64. 15
  65. 0.8397642092
  66. 0.8636246576
  67. 16
  68. 0.8636246576
  69. 0.8432620820
  70. 17
  71. 0.8432620820
  72. 0.8606090854
  73. 18
  74. 0.8606090854
  75. 0.8458088378
  76. 19
  77. 0.8458088378
  78. 0.8584201128
  79. 20
  80. 0.8584201128
  81. 0.8476623182
  82. 21
  83. 0.8476623182
  84. 0.8568305216
  85. 22
  86. 0.8568305216
  87. 0.8490108262
  88. 23
  89. 0.8490108262
  90. 0.8556758574
  91. 24
  92. 0.8556758574
  93. 0.8499917146
  94. 25
  95. 0.8499917146
  96. 0.8548369466
  97. 26
  98. 0.8548369466
  99. 0.8507050812
  100. 27
  101. 0.8507050812
  102. 0.8542273518
  103. 28
  104. 0.8542273518
  105. 0.8512238248
  106. 29
  107. 0.8512238248
  108. 0.8537843416
  109. 30
  110. 0.8537843416
  111. 0.8516010092
  112. 31
  113. 0.8516010092
  114. 0.8534623684
  115. 32
  116. 0.8534623684
  117. 0.8518752460
  118. 33
  119. 0.8518752460
  120. 0.8532283496
  121. 34
  122. 0.8532283496
  123. 0.8520746242
  124. 35
  125. 0.8520746242
  126. 0.8530582516
  127. 36
  128. 0.8530582516
  129. 0.8522195726
  130. 37
  131. 0.8522195726
  132. 0.8529346112
  133. 38
  134. 0.8529346112
  135. 0.8523249480
  136. 39
  137. 0.8523249480
  138. 0.8528447376
  139. 40
  140. 0.8528447376
  141. 0.8524015528
  142. 41
  143. 0.8524015528
  144. 0.8527794078
  145. 42
  146. 0.8527794078
  147. 0.8524572418
  148. 43
  149. 0.8524572418
  150. 0.8527319190
  151. 44
  152. 0.8527319190
  153. 0.8524977252
  154. 45
  155. 0.8524977252
  156. 0.8526973980
  157. 46
  158. 0.8526973980
  159. 0.8525271546
  160. 47
  161. 0.8525271546
  162. 0.8526723040
  163. 48
  164. 0.8526723040
  165. 0.8525485482
  166. 49
  167. 0.8525485482
  168. 0.8526540626
  169. 50
  170. 0.8526540626
  171. 0.8525640998
  172. 51
  173. 0.8525640998
  174. 0.8526408024
  175. 52
  176. 0.8526408024
  177. 0.8525754054
  178. 53
  179. 0.8525754054
  180. 0.8526311628
  181. 54
  182. 0.8526311628
  183. 0.8525836236
  184. 55
  185. 0.8525836236
  186. 0.8526241560
  187. 56
  188. 0.8526241560
  189. 0.8525895976
  190. 57
  191. 0.8525895976
  192. 0.8526190624
  193. 58
  194. 0.8526190624
  195. 0.8525939406
  196. 59
  197. 0.8525939406
  198. 0.8526153594
  199. 60
  200. 0.8526153594
  201. 0.8525970976
  202. 61
  203. 0.8525970976
  204. 0.8526126678
  205. 62
  206. 0.8526126678
  207. 0.8525993924
  208. 63
  209. 0.8525993924
  210. 0.8526107112
  211. 64
  212. 0.8526107112
  213. 0.8526010608
  214. 65
  215. 0.8526010608
  216. 0.8526092886
  217. 66
  218. 0.8526092886
  219. 0.8526022736
  220. 67
  221. 0.8526022736
  222. 0.8526082546
  223. 68
  224. 0.8526082546
  225. 0.8526031550
  226. 69
  227. 0.8526031550
  228. 0.8526075032
  229. 70
  230. 0.8526075032
  231. 0.8526037958
  232. 71
  233. 0.8526037958
  234. 0.8526069566
  235. 72
  236. 0.8526069566
  237. 0.8526042618
  238. 73
  239. 0.8526042618
  240. 0.8526065594
  241. 74
  242. 0.8526065594
  243. 0.8526046006
  244. 75
  245. 0.8526046006
  246. 0.8526062708
  247. 76
  248. 0.8526062708
  249. 0.8526048466
  250. 77
  251. 0.8526048466
  252. 0.8526060610
  253. 78
  254. 0.8526060610
  255. 0.8526050254
  256. 79
  257. 0.8526050254
  258. 0.8526059084
  259. 80
  260. 0.8526059084
  261. 0.8526051556
  262. 81
  263. 0.8526051556
  264. 0.8526057974
  265. 82
  266. 0.8526057974
  267. 0.8526052502
  268. 83
  269. 0.8526052502
  270. 0.8526057168
  271. 84
  272. 0.8526057168
  273. 0.8526053188
  274. 85
  275. 0.8526053188
  276. 0.8526056582
  277. 86
  278. 0.8526056582
  279. 0.8526053690
  280. 87
  281. 0.8526053690
  282. 0.8526056154
  283. 88
  284. 0.8526056154
  285. 0.8526054054
  286. 89
  287. 0.8526054054
  288. 0.8526055844
  289. 90
  290. 0.8526055844
  291. 0.8526054318
  292. 91
  293. 0.8526054318
  294. 0.8526055620
  295. 92
  296. 0.8526055620
  297. 0.8526054508
  298. 93
  299. 0.8526054508
  300. 0.8526055456
  301. 94
  302. 0.8526055456
  303. 0.8526054650
  304. 95
  305. 0.8526054650
  306. 0.8526055336
  307. 96
  308. 0.8526055336
  309. 0.8526054750
  310. 97
  311. 0.8526054750
  312. 0.8526055248
  313. 98
  314. 0.8526055248
  315. 0.8526054828
  316. 99
  317. 0.8526054828
  318. 0.8526055184
  319. 100
  320. 0.8526055184
  321. 0.8526054882
  322. 101
  323. 0.8526054882
  324. 0.8526055136
  325. 102
  326. 0.8526055136
  327. 0.8526054922
  328. 103
  329. 0.8526054922
  330. 0.8526055104
  331. 104
  332. 0.8526055104
  333. 0.8526054948
  334. 105
  335. 0.8526054948
  336. 0.8526055082
  337. 106
  338. 0.8526055082
  339. 0.8526054968
  340. 107
  341. 0.8526054968
  342. 0.8526055064
  343.  
  344. r1;
  345. x exp(x) - 2 = 0
  346. r2;
  347. /2\
  348. x = ln|-|
  349. \x/
  350. r3 := (r2-x*(diff(r2, x)))/(1-(diff(r2, x)));
  351. /2\
  352. ln|-| + 1
  353. \x/
  354. 0 = ---------
  355. 1
  356. - + 1
  357. x
  358. r4 := x = rhs(r3);
  359. /2\
  360. ln|-| + 1
  361. \x/
  362. x = ---------
  363. 1
  364. - + 1
  365. x
  366. plot([abs(diff(rhs(r4), x)), 1], x = 0 .. 1);
  367.  
  368. xs := 1.0;
  369. 1.0
  370. xn := (ln(2/xs)+1)/(1/xs+1);
  371. 0.8465735905
  372.  
  373. i := 0;
  374. 0
  375. while abs(xn-xs) > eps do i := i+1; xs := xn; xn := (ln(2/xs)+1)/(1/xs+1) end do;
  376. 1
  377. 0.8465735905
  378. 0.8525939199
  379. 2
  380. 0.8525939199
  381. 0.8526055019
  382. 3
  383. 0.8526055019
  384. 0.8526055018
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement