MATH 5 - Differenial Equaions Sepember 15, 8 Projec 1, Fall 8 Due: Sepember 4, 8 Lab 1.3 - Logisics Populaion Models wih Harvesing For his projec we consider lab 1.3 of Differenial Equaions pages 146 o 147. Afer reading his maerial consruc a lab repor addressing each of he following quesions for cases 3, 4, 6, 8 of able 1.1 on page 147: 1 1. Given a logisics growh model wih consan harvesing, dp (1 d = kp p ) a, k, N, a R +. (1) N (a) Consruc a lis of variables and parameers associaed wih (3) and describe he meaning of each. (b) Analicall solve (3) using he mehods discussed in secion 1. of he ex. (c) Discuss qualiaive behavior of he soluions o (3) hrough he equaion s: i. Equilibrium Poins ii. Phase Line (d) Using Euler s mehod and a slope field diagram address he following quesion: For a = a 1, wha will happen o he fish populaion for various iniial condiions?. Given a logisics growh model wih periodic harvesing, dp (1 d = kp p ) a(1 + sin(b)), k, N, a, b R +. () N (a) In his case wha do he parameers a and b represen? (b) Is i possible o solve (3) using he mehods discussed in secion 1. of he ex? (c) Using Euler s mehod and a slope field diagram address he following quesions: i. For a = a 1 and b = 1 wha will happen o he fish populaion for various iniial condiions? ii. For a = a and b = 1 wha will happen o he fish populaion for various iniial condiions? Explain wh here are no equilibrium poins and hus no phase line for his problem. 3. Summar and conclusions. In a shor essa forma summarize our resuls from he previous quesions. Compare and conras each of he wo models. Be sure o jusif our conclusions b referencing our previous summar and analsis. 1 Your repor should be well organized and clearl presened. If seps are unclear hen include more seps or make annoaions clarifing he procedure and purpose. Be sure o label and ile an included graphs or ables of daa. 1
Lab Repor - Logisics Populaion Models wih Harvesing In he following we respond o quesions associaed wih Lab 1.3 of Differenial Equaions. This repor is organized ino hree secions. In he firs secion we address a logisics growh model equaion wih consan harvesing hrough quaniaive, qualiaive and numerical analses. A similar analsis is conduced on he logisics growh model wih periodic harvesing in he second secion. The hird secion will consis of a summar of resuls and repor of conclusions associaed wih a comparison of hese wo models. 1. Logisics Growh Wih Consan Harvesing (a) Variable and Parameer Lising: dp ( d = f 1(p) = kp 1 p ) a, k, N, a R +. (3) N i. k growh parameer : This parameer describes he rae of growh of he populaion p. ii. N carring capaci : This parameer describes he oal amoun of p ha he resources can suppor. iii. a harvesing rae : This parameer describes he rae ha p will be aken from he ssem. iv. p populaion : This dependen variable describes he populaion as a funcion of ime. v. ime : This independen variable parameerizes he evoluion of he populaion, p. (b) To solve his ODE analicall we noe ha (3) is auonomous and herefore separable. Appling separaion of variables o (3) gives he following: which implies, dp ( kp 1 p n ) a = d = + C, C R (4) = ( A + B ) dp p p 1 p p (5) = A ln p p 1 + B ln p p (6) = ln (p p 1 ) A (p p ) B, (7) (p p 1 ) A (p p ) B = Ce, (8) where p 1 and p are roos o he quadraic polnomial in p and A = (p 1 p ) 1, B = (p p 1 ) 1 are found b parial fracions. If we assume he iniial populaion p() = p is given hen we find ha C = (p p 1 ) A (p p ) B. I is no, in general, clear how we should solve for p explicil. To do his we would need values for N, k, a o find p 1 and p and hus A and B. If hese numbers were known hen polnomial roo finding would give explici formula for p. (c) We now address he qualiaive informaion ha is given b he differenial equaion iself. To do his we firs find he equilibrium soluions of (3) b solving, o ge he equilibrium soluions, dp (1 d = = kp p ) a p Np + an N k p 1 () = p () = N + N N 4 an k N 4 an k =, (9), (1), (11)
assuming ha kn 4a. Oherwise, p 1 = p. We also noe ha p is no phsical for N 4 an k > N or N 4 a k < and ha for phsicall relevan cases p 1() > p () for all. To classif hese equilibria we define p 1 = p + and p = p appl linearizaion o ge, ( df dp = k 1 p ) ± (1) p=p± N ( ) = k ± 1 4a, (13) kn which implies ha p 1 = p + is a sink and p = p is a source. Figure 1.1 shows he phase line for he ssem for kn 4a and he special case where p 1 = p. (d) The following able correlaes figures o parameer choices in able 1.1 page 147. Figure Label Figure Tpe Choice Fig. 1. Slope Fields 3 Fig. 1.3 Euler s Mehod 3 Fig. 1.4 Slope Fields 4 Fig. 1.5 Euler s Mehod 4 Fig. 1.6 Slope Fields 6 Fig. 1.7 Euler s Mehod 6 Fig. 1.8 Slope Fields 8 Fig. 1.9 Euler s Mehod 8 From hese figures we can conclude ha as we increase he harvesing parameer he rajecories change giving rise o deca for paricular iniial populaions. This is due o he fac ha as a increases he equilibrium soluions of he ssem ge closer ogeher. As hese equilibria ge closer ogeher he source, p, moves up he p-axis and consequenl decaing rajecories wih iniial populaions beween zero and p are creaed. From his we can also conclude ha as he wo equilibria ge closer and closer o each oher he amoun of rajecories, which grow in ime decreases. Thus i is possible o increase he harvesing parameer o a poin where here are no iniial populaions, which grow in ime. Moreover, i is possible harves o a poin where here are no iniial populaions whose rajecor is viable in he long-erm. Using his model for a biological ssem gives insigh ino how much one could harves wihou desroing he populaion in finie-ime. This can be seen mahemaicall hrough he phase line in figure 1.1.. Logisics Growh Wih Periodic Harvesing dp ( d = f (p, ) = kp 1 p ) a(1 + sin(b)), k, N, a, b R +. (14) N (a) In his case we have ha b is a parameer, which conrols he frequenc of he periodic harvesing and a represens he overall ampliude of he periodic harvesing. We noe ha in his case he harvesing can be as much as a and as lile as. (b) Noing f (p, ) h(p)g() implies ha (14) is no separable. There is no clear wa o solve his differenial equaion analicall. 3
(c) Though here are no clear analic soluions we can sill use qualiaive and numerical echniques o he problem. The following able correlaes figures o parameer choices in able 1.1 page 147 for a = a 1 and b = 1. Figure Label Figure Tpe Choice Fig. 1.1 Slope Fields 3 Fig. 1.11 Euler s Mehod 3 Fig. 1.1 Slope Fields 4 Fig. 1.13 Euler s Mehod 4 Fig. 1.14 Slope Fields 6 Fig. 1.15 Euler s Mehod 6 Fig. 1.16 Slope Fields 8 Fig. 1.17 Euler s Mehod 8 Figures 1.18 and 1.19 show slope fields and numerical approximaions for a = a and b = 1. In general we noice similar behavior of he populaion as we did in par (1). Specificall, as we increase he ampliude of he harvesing parameer we creae rajecories are shifed in he negaive direcion. However, in his case, since he ssem is non-auonomous here are no classical equilibrium soluions. Insead, for his ssem, hose rajecories, which do no go exinc end o a sead long-erm oscillaor behavior. 3. Conclusions and Summar In his lab we are presened wo possible models for populaion harvesing. Wih hese models one can conclude ha here is a relaionship beween harvesing rae and he long-erm populaion. We, in general, noe ha increased harvesing leads o more exincion-rajecories. Moreover, i is possible o use he auonomous ssem o infer behaviors of he non-auonomous ssem. While he auonomous ssem has equilibrium soluions he non-auonomous seems o have wha could be considered sead-sae soluions. These sead-sae soluion can be used, like equilibrium soluions, as reference poins o describe he behavior of neighboring soluions. While i is clear how he presence of consan harvesing can shif an auonomous ssem s equilibrium soluion i is no obvious ha his should be rue of a non-auonomous ssem. If we compare Figure 1.1 and Figure 1.6 we see ha for paricular harvesing values boh populaions sele o a rajecor ha remains viable in he long-erm. In conras o his we have Figures 1.8 and 1.16, which show ha he effec of increased harvesing ends o lead o he loss of he populaion in finie-ime. If we hink of hese long-erm paerns as sead-sae soluion hen we can conclude ha he effec of harvesing is o shif hese sead-sae soluions, which creaes iniial populaions desined for exincion similar o he auonomous ssem. Equilibrium soluions as well as sead-sae soluions make i possible o deermine he long-erm behavior of a paricular soluion saisfing an iniial condiion. This is imporan o he phsical reali of he mahemaical model, because one would naurall like o know how much harvesing can ake place wihou decimaing a paricular populaion. Since i is unlikel ha consan harvesing is possible for a paricular populaion i makes sense o consider periodic harvesing insead. While here are clearl qualiaive similariies in he models here are some ineresing quaniaive differences. If we compare Figure 1.13 o Figure 1.5 we see ha oscillaor harvesing has a derimenal effec causing he purple rajecor o go exinc more quickl han in he consan harvesing case. l, sinusoidal harvesing, while more complicaed, does no guaranee he long-erm viabili of all iniial populaions. This is because he harvesing does no ake ino accoun he curren populaion a he ime of harves. A wa o avoid his deca o exincion is o harves he populaion a a rae proporional o he populaion iself. In his wa he harvesing rae will decrease when he populaion is low and possibl migh save a populaion from exincion. This mehod is oulined in problems and 1 of chaper one secion seven of he ex and shows ha his more complicaed auonomous model gives beer(safer) resuls when a populaion ges ver low. 4
Figures Figure 1.1 Slope Line for kn 4a. N + N 4 an k p 1 () = N N 4 an k p () =, sink, source Slope Line for kn = 4a. p() = N, node In general, i is no obvious wh his case ields a node. The bes wa o hink abou i is he following: 1. Think of f 1 (p) as a parabola ha opens downward in he f 1 plane and has roos p 1 and p.. As a increases he parabola is shifed in he negaive f 1 direcion. This causes he roos o ge closer ogeher. 3. When kn = 4a he roos join and become a double roo and f 1 is negaive on boh sides of his roo. This implies ha he roo (aka equilibrium soluion) is a node. You ma wan o draw a picure of his o convince ourself. I is ineresing o noe ha his sor of change in he equilibria of an auonomous ODE wih respec o a parameer is called a bifurcaion, specificall a saddle-node bifurcaion, and is highl imporan in he sud of ODE s and heir connecions o mahemaical chaos. We do no have ime for his secion, 1.7 of our ex. I is accessible wih he background we have and hose ineresed should check i ou. 5
Figure 1. 5 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.1 max 5 max 1.5 dela.5 5 1 = 1.134 = -.74.**(1-/5)-.1 d/d= -.38 Equaions Runge Kua 4 Draw s max 5 max 1.5 dela.5 6
Figure 1.3 6 5 4 3 1 4 6 8 1 1 - -3-4 -5 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.1 3.3.3 1 -.5.5 -.1 1.15 4.97375 -.478.15 3.375.9849.15.375.3149.15.99375 -.575.15.485 -.141.5 4.948158 -.19974.5 3.7481.9699.5.7519.398.5.98746 -.515.5.469699 -.1488 3.375 4.9318554 -.19487.375 3.11193.9547.375.1136.3447.375.98966 -.53.375.45488 -.1743 4.5 4.898864 -.1917.5 3.14887.9396.5.1511.3595.5.97449 -.539.5.438159 -.135 5.65 4.8755456 -.18564.65 3.18561.944.65.18936.3743.65.96779 -.5391.65.4194 -.1374 6.75 4.8518514 -.1815.75 3.17.99.75.779.389.75.96154 -.5473.75.45311 -.13551 7.875 4.8919415 -.1771.875 3.5853.8939.875.664.3137.875.9541 -.5558.875.388373 -.13836 8 1 4.87684 -.179 1 3.9471.8786 1.35.31184 1.94765 -.5644 1.37178 -.1419 9 1.15 4.7854556 -.16893 1.15 3.3369.8633 1.15.34418.3139 1.15.941 -.573 1.15.353416 -.14431 1 1.5 4.7643396 -.1659 1.5 3.36648.848 1.5.38334.31475 1.5.93345 -.581 1.5.335377 -.1474 11 1.375 4.74379 -.16137 1.375 3.48.837 1.375.469.31619 1.375.95768 -.5913 1.375.316949 -.1563 1 1.5 4.7353189 -.15776 1.5 3.43749.8173 1.5.461.31763 1.5.918377 -.66 1.5.981 -.15393 13 1.65 4.738114 -.1547 1.65 3.4771.8 1.65.5191.3197 1.65.9187 -.611 1.65.78879 -.15734 14 1.75 4.6845749 -.1589 1.75 3.5773.7866 1.75.5418.35 1.75.9343 -.6199 1.75.591 -.1685 15 1.875 4.66566669 -.1476 1.875 3.5456.771 1.875.58186.319 1.875.895495 -.698 1.875.3917 -.16447 16 4.6471613 -.1444 3.577.7558.61.3334.88763 -.6399.18548 -.168 17.15 4.6916345 -.14133.15 3.61165.744.15.665.3474.15.87964 -.65.15.19753 -.176 18.5 4.61149676 -.13834.5 3.64591.749.5.7311.3615.5.871496 -.668.5.17616 -.1764 19.375 4.594467 -.13543.375 3.67997.795.375.74388.3754.375.86336 -.6716.375.1541 -.1815.5 4.57776 -.136.5 3.71384.6941.5.7848.3893.5.854841 -.686.5.131493 -.18439 1.65 4.5678 -.1986.65 3.74751.6786.65.8594.3331.65.84638 -.6939.65.18444 -.18878.75 4.54446838 -.1719.75 3.781.663.75.8673.33168.75.837634 -.754.75.84847 -.1933 3.875 4.5856913 -.146.875 3.8149.6478.875.9869.3334.875.88817 -.7171.875.668 -.1981 4 3 4.5199366 -.19 3 3.84738.633 3.953.3344 3.819853 -.79 3.3593 -.87 5 3.15 4.4977395 -.11964 3.15 3.889.6169 3.15.991.33575 3.15.81738 -.7414 3.15.157 -.789 6 3.5 4.487786 -.1176 3.5 3.913.615 3.5.1349.3379 3.5.8147 -.754 3.5 -.1541 -.139 7 3.375 4.468111 -.11494 3.375 3.9455.586 3.375.176.3384 3.375.7945 -.7668 3.375 -.45 -.1848 8 3.5 4.45375371 -.1169 3.5 3.97784.576 3.5.11185.33974 3.5.7846 -.78 3.5 -.6936 -.46 9 3.65 4.43966794 -.1149 3.65 3.1998.555 3.65.11699.3415 3.65.7771 -.7934 3.65 -.9737 -.985 3 3.75 4.4585638 -.1836 3.75 3.1419.5398 3.75.136.3435 3.75.76793 -.87 3.75 -.161 -.3586 31 3.875 4.4131177 -.168 3.875 3.17366.544 3.875.14641.34364 3.875.7573 -.81 3.875 -.15558 -.48 3 4 4.39979 -.145 4 3.115.591 4.18937.3449 4.74438 -.8356 4 -.18584 -.4855 33 4.15 4.38599557 -.18 4.15 3.113658.4937 4.15.13349.346 4.15.731993 -.853 4.15 -.1691 -.556 34 4.5 4.373167 -.136 4.5 3.116775.4784 4.5.137576.34746 4.5.71363 -.8654 4.5 -.488 -.64 35 4.375 4.366668 -.9848 4.375 3.119873.463 4.375.141919.34871 4.375.71546 -.889 4.375 -.816 -.6949 36 4.5 4.34835569 -.9666 4.5 3.195.4477 4.5.14678.34995 4.5.699535 -.8967 4.5 -.3159 -.773 37 4.65 4.3367359 -.9488 4.65 3.161.434 4.65.15653.35118 4.65.68837 -.919 4.65 -.3499 -.8488 38 4.75 4.344149 -.9314 4.75 3.195.417 4.75.1554.354 4.75.676916 -.995 4.75 -.38553 -.935 39 4.875 4.3177166 -.9145 4.875 3.1374.419 4.875.159447.35361 4.875.66598 -.9465 4.875 -.416 -.3156 4 5 4.313495 -.8979 5 3.13576.3867 5.163867.35481 5.653467 -.9639 5 -.45985 -.3143 41 5.15 4.91168 -.8818 5.15 3.13859.3715 5.15.16833.35599 5.15.641419 -.9817 5.15 -.49866 -.31968 4 5.5 4.79941 -.8661 5.5 3.1414.3564 5.5.1775.35716 5.5.69147 -.1 5.5 -.5386 -.3933 43 5.375 4.686833 -.857 5.375 3.143969.341 5.375.17717.3583 5.375.616646 -.1188 5.375 -.57978 -.3394 44 5.5 4.5763447 -.8357 5.5 3.146896.361 5.5.181696.35947 5.5.63911 -.1381 5.5 -.61 -.34993 45 5.65 4.471887 -.811 5.65 3.14983.311 5.65.186189.3661 5.65.59936 -.1578 5.65 -.66595 -.3693 46 5.75 4.369474 -.868 5.75 3.1569.96 5.75.19697.36173 5.75.577713 -.1781 5.75 -.7116 -.3744 47 5.875 4.684 -.798 5.875 3.15556.81 5.875.19519.3684 5.875.56437 -.1989 5.875 -.7576 -.38448 48 6 4.16933 -.7791 6 3.158413.66 6.199754.36394 6.5551 -.11 6 -.8568 -.3971 49 6.15 4.71917 -.7658 6.15 3.16146.51 6.15.434.3653 6.15.536498 -.1141 6.15 -.8553 -.4133 5 6.5 4.19761856 -.758 6.5 3.1646.361 6.5.8866.3661 6.5.5 -.11646 6.5 -.9661 -.44 51 6.375 4.1889 -.74 6.375 3.166855.1 6.375.13443.36715 6.375.57664 -.11878 6.375 -.95963 -.43876 5 6.5 4.17895877 -.776 6.5 3.169631.64 6.5.183.368 6.5.49817 -.1115 6.5.1448 -.4546 7
Figure 1.4 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.16 max 1 max 1.5 dela.5 1 = =.5 1.**(1-/5)-.16 d/d= -.7 Equaions Runge Kua 4 Draw s max 1 max 1.5 dela.5 8
Figure 1.5 6 5 4 3 1 4 6 8 1 1 - DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.16 3.8.8 1.5 -.7 1.15 4.98 -.156.15 3.1.79596.15.1.8396.15 1.15.4915 -.714.5 4.96498 -.1516.5 3.1995.79186.5.5.8786.5 1.5.4835 -.784 3.375 4.94147775 -.14843.375 3.9848.7877.375.3148.8117.375 1.375.473 -.7431 4.5 4.99368 -.1448.5 3.39694.78349.5.494.81547.5 1.5.46393 -.758 5.65 4.94888 -.14133.65 3.49488.7793.65.5487.81918.65 1.65.45445 -.7737 6.75 4.8871557 -.13794.75 3.598.77491.75.677.88.75 1.75.444781 -.7896 7.875 4.869915 -.13466.875 3.68914.7753.875.711.8639.875 1.875.43491 -.858 8 1 4.85389 -.13148 1 3.78546.76611 1.8134.8989 1 1 1.44839 -.85 9 1.15 4.83664516 -.184 1.15 3.881.76164 1.15.91716.8333 1.15 1 1.15.414557 -.8396 1 1.5 4.859561 -.1541 1.5 3.97643.75713 1.5.113.83668 1.5 1 1.5.446 -.857 11 1.375 4.8491979 -.151 1.375 3.1717.7557 1.375.11591.83997 1.375 1 1.375.393347 -.875 1 1.5 4.7896651 -.11969 1.5 3.116514.74796 1.5.139.84318 1.5 1 1.5.3847 -.8937 13 1.65 4.774645 -.11696 1.65 3.15864.7433 1.65.13363.84631 1.65 1 1.65.37136 -.917 14 1.75 4.76497 -.11431 1.75 3.135155.73863 1.75.1449.84937 1.75 1 1.75.35988 -.931 15 1.875 4.74573641 -.11173 1.875 3.144388.73391 1.875.15486.8534 1.875 1 1.875.348176 -.951 16 4.73176975 -.193 3.15356.7914.16548.8554 1.33674 -.977 17.15 4.71811577 -.168.15 3.16676.7434.15.176171.8585.15 1.15.34116 -.9938 18.5 4.7476558 -.1444.5 3.171731.71951.5.186896.8679.5 1.5.311693 -.1155 19.375 4.691716 -.114.375 3.1874.71465.375.197656.86344.375 1.375.99 -.1378.5 4.6789465 -.9991.5 3.189658.7975.5.8449.866.5 1.5.868 -.167 1.65 4.66645369 -.9774.65 3.19859.748.65.1974.86848.65 1.65.777 -.184.75 4.654368 -.9563.75 3.734.69987.75.313.8787.75 1.75.5917 -.1184 3.875 4.6484 -.9357.875 3.1688.69489.875.4116.87317.875 1.875.45361 -.11334 4 3 4.6358555 -.9158 3 3.4774.68988 3.51931.87538 3 1 3.31194 -.1159 5 3.15 4.61913857 -.8963 3.15 3.33398.68485 3.15.6873.87751 3.15 1 3.15.1677 -.11854 6 3.5 4.6793483 -.8774 3.5 3.41958.6798 3.5.7384.87954 3.5 1 3.5.189 -.115 7 3.375 4.59696789 -.8589 3.375 3.5456.67473 3.375.84836.88148 3.375 1 3.375.186733 -.145 8 3.5 4.586315 -.849 3.5 3.5889.66963 3.5.95854.88333 3.5 1 3.5.1717 -.1693 9 3.65 4.57571971 -.834 3.65 3.676.6645 3.65.36896.8858 3.65 1 3.65.155361 -.1989 3 3.75 4.5654664 -.864 3.75 3.75567.6594 3.75.31796.88674 3.75 1 3.75.13915 -.1395 31 3.875 4.55534671 -.7898 3.875 3.8389.6546 3.875.3944.88831 3.875 1 3.875.156 -.1361 3 4 4.54547446 -.7736 4 3.91988.6491 4.34148.88978 4 1 4.15494 -.13935 33 4.15 4.5358463 -.7578 4.15 3.311.64394 4.15.3517.89115 4.15 1 4.15.8875 -.147 34 4.5 4.563313 -.744 4.5 3.3815.63876 4.5.3649.8943 4.5 1 4.5.739 -.14615 35 4.375 4.5175 -.774 4.375 3.316135.63357 4.375.373565.89361 4.375 1 4.375.5197 -.14971 36 4.5 4.579595 -.718 4.5 3.3455.6837 4.5.384735.89469 4.5 1 4.5.3356 -.15339 37 4.65 4.4995 -.6985 4.65 3.33199.6317 4.65.395918.89567 4.65 1 4.65.1481 -.15719 38 4.75 4.493191 -.6845 4.75 3.339699.61796 4.75.47114.89655 4.75 1 4.75 -.557 -.16111 39 4.875 4.4817616 -.671 4.875 3.34743.6175 4.875.41831.89733 4.875 1 4.875 -.571 -.16517 4 5 4.4733756 -.6577 5 3.35583.6753 5.49538.8981 5 1 5 -.4635 -.16936 41 5.15 4.4651541 -.6447 5.15 3.36677.63 5.15.44763.8986 5.15 1 5.15 -.675 -.17369 4 5.5 4.457955 -.631 5.5 3.376.5971 5.5.451995.8998 5.5 1 5.5 -.893 -.17817 43 5.375 4.44919395 -.6197 5.375 3.37767.59188 5.375.46334.89946 5.375 1 5.375 -.1115 -.188 44 5.5 4.44144716 -.677 5.5 3.38568.58666 5.5.474477.89974 5.5 1 5.5 -.13435 -.18759 45 5.65 4.4338518 -.5959 5.65 3.3941.58145 5.65.48574.8999 5.65 1 5.65 -.1578 -.1956 46 5.75 4.46418 -.5844 5.75 3.399669.5764 5.75.496973.9 5.75 1 5.75 -.18187 -.1977 47 5.875 4.419975 -.573 5.875 3.4687.5713 5.875.583.89997 5.875 1 5.875 -.658 -.3 48 6 4.4119338 -.56 6 3.4141.56583 6.51947.89985 6 1 6 -.3196 -.854 49 6.15 4.449495 -.5515 6.15 3.4183.5664 6.15.5371.8996 6.15 1 6.15 -.583 -.147 5 6.5 4.3981164 -.541 6.5 3.4891.55546 6.5.541966.8993 6.5 1 6.5 -.8481 -.1 51 6.375 4.391494 -.537 6.375 3.43535.558 6.375.5537.89887 6.375 1 6.375 -.3134 -.637 5 6.5 4.3846158 -.57 6.5 3.441913.5451 6.5.564443.89834 6.5 1 6.5 -.3464 -.377 9
Figure 1.6 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.9 max 1 max 1.5 dela.5 1 1 =.718 =.83.**(1-/5)-.9 d/d= -.37 Equaions Runge Kua 4 Draw s max 1 max 1.5 dela.5 1
Figure 1.5 6 5 4 3 1 4 6 8 1 1 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.9 3.15.15 1.7.5 1.15 4.98875 -.8776.15 3.1875.14936.15.1875.15736.15 1.875.7147.15.5.5 4.977786 -.8558.5 3.3744.148448.5.3759.151447.5 1.17631.713.5.5 3.375 4.9678363 -.8346.375 3.5596.147636.375.5653.15133.375 1.6644.73169.375.5 4.5 4.9566511 -.8141.5 3.74415.1468.5.7554.15793.5 1.3579.7444.5.5 5.65 4.94647545 -.7941.65 3.9765.145945.65.94639.15347.65 1.457.7537.65.5 6.75 4.9365494 -.7747.75 3.1118.14567.75.113817.15435.75 1.54486.764.75.5 7.875 4.9686538 -.7559.875 3.1914.144167.875.13371.154615.875 1.6439.7751.875.5 8 1 4.917417 -.7376 1 3.147163.14347 1.15398.155167 1 1.7379.7863 1.5 9 1.15 4.9819748 -.7198 1.15 3.16569.1437 1.15.171794.155691 1.15 1.83558.79748 1.15.5 1 1.5 4.8994 -.75 1.5 3.18857.141348 1.5.19156.156187 1.5 1.9356.8873 1.5.5 11 1.375 4.8941959 -.6856 1.375 3.55.14371 1.375.1779.156654 1.375 1.13635.87 1.375.5 1 1.5 4.8818496 -.6693 1.5 3.187.139375 1.5.3361.1579 1.5 1.113886.83148 1.5.5 13 1.65 4.8734833 -.6534 1.65 3.35494.13836 1.65.49997.1575 1.65 1.1479.8496 1.65.5 14 1.75 4.86531593 -.6379 1.75 3.5789.13733 1.75.69685.157878 1.75 1.134816.85451 1.75.5 15 1.875 4.8573433 -.68 1.875 3.69955.13687 1.875.89419.1586 1.875 1.145498.86613 1.875.5 16 4.8495571 -.68 3.86991.1356.39198.158544 1.15634.87781.5 17.15 4.8419549 -.5939.15 3.33895.13415.15.3916.158831.15 1.16797.88956.15.5 18.5 4.83453116 -.58.5 3.3663.1336.5.348869.15986.5 1.178417.9137.5.5 19.375 4.87898 -.5665.375 3.33796.131957.375.368755.159311.375 1.189684.9133.375.5.5 4.81998 -.5533.5 3.353791.1384.5.388669.15954.5 1.199.9514.5.5 1.65 4.8138316 -.545.65 3.37146.19714.65.4867.159666.65 1.1663.93711.65.5.75 4.865677 -.58.75 3.38636.18575.75.48565.159796.75 1.4377.94911.75.5 3.875 4.7999644 -.5159.875 3.443.1745.875.44854.159894.875 1.3641.96117.875.5 4 3 4.79347813 -.54 3 3.41836.1665 3.46857.15996 3 1.4856.9735 3.5 5 3.15 4.7871779 -.495 3.15 3.434143.1595 3.15.4885.159995 3.15 1.641.98538 3.15.5 6 3.5 4.7811 -.481 3.5 3.44978.13917 3.5.5851.159997 3.5 1.7738.99753 3.5.5 7 3.375 4.77567 -.473 3.375 3.46569.173 3.375.5851.159967 3.375 1.858.1971 3.375.5 8 3.5 4.769184 -.4596 3.5 3.48611.11536 3.5.548517.15996 3.5 1.9789.1191 3.5.5 9 3.65 4.7633837 -.449 3.65 3.49583.1335 3.65.56855.15981 3.65 1.3163.13413 3.65.5 3 3.75 4.75776917 -.439 3.75 3.51845.11918 3.75.588481.159687 3.75 1.3353.14637 3.75.5 31 3.875 4.758157 -.491 3.875 3.55736.117915 3.875.6844.15953 3.875 1.33669.15861 3.875.5 3 4 4.7469177 -.4195 4 3.54475.116696 4.68383.159341 4 1.34984.1785 4.5 33 4.15 4.74167451 -.41 4.15 3.5556.115474 4.15.64831.1591 4.15 1.3637.1831 4.15.5 34 4.5 4.73654899 -.49 4.5 3.569496.11447 4.5.668191.158868 4.5 1.376766.19534 4.5.5 35 4.375 4.7315383 -.3919 4.375 3.583777.11317 4.375.6885.158585 4.375 1.39458.11757 4.375.5 36 4.5 4.7663941 -.383 4.5 3.59794.111784 4.5.77873.1587 4.5 1.4433.111978 4.5.5 37 4.65 4.718498 -.3746 4.65 3.611877.11549 4.65.77657.15797 4.65 1.4183.113197 4.65.5 38 4.75 4.7171667 -.3663 4.75 3.65696.1931 4.75.747398.15755 4.75 1.43449.114413 4.75.5 39 4.875 4.7158758 -.358 4.875 3.63936.1874 4.875.7679.157146 4.875 1.446751.11567 4.875.5 4 5 4.781986 -.353 5 3.65869.16836 5.786735.156711 5 1.4614.116836 5.5 41 5.15 4.7373111 -.346 5.15 3.6664.15597 5.15.8634.15647 5.15 1.47589.11841 5.15.5 4 5.5 4.69944896 -.335 5.5 3.67943.14358 5.5.85855.155753 5.5 1.49564.1194 5.5.5 43 5.375 4.695618 -.377 5.375 3.69468.1311 5.375.84534.1553 5.375 1.55469.1436 5.375.5 44 5.5 4.691165 -.35 5.5 3.75358.11884 5.5.86477.154679 5.5 1.554.1165 5.5.5 45 5.65 4.687159 -.3135 5.65 3.71894.165 5.65.8846.1541 5.65 1.53577.187 5.65.5 46 5.75 4.683487 -.366 5.75 3.73675.99418 5.75.9335.153493 5.75 1.55178.1398 5.75.5 47 5.875 4.6794817 -.999 5.875 3.7431.98188 5.875.9511.15859 5.875 1.566576.15149 5.875.5 48 6 4.6756597 -.934 6 3.755376.96961 6.941619.15199 6 1.5819.1637 6.5 49 6.15 4.6719916 -.87 6.15 3.767496.95738 6.15.96644.15151 6.15 1.5988.17456 6.15.5 5 6.5 4.6684387 -.88 6.5 3.779463.94519 6.5.979583.158 6.5 1.61394.18596 6.5.5 51 6.375 4.66489399 -.747 6.375 3.79178.9334 6.375.998433.1563 6.375 1.6314.1975 6.375.5 5 6.5 4.6614616 -.688 6.5 3.8941.994 6.5 3.17191.14931 6.5 1.6463.13843 6.5.5 11
Figure 1.8 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.4 max 1 max 1.5 dela.5 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.4 max 1 max 1.5 dela.5 1
Figure 1.9 6 4 4 6 8 1 1 - -4-6 -8 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.4 3 1 -.8.5 -.15 1.15 4.97 -.344.15 3.15.15.99 -.81.15.4815 -.1531.5 4.947455 -.89.5 3.5.5.97985 -.843.5.4613 -.1561 3.375 4.91931 -.75.375 3.375.375.969545 -.8369.375.4469 -.15931 4.5 4.88436554 -.1741.5 3.5.5.95984 -.8498.5.4694 -.1661 5.65 4.85718954 -.15.65 3.65.65.94846 -.869.65.4368 -.166 6.75 4.8365783 -.78.75 3.75.75.937675 -.8763.75.381618 -.1695 7.875 4.84748 -.47.875 3.875.875.9671 -.891.875.3643 -.17311 8 1 4.77943868 -.19783 1 3 1 1.915595 -.941 1.338791 -.17683 9 1.15 4.7547948 -.19335 1.15 3 1.15 1.15.9493 -.9185 1.15.316687 -.1867 1 1.5 4.735491 -.1891 1.5 3 1.5 1.5.8981 -.933 1.5.9413 -.18464 11 1.375 4.7691434 -.1848 1.375 3 1.375 1.375.881147 -.9483 1.375.713 -.18873 1 1.5 4.68381199 -.1876 1.5 3 1.5 1.5.86993 -.9637 1.5.47431 -.1996 13 1.65 4.6611681 -.17683 1.65 3 1.65 1.65.85747 -.9795 1.65.3311 -.19733 14 1.75 4.639115 -.1733 1.75 3 1.75 1.75.8454 -.9956 1.75.198645 -.185 15 1.875 4.61748351 -.16935 1.875 3 1.875 1.875.83559 -.111 1.875.173413 -.65 16 4.59631483 -.16578 3.81997 -.191.147598 -.1135 17.15 4.5755915 -.163.15 3.15.15.8744 -.1464.15.11179 -.1635 18.5 4.5553174 -.15897.5 3.5.5.793963 -.164.5.94135 -.153 19.375 4.5354341 -.1557.375 3.375.375.7866 -.185.375.66445 -.689.5 4.51596553 -.1556.5 3.5.5.76713 -.1111.5.3884 -.344 1.65 4.49689494 -.1495.65 3.65.65.753365 -.113.65.98 -.38.75 4.478699 -.14653.75 3.75.75.73936 -.11399.75 -.75 -.4417 3.875 4.4598948 -.14365.875 3.875.875.7511 -.1161.875 -.517 -.536 4 3 4.44193463 -.1484 3 3 3 3.71611 -.1188 3 -.856 -.5679 5 3.15 4.44397 -.1381 3.15 3 3.15 3.15.69585 -.1 3.15 -.11466 -.6346 6 3.5 4.476386 -.13548 3.5 3 3.5 3.5.6887 -.138 3.5 -.14759 -.739 7 3.375 4.39194 -.139 3.375 3 3.375 3.375.66553 -.1461 3.375 -.18139 -.7759 8 3.5 4.37351645 -.134 3.5 3 3.5 3.5.649954 -.1691 3.5 -.169 -.859 9 3.65 4.3571613 -.1797 3.65 3 3.65 3.65.6349 -.196 3.65 -.5173 -.988 3 3.75 4.3411988 -.156 3.75 3 3.75 3.75.61793 -.13169 3.75 -.8834 -.399 31 3.875 4.355194 -.133 3.875 3 3.875 3.875.61471 -.13418 3.875 -.3596 -.3944 3 4 4.31168 -.116 4 3 4 4.584699 -.13674 4 -.36464 -.3185 33 4.15 4.9497438 -.11888 4.15 3 4.15 4.15.56767 -.13937 4.15 -.444 -.3743 34 4.5 4.81147 -.11675 4.5 3 4.5 4.5.55187 -.147 4.5 -.44535 -.337 35 4.375 4.65568 -.11468 4.375 3 4.375 4.375.5348 -.14485 4.375 -.48747 -.347 36 4.5 4.5118536 -.1167 4.5 3 4.5 4.5.51431 -.1477 4.5 -.5385 -.35744 37 4.65 4.37111 -.117 4.65 3 4.65 4.65.495857 -.1566 4.65 -.57553 -.36836 38 4.75 4.36449 -.1879 4.75 3 4.75 4.75.4774 -.1537 4.75 -.6157 -.37977 39 4.875 4.96669 -.169 4.875 3 4.875 4.875.457811 -.1568 4.875 -.6695 -.39171 4 5 4.1963149 -.151 5 3 5 5.4389 -.164 5 -.7181 -.44 41 5.15 4.1831643 -.133 5.15 3 5.15 5.15.4184 -.16336 5.15 -.76854 -.41733 4 5.5 4.17499 -.1159 5.5 3 5.5 5.5.397784 -.16677 5.5 -.87 -.4318 43 5.375 4.1575543 -.999 5.375 3 5.375 5.375.376938 -.173 5.375 -.87459 -.44551 44 5.5 4.145636 -.985 5.5 3 5.5 5.5.355651 -.17393 5.5 -.938 -.4667 45 5.65 4.137819 -.9664 5.65 3 5.65 5.65.33391 -.17768 5.65 -.98786 -.47661 46 5.75 4.17 -.957 5.75 3 5.75 5.75.3117 -.18155 5.75.4744 -.49337 47 5.875 4.1881864 -.9353 5.875 3 5.875 5.875.897 -.18554 5.875.1911 -.5113 48 6 4.971715 -.93 6 3 6 6.65814 -.18966 6.1799 -.5963 49 6.15 4.85638 -.957 6.15 3 6.15 6.15.416 -.1939 6.15.3919 -.5496 5 6.5 4.7438 -.8914 6.5 3 6.5 6.5.17866 -.19833 6.5.3785 -.56999 51 6.375 4.6315994 -.8774 6.375 3 6.375 6.375.19375 -.88 6.375.3791 -.5919 5 6.5 4.5196 -.8637 6.5 3 6.5 6.5.167716 -.758 6.5.4539 -.6158 13
Figure 1.1 5 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.1*(1+sin()) max 5 max 1.5 dela.5 5 = 9.56 =.7 1.**(1-/5)-.1*(1+sin()) d/d= -.9 Equaions Runge Kua 4 Draw s max 5 max 1.5 dela.5 14
Figure 1.11 6 4 4 6 8 1 1 - -4-6 -8 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.1 3.3.3 1 -.5.5 -.1 1.15 4.97375 -.396.15 3.375.3668.15.375.3968.15.99375 -.7693.15.485 -.14859.5 4.944889 -.515.5 3.48 -.1.5.446 -.179.5.984133 -.1387.5.46646 -.17737 3.375 4.91349855 -.699.375 3.1443 -.4698.375.153 -.4686.375.97115 -.1341.375.44455 -.596 4.5 4.879759 -.871.5.995571 -.75.5 1.995666 -.785.5.954848 -.15618.5.41851 -.3398 5.65 4.84385783 -.36.65.986758 -.935.65 1.98689 -.9341.65.93536 -.188.65.3896 -.618 6.75 4.86368 -.31586.75.97515 -.1118.75 1.975133 -.11416.75.9176 -.39.75.35667 -.8691 7.875 4.76654877 -.3667.875.961193 -.1969.875 1.96863 -.1381.875.88736 -.5.875.3764 -.31115 8 1 4.7571454 -.33486 1.944981 -.14463 1 1.9446 -.1496 1.85983 -.4441 1.8187 -.33351 9 1.15 4.6838569 -.345 1.15.969 -.15677 1.15 1.9569 -.1667 1.15.88531 -.613 1.15.4181 -.35375 1 1.5 4.64136 -.347 1.5.9737 -.1659 1.5 1.9595 -.17343 1.5.795877 -.7545 1.5.195963 -.37163 11 1.375 4.59848899 -.3413 1.375.886566 -.17196 1.375 1.883615 -.18118 1.375.761446 -.8689 1.375.14959 -.38698 1 1.5 4.55576 -.33851 1.5.86571 -.17481 1.5 1.86967 -.18581 1.5.75585 -.954 1.5.11136 -.39966 13 1.65 4.513485 -.33184 1.65.843 -.1744 1.65 1.837741 -.1874 1.65.688658 -.393 1.65.51179 -.4956 14 1.75 4.4719754 -.318 1.75.814 -.1777 1.75 1.814337 -.18544 1.75.6514 -.3338 1.75.6E-5 -.41664 15 1.875 4.4316554 -.3961 1.875.873 -.16396 1.875 1.791157 -.1846 1.875.613119 -.377 1.875 -.51 -.489 16 4.39954 -.948.779578 -.1548 1.7686 -.1735.57573 -.9914 -.1471 -.433 17.15 4.35616877 -.7638.15.76319 -.1418.15 1.74756 -.1614.15.537881 -.956.15 -.1575 -.416 18.5 4.31616 -.5613.5.74659 -.1575.5 1.769 -.1473.5.5131 -.8319.5 -.113 -.41719 19.375 4.89651 -.3378.375.7694 -.1773.375 1.78487 -.1373.375.46591 -.7117.375 -.68 -.4188.5 4.63846 -.964.5.713473 -.875.5 1.69146 -.11178.5.4315 -.5674.5 -.31364 -.434 1.65 4.3417783 -.184.65.7536 -.6536.65 1.678173 -.974.65.39993 -.414.65 -.36393 -.39181.75 4.1117556 -.1577.75.694365 -.4166.75 1.66683 -.679.75.36996 -.164.75 -.4191 -.37955 3.875 4.19151635 -.1977.875.689158 -.1675.875 1.658341 -.4366.875.34 -.157.875 -.4635 -.36587 4 3 4.1759476 -.119 3.68763.8965 3 1.65883 -.1834 3.3174 -.185 3 -.569 -.3511 5 3.15 4.165579 -.745 3.15.688184.3599 3.15 1.65591.7656 3.15.94473 -.1586 3.15 -.54997 -.33558 6 3.5 4.153317 -.466 3.5.69571.6138 3.5 1.651548.3396 3.5.74715 -.13535 3.5 -.5919 -.31968 7 3.375 4.14747437 -. 3.375.76.86968 3.375 1.655789.664 3.375.57796 -.1153 3.375 -.63188 -.3378 8 3.5 4.14497497.547 3.5.71197.11188 3.5 1.66397.8566 3.5.4373 -.8997 3.5 -.66985 -.885 9 3.65 4.1456533.981 3.65.758.135581 3.65 1.6744.11317 3.65.3484 -.686 3.65 -.7588 -.735 3 3.75 4.14931341.5118 3.75.743.157685 3.75 1.687794.133641 3.75.3977 -.4718 3.75 -.747 -.5989 31 3.875 4.1557157.7919 3.875.76174.177835 3.875 1.74499.1556 3.875.1879 -.771 3.875 -.7756 -.4781 3 4 4.1645863.8895 4.78397.19573 4 1.7397.174836 4.14615 -.999 4 -.8353 -.376 33 4.15 4.1755955.1498 4.15.88433.1997 4.15 1.745761.1947 4.15.13367.5654 4.15 -.8333 -.96 34 4.5 4.1884475.113919 4.5.83487.3464 4.5 1.769767.6618 4.5.1473.1899 4.5 -.86194 -.416 35 4.375 4.64467.11 4.375.8674.3897 4.375 1.795595.18313 4.375.1644.9575 4.375 -.88996 -.151 36 4.5 4.1791976.1731 4.5.89185.39139 4.5 1.8884.694 4.5.136.3737 4.5 -.91764 -.193 37 4.65 4.338368.18953 4.65.91745.484 4.65 1.85151.3364 4.65.488.4139 4.65 -.94539 -.563 38 4.75 4.499479.17359 4.75.955.41679 4.75 1.8897.3449 4.75.375.43749 4.75 -.97359 -.378 39 4.875 4.65867.1499 4.875.9815.3798 4.875 1.9968.3387 4.875.35544.4119 4.875.69 -.435 4 5 4.811758.114471 5 3.11956.389 5 1.938769.8775 5.489.3716 5.3313 -.5795 41 5.15 4.954839.1346 5.15 3.4818.677 5.15 1.967366.18 5.15.45461.958 5.15.6537 -.761 4 5.5 4.38411.8956 5.5 3.684.7453 5.5 1.994994.1175 5.5.4993.17713 5.5.9988 -.9799 43 5.375 4.3196738.7313 5.375 3.94334.19143 5.375.166.196394 5.375.5137.397 5.375.13713 -.3359 44 5.5 4.3874773.54391 5.5 3.11863.17874 5.5.45816.17991 5.5.51719 -.143 5.5.17758 -.358 45 5.65 4.33554657.33684 5.65 3.13987.1576 5.65.6835.16999 5.65.49966 -.345 5.65.168 -.38558 46 5.75 4.3397577.1135 5.75 3.158881.19374 5.75.8849.139963 5.75.4579 -.5653 5.75.6988 -.4174 47 5.875 4.34117587 -.14 5.875 3.17553.1513 5.875.1595.117147 5.875.3864 -.8119 5.875.36 -.46113 48 6 4.339646 -.3669 6 3.188194.79733 6.1568.9919 6.8493 -.1771 6.384 -.5357 49 6.15 4.3355919 -.616 6.15 3.198161.53583 6.15.13183.67669 6.15.159 -.13576 6.15.44319 -.54887 5 6.5 4.373574 -.866 6.5 3.4859.795 6.5.1464.418 6.5.19859 -.16499 6.5.51179 -.59681 51 6.375 4.31653183 -.1115 6.375 3.846.68 6.375.145867.1573 6.375.177435 -.1953 6.375.5864 -.647 5 6.5 4.36614 -.13515 6.5 3.8331 -.54 6.5.147833 -.114 6.5.15357 -.55 6.5.6673 -.69983 15
Figure 1.1 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.16*(1+sin()) max 1 max 1.5 dela.5 1 = =.5 1.**(1-/5)-.16*(1+sin()) d/d= -.7 Equaions Runge Kua 4 Draw s max 1 max 1.5 dela.5 16
Figure 1.13 6 5 4 3 1 4 6 8 1 1 - -3-4 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.16 3.8.8 1.5 -.7 1.15 4.98 -.17596.15 3.1.59648.15.1.6448.15 1 -.1995.15.4915 -.9135.5 4.958451 -.1916.5 3.17456.3975.5.17556.4115.5.99757 -.3988.5.479831 -.1183 3.375 4.934975 -.56.375 3.419.48.375.694.84.375.9951 -.595.375.46578 -.13413 4.5 4.9839789 -.187.5 3.4979.68.5.548.485.5.98583 -.7851.5.448961 -.15498 5.65 4.8815748 -.339.65 3.563 -.1465.65.615 -.16.65.9757 -.9661.65.49589 -.1758 6.75 4.855836 -.439.75 3.3431 -.3.75.444 -.811.75.963194 -.11353.75.4774 -.19417 7.875 4.8999 -.4851.875 3.19678 -.4361.875.96 -.4199.875.949 -.193.875.38343 -.1 8 1 4.7911458 -.5461 1 3.147 -.551 1.15678 -.54 1.93873 -.1487 1.35693 -.834 9 1.15 4.75931966 -.5854 1.15 3.736 -.6466 1.15.896 -.641 1.15.91514 -.15485 1.15.38389 -.43 1 1.5 4.77168 -.6 1.5.99943 -.7181 1.5.95 -.718 1.5.895658 -.16479 1.5.9814 -.5579 11 1.375 4.6944743 -.5957 1.375.9967 -.7656 1.375 1.99195 -.777 1.375.87559 -.1756 1.375.6641 -.6657 1 1.5 4.66785 -.5657 1.5.98698 -.7884 1.5 1.9891 -.83 1.5.853489 -.1784 1.5.37 -.75 13 1.65 4.6995613 -.513 1.65.97843 -.7863 1.65 1.9751 -.891 1.65.83134 -.18116 1.65.198317 -.8167 14 1.75 4.59855194 -.4359 1.75.96113 -.7594 1.75 1.96138 -.791 1.75.88589 -.18187 1.75.16318 -.8588 15 1.875 4.56816 -.3374 1.875.95151 -.781 1.875 1.956 -.7465 1.875.785855 -.1819 1.875.17373 -.8783 16 4.53888566 -.177.9467 -.6333 1.9493 -.679.76333 -.17613.91394 -.8754 17.15 4.51116444 -.784.15.934754 -.5361.15 1.93444 -.5885.15.741316 -.16977.15.5545 -.858 18.5 4.48518413 -.1913.5.9853 -.418.5 1.9787 -.476.5.795 -.1611.5.19816 -.854 19.375 4.46116789 -.17484.375.985 -.814.375 1.91134 -.3439.375.699943 -.156.375 -.155 -.745.5 4.4393139 -.15619.5.91938 -.179.5 1.916835 -.1936.5.681118 -.1389.5 -.4951 -.6576 1.65 4.41978917 -.13645.65.91779.3994.65 1.914415 -.74.65.663857 -.1388.65 -.873 -.5585.75 4.47381 -.11588.75.9188.1938.75 1.9147.15.75.64837 -.181.75 -.11471 -.4453 3.875 4.3884763 -.9477.875.9951.4761.875 1.91597.345.875.634846 -.913.875 -.1457 -.35 4 3 4.3764131 -.7341 3.9646.616 3 1.948.53976 3.63433 -.7344 3 -.1748 -.1865 5 3.15 4.36745 -.51 3.15.933566.7986 3.15 1.96995.741 3.15.61453 -.549 3.15 -.161 -.46 6 3.5 4.36713 -.3118 3.5.943544.9944 3.5 1.93671.946 3.5.67391 -.3597 3.5 -.719 -.1919 7 3.375 4.356817 -.19 3.375.955974.118691 3.375 1.94896.11483 3.375.6895 -.1695 3.375 -.596 -.1757 8 3.5 4.35544983.8418 3.5.9781.13759 3.5 1.96449.134567 3.5.6776.1843 3.5 -.79 -.16144 9 3.65 4.3565.654 3.65.987968.154843 3.65 1.9797.15351 3.65.616.11 3.65 -.931 -.14769 3 3.75 4.359815.4393 3.75 3.733.171155 3.75 1.99846.171388 3.75.635.37585 3.75 -.31157 -.13475 31 3.875 4.365167.57946 3.875 3.8718.18593 3.875.19884.187884 3.875.68.53951 3.875 -.3841 -.189 3 4 4.3744488.7846 4 3.51958.1989 4.43369.748 4.614964.68954 4 -.34377 -.1139 33 4.15 4.381368.81611 4.15 3.7681.9874 4.15.68713.1574 4.15.63583.8345 4.15 -.3578 -.135 34 4.5 4.39153.987 4.5 3.1355.18651 4.5.95681.6659 4.5.633876.939 4.5 -.3776 -.9645 35 4.375 4.4769.96159 4.375 3.13386.584 4.375.1413.3535 4.375.645614.1343 4.375 -.3881 -.9145 36 4.5 4.4147883.99749 4.5 3.1585.959 4.5.15349.416 4.5.65854.11766 4.5 -.3944 -.8866 37 4.65 4.475146.1818 4.65 3.187154.35 4.65.18369.45386 4.65.67388.115783 4.65 -.4533 -.885 38 4.75 4.43985365.99366 4.75 3.15967.938 4.75.143.466 4.75.686861.118388 4.75 -.41636 -.93 39 4.875 4.457435.95434 4.875 3.4464.571 4.875.4513.4591 4.875.71659.11858 4.875 -.4765 -.9496 4 5 4.464363.8914 5 3.7854.19536 5.75791.41417 5.716475.116189 5 -.4395 -.1 41 5.15 4.47534164.8493 5.15 3.395.1953 5.15.35968.3566 5.15.73999.111398 5.15 -.459 -.117 4 5.5 4.485433.69756 5.5 3.36665.95 5.5.33535.6345 5.5.74494.1418 5.5 -.4663 -.1453 43 5.375 4.494186.5781 5.375 3.351676.18718 5.375.363645.15398 5.375.757951.94753 5.375 -.48187 -.1395 44 5.5 4.5158.4685 5.5 3.37567.1757 5.5.39569.47 5.5.769795.8314 5.5 -.49931 -.15695 45 5.65 4.5659364.681 5.65 3.3966.155714 5.65.41587.187586 5.65.78188.69559 5.65 -.51893 -.17669 46 5.75 4.59945.979 5.75 3.41664.137758 5.75.439319.171177 5.75.788883.548 5.75 -.5411 -.19859 47 5.875 4.51116138 -.88 5.875 3.43383.11867 5.875.46716.153449 5.875.795659.373 5.875 -.56584 -.46 48 6 4.51164 -.69 6 3.448117.98749 6.479897.13469 6.834.1915 6 -.59364 -.481 49 6.15 4.5676169 -.4588 6.15 3.46461.7835 6.15.496733.1154 6.15.8718 -.6E-5 6.15 -.6466 -.7534 5 6.5 4.51675 -.6486 6.5 3.4749.57653 6.5.511134.9534 6.5.8714 -.199 6.5 -.6597 -.3388 51 6.375 4.499198 -.8354 6.375 3.477456.37113 6.375.5347.7539 6.375.84 -.44 6.375 -.6976 -.3335 5 6.5 4.4847743 -.1163 6.5 3.4895.17 6.5.5346.55539 6.5.795194 -.667 6.5 -.73875 -.364 17
Figure 1.14 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.9*(1+sin()) max 1 max 1.5 dela.5 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.9*(1+sin()) max 1 max 1.5 dela.5 18
Figure 1.15 6 5 4 3 1 4 6 8 1 1 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.9 3.15.15 1.7.5 1.15 4.98875 -.9898.15 3.1875.13815.15.1875.139515.15 1.875.5986.15.5 -.11.5 4.9763783 -.1756.5 3.36.164.5.36189.1919.5 1.168.49671.5.498597 -.49 3.375 4.969349 -.11561.375 3.5178.114857.375.5331.11919.375 1.437.3978.375.495786 -.3364 4.5 4.94848174 -.195.5 3.66139.1431.5.678.19359.5 1.741.311.5.491581 -.445 5.65 4.9331189 -.1946.65 3.79143.9395.65.8878.1315.65 1.31164.14.65.48619 -.549 6.75 4.9169336 -.1351.75 3.9884.84687.75.93417.94.75 1.33795.166.75.479156 -.647 7.875 4.95416 -.13949.875 3.1147.7645.875.149.84678.875 1.35377.5116.875.47168 -.7374 8 1 4.8861799 -.1481 1 3.1116.69334 1.11557.78354 1 1.3617 -.146 1.461851 -.8189 9 1.15 4.8647671 -.14489 1.15 3.119693.63435 1.15.1531.7318 1.15 1.35834 -.696 1.15.451614 -.894 1 1.5 4.84665597 -.14568 1.5 3.176.58835 1.5.134449.6946 1.5 1.34965 -.116 1.5.44484 -.957 11 1.375 4.8844593 -.14515 1.375 3.134976.5559 1.375.14314.6665 1.375 1.33557 -.143 1.375.486 -.9991 1 1.5 4.81358 -.1437 1.5 3.14195.53743 1.5.151433.65365 1.5 1.3177 -.16 1.5.416111 -.1348 13 1.65 4.793937 -.147 1.65 3.148643.5333 1.65.15963.65497 1.65 1.9769 -.1633 1.65.43177 -.1573 14 1.75 4.77488446 -.13556 1.75 3.15536.5464 1.75.16779.677 1.75 1.778 -.156 1.75.38996 -.1665 15 1.875 4.7579391 -.198 1.875 3.1689.56598 1.875.176169.69938 1.875 1.58 -.18 1.875.37669 -.16 16 4.7417141 -.185 3.169164.65.184911.7419 1.41 -.895.36335 -.1445 17.15 4.76358 -.1148.15 3.176695.65155.15.194185.7973.15 1.311 -.378.15.3596 -.1138 18.5 4.71875 -.1575.5 3.18484.7113.5.4151.8647.5 1.69.668.5.33763 -.976 19.375 4.6987951 -.958.375 3.193741.78317.375.1496.94318.375 1.963.133.375.35491 -.9157.5 4.68681316 -.8515.5 3.3531.86339.5.675.13151.5 1.451.194.5.31444 -.85 1.65 4.67616958 -.7388.65 3.1433.95137.65.39644.11836.65 1.669.8714.65.33419 -.7745.75 4.66693441 -.617.75 3.615.14555.75.53748.135.75 1.318.394.75.93738 -.695 3.875 4.659167 -.519.875 3.3985.11448.875.69151.134158.875 1.3513.5455.875.8516 -.5994 4 3 4.65889 -.381 3 3.53588.14583 3.8591.145466 3 1.4143.6 3.77614 -.56 5 3.15 4.6481676 -.67 3.15 3.69161.13484 3.15.3415.15697 3.15 1.495.74314 3.15.71331 -.417 6 3.5 4.64486786 -.148 3.5 3.8617.1455 3.5.3376.168495 3.5 1.58494.866 3.5.6631 -.984 7 3.375 4.643876 -.9 3.375 3.34145.154951 3.375.344788.179853 3.375 1.693.98943 3.375.658 -.1943 8 3.5 4.6478.79 3.5 3.33513.164444 3.5.36769.19866 3.5 1.8169.11116 3.5.615 -.911 9 3.65 4.6437185.181 3.65 3.34469.173334 3.65.39118.1358 3.65 1.95578.1936 3.65.5914.951 3 3.75 4.6459635.734 3.75 3.365736.181461 3.75.41697.1116 3.75 1.11945.1346 3.75.59133.1581 31 3.875 4.64936665.35455 3.875 3.388418.188675 3.875.4469.115 3.875 1.1778.14491 3.875.6455.1964 3 4 4.65379857.4558 4 3.413.19484 4.477.877 4 1.145843.154763 4.698.799 33 4.15 4.65911835.48443 4.15 3.436358.199845 4.15.498716.34915 4.15 1.165188.163646 4.15.66399.35356 34 4.5 4.66517379.533 4.5 3.461338.358 4.5.5881.4518 4.5 1.185644.171448 4.5.7819.41779 35 4.375 4.671853.5657 4.375 3.486786.5976 4.375.558145.44791 4.375 1.775.1786 4.375.7641.4786 36 4.5 4.67883465.5885 4.5 3.51533.6969 4.5.588744.47663 4.5 1.933.183394 4.5.8197.51184 37 4.65 4.686956.58496 4.65 3.53844.655 4.65.6197.4983 4.65 1.557.18738 4.65.8835.53996 38 4.75 4.693477.57495 4.75 3.564.4634 4.75.65838.496 4.75 1.75679.189978 4.75.9574.55468 39 4.875 4.759413.5518 4.875 3.589799.136 4.875.681966.47488 4.875 1.9947.191158 4.875.38.55566 4 5 4.774865.51384 5 3.61496.196578 5.719.4449 5 1.3331.199 5.38954.5476 41 5.15 4.7139565.4639 5.15 3.639535.1955 5.15.743463.476 5.15 1.347186.18988 5.15.315738.5167 4 5.5 4.719746.41 5.5 3.663348.183169 5.5.773473.34313 5.5 1.37847.18635 5.5.3189.4759 43 5.375 4.747319.3978 5.375 3.68644.174668 5.375.876.788 5.375 1.394135.1838 5.375.38138.476 44 5.5 4.788544.4787 5.5 3.7877.16511 5.5.831173.1911 5.5 1.41689.176574 5.5.3334.35736 45 5.65 4.7319577.15787 5.65 3.78717.15466 5.65.85856.999 5.65 1.43896.1719 5.65.337889.863 46 5.75 4.7339613.618 5.75 3.7485.14344 5.75.8848.1998 5.75 1.4614.16499 5.75.341397.1936 47 5.875 4.7346914 -.43 5.875 3.76598.131617 5.875.99778.1898 5.875 1.4857.15415 5.875.343817.976 48 6 4.7341885 -.145 6 3.7843.11936 6.93344.177634 6 1.499796.145131 6.34537 -.61 49 6.15 4.7337395 -.516 6.15 3.797353.1685 6.15.95568.165874 6.15 1.517937.135599 6.15.34496 -.1159 5 6.5 4.79866 -.3579 6.5 3.8179.9468 6.5.976343.15391 6.5 1.534887.1578 6.5.343513 -.33 51 6.375 4.747546 -.463 6.375 3.8493.81789 6.375.995581.14194 6.375 1.5563.115694 6.375.34634 -.3477 5 6.5 4.71897549 -.5631 6.5 3.83716.69594 6.5 3.133.1399 6.5 1.56565.15675 6.5.33688 -.4663 19
Figure 1.16 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.4*(1+sin()) max 1 max 1.5 dela.5 1 1 = 4.918 = 1.37.**(1-/5)-.4*(1+sin()) d/d= -.54 Equaions Runge Kua 4 Draw s max 1 max 1.5 dela.5
Figure 1.17 6 4-4 6 8 1 1-4 -6-8 4 6 DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.4 3 1 -.8.5 -.15 1.15 4.97 -.6396.15 3 -.99.15 -.99.15.99 -.11113.15.4815 -.1894.5 4.93756 -.8694.5.9966 -.593.5 1.9966 -.5953.5.97619 -.147.5.458383 -.161 3.375 4.9113817 -.385.375.988856 -.8746.375 1.988819 -.8836.375.95836 -.1798.375.43137 -.497 4.5 4.865767 -.3833.5.97793 -.114.5 1.977774 -.11597.5.93674 -.8.5.436 -.814 5.65 4.8153115 -.346.65.963648 -.139.65 1.96378 -.14195.65.91135 -.3138.65.36558 -.3174 6.75 4.77887 -.361.75.94671 -.16156.75 1.945535 -.16589.75.8843 -.585.75.35965 -.3465 7.875 4.733187 -.37368.875.9676 -.18147.875 1.94798 -.18744.875.85148 -.839.875.83134 -.3779 8 1 4.6864181 -.38317 1.9339 -.19846 1 1.91368 -.69 1.81476 -.3555 1.36785 -.39684 9 1.15 4.6385185 -.38948 1.15.878584 -.18 1.15 1.87558 -.14 1.15.776567 -.3535 1.15.18718 -.451 1 1.5 4.58983744 -.3945 1.5.8549 -.71 1.5 1.847814 -.3477 1.5.735898 -.344 1.5.134616 -.44156 11 1.375 4.54788 -.391 1.375.841 -.96 1.375 1.818468 -.4399 1.375.693118 -.3561 1.375.794 -.45978 1 1.5 4.491788 -.3889 1.5.79558 -.389 1.5 1.787969 -.4968 1.5.648617 -.3665 1.5.1949 -.4753 13 1.65 4.443699 -.387 1.65.766396 -.349 1.65 1.756759 -.5174 1.65.685 -.3736 1.65 -.3743 -.48719 14 1.75 4.39568187 -.3699 1.75.737335 -.841 1.75 1.7591 -.516 1.75.5561 -.37731 1.75 -.9833 -.4961 15 1.875 4.3494444 -.3558 1.875.78784 -.7 1.875 1.694 -.4496 1.875.58939 -.37755 1.875 -.1635 -.58 16 4.3496945 -.33855.681194 -.954 1.6634 -.363.461744 -.37441 -.311 -.5485 17.15 4.66595 -.31835.15.6551 -.1954.15 1.633871 -.48.15.414943 -.36798.15 -.86 -.546 18.5 4.85667 -.9547.5.6361 -.1774.5 1.65861 -.87.5.368946 -.35839.5 -.3493 -.5148 19.375 4.185933 -.718.375.68443 -.15695.375 1.579771 -.1936.375.34147 -.34586.375 -.41198 -.49567.5 4.151516 -.48.5.58884 -.13395.5 1.555977 -.1698.5.8915 -.3361.5 -.47394 -.48741 1.65 4.1179888 -.1375.65.578 -.1875.65 1.534817 -.1458.65.39589 -.319.65 -.53486 -.47696.75 4.95811 -.18337.75.558487 -.8174.75 1.516591 -.18.75.474 -.9311.75 -.59448 -.46463 3.875 4.715888 -.159.875.5487 -.533.875 1.51556 -.931.875.163835 -.7153.875 -.6556 -.4577 4 3 4.531478 -.136 3.54165 -.394 3 1.489918 -.6468 3.19893 -.4857 3 -.7891 -.43575 5 3.15 4.38115 -.8861 3.15.538613.5958 3.15 1.481833 -.3545 3.15.9883 -.461 3.15 -.76338 -.41997 6 3.5 4.756 -.5731 3.5.539357.3595 3.5 1.4774 -.586 3.5.7747 -.8 3.5 -.81587 -.4383 7 3.375 4.19864 -.689 3.375.543845.65434 3.375 1.476669.36 3.375.45736 -.17543 3.375 -.86635 -.38778 8 3.5 4.165133.197 3.5.555.948 3.5 1.4796.5541 3.5.388 -.1517 3.5 -.91483 -.375 9 3.65 4.1677594.957 3.65.563785.11389 3.65 1.48619.8439 3.65.494 -.1746 3.65 -.96136 -.35769 3 3.75 4.46686.547 3.75.578958.14695 3.75 1.49645.16874 3.75 -.111 -.153 3.75.67 -.34453 31 3.875 4.7346.77351 3.875.59734.1778 3.875 1.5964.1314 3.875 -.414 -.8419 3.875.4913 -.333 3 4 4.3697345.97141 4.61869.1917 4 1.563.153688 4 -.3466 -.6535 4.978 -.341 33 4.15 4.491169.113783 4.15.6449.896 4.15 1.54543.173311 4.15 -.483 -.4887 4.15.1313 -.31768 34 4.5 4.63339.1736 4.5.668613.366 4.5 1.56696.189971 4.5 -.4894 -.359 4.5.1711 -.3145 35 4.375 4.791854.13671 4.375.69657.3494 4.375 1.59653.3393 4.375 -.5333 -.431 4.375.19 -.31418 36 4.5 4.963755.14679 4.5.75936.4565 4.5 1.61677.13354 4.5 -.5636 -.1679 4.5.4956 -.31776 37 4.65 4.1141446.144866 4.65.75656.46457 4.65 1.64746.19689 4.65 -.5846 -.175 4.65.898 -.356 38 4.75 4.1357.14361 4.75.78764.46534 4.75 1.677.88 4.75 -.66 -.133 4.75.3994 -.33691 39 4.875 4.151588.137913 4.875.81788.479 4.875 1.697993.115 4.875 -.616 -.1564 4.875.376 -.3588 4 5 4.1673974.18933 5.8489.3591 5 1.7563.16156 5 -.6355 -.73 5.41616 -.37331 41 5.15 4.1835146.11649 5.15.877641.4154 5.15 1.75651.7517 5.15 -.6639 -.336 5.15.4683 -.3983 4 5.5 4.198757.186 5.5.9566.956 5.5 1.778591.19537 5.5 -.759 -.4817 5.5.516 -.479 43 5.375 4.16761.8157 5.375.931855.191753 5.375 1.837.179781 5.375 -.7661 -.6634 5.375.5661 -.4611 44 5.5 4.9456.6864 5.5.95584.17119 5.5 1.85479.161131 5.5 -.8491 -.8794 5.5.6387 -.59 45 5.65 4.855355.3788 5.65.977.147695 5.65 1.84561.139675 5.65 -.959 -.1174 5.65.68648 -.5446 46 5.75 4.331453.1186 5.75.995664.116 5.75 1.8638.11576 5.75 -.1999 -.145 5.75.7545 -.595 47 5.875 4.346974 -.151 5.875 3.1933.9485 5.875 1.87755.89769 5.875 -.1755 -.1789 5.875.885 -.64418 48 6 4.3857 -.434 6 3.787.6617 6 1.888771.6116 6 -.1489 -.361 6.994 -.753 49 6.15 4.74467 -.7155 6.15 3.315.3656 6.15 1.896536.334 6.15 -.17437 -.388 6.15.99661 -.7697 5 6.5 4.184848 -.116 6.5 3.35618.6488 6.5 1.969.3596 6.5 -.415 -.7453 6.5 -.9173 -.854 51 6.375 4.595998 -.184 6.375 3.3649 -.351 6.375 1.9114 -.635 6.375 -.3847 -.31197 6.375 -.19491 -.89369 5 6.5 4.1899791 -.15586 6.5 3.3349 -.531 6.5 1.897846 -.5613 6.5 -.7747 -.35 6.5 -.366 -.96577 1
Figure 1.18 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.5*(1+sin()) max 1 max 1.5 dela.5 1 1 Equaions Runge Kua 4 Draw s.**(1-/5)-.5*(1+sin()) max 1 max 1.5 dela.5 Figure 1.19
1 5 4 6 8 1 1-5 5 - DELTA T.15 k _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) _k _k f(_k,_k) 5 -.5 3 -.1 -.1 1 -.9.5 -.16 1.15 4.96875 -.7496.15.99875 -.411.15 1.99875 -.41.15.98875 -.15.15.48 -.19438.5 4.934388 -.989.5.99361 -.716.5 1.993598 -.711.5.973435 -.1557.5.4557 -.9 3.375 4.8971787 -.314.375.984661 -.196.375 1.984584 -.119.375.95451 -.18717.375.4775 -.6345 4.5 4.85684338 -.344.5.974 -.1877.5 1.97181 -.131.5.93655 -.1837.5.394144 -.974 5.65 4.8148778 -.3647.65.955944 -.15459.65 1.955433 -.15814.65.93359 -.484.65.356988 -.3997 6.75 4.769848 -.37635.75.9366 -.1784.75 1.935666 -.18315.75.8738 -.7638.75.31574 -.3615 7.875 4.719848 -.38937.875.914366 -.19875.875 1.9177 -.568.875.837781 -.34.875.7586 -.397 8 1 4.673313 -.3993 1.8895 -.1644 1 1.8876 -.54 1.79998 -.3597 1.1748 -.41798 9 1.15 4.63458 -.459 1.15.86467 -.38 1.15 1.858888 -.41 1.15.75934 -.34678 1.15.1695 -.448 1 1.5 4.576657 -.498 1.5.833614 -.417 1.5 1.88637 -.558 1.5.715887 -.36457 1.5.114148 -.46494 11 1.375 4.515519 -.4869 1.375.834 -.4891 1.375 1.79677 -.651 1.375.67316 -.37913 1.375.5631 -.48414 1 1.5 4.4743946 -.4468 1.5.7789 -.534 1.5 1.7636 -.717 1.5.694 -.3931 1.5 -.449 -.57 13 1.65 4.41985459 -.3977 1.65.74746 -.5195 1.65 1.79718 -.7337 1.65.574135 -.39799 1.65 -.67 -.513 14 1.75 4.3717 -.38591 1.75.795 -.4775 1.75 1.695548 -.7188 1.75.54386 -.41 1.75 -.13117 -.59 15 1.875 4.319838 -.37131 1.875.67884 -.3979 1.875 1.66156 -.6664 1.875.47411 -.469 1.875 -.19654 -.5937 16 4.7556978 -.35343.64831 -.8 1.683 -.577.43785 -.39975 -.671 -.5363 17.15 4.31391 -.3349.15.619784 -.1315.15 1.59617 -.457.15.373816 -.39341.15 -.399 -.5377 18.5 4.1898993 -.3874.5.59314 -.19487.5 1.565358 -.946.5.34641 -.38381.5 -.39588 -.5996 19.375 4.151375 -.848.375.56878 -.17361.375 1.536676 -.154.375.76665 -.37115.375 -.4613 -.5439.5 4.1159693 -.547.5.54781 -.14971.5 1.51358 -.18879.5.371 -.35568.5 -.5768 -.5169 1.65 4.8416858 -.386.65.58367 -.1351.65 1.486759 -.16455.65.185811 -.3377.65 -.59 -.5595.75 4.561856 -.198.75.5198 -.954.75 1.466191 -.13817.75.143598 -.3175.75 -.65546 -.49369 3.875 4.31515 -.15976.875.51 -.6586.875 1.4489 -.115.875.1398 -.9551.875 -.71717 -.47987 4 3 4.1184 -.1675 3.49768 -.358 3 1.435164 -.864 3.66969 -.77 3 -.77715 -.46487 5 3.15 3.9963366 -.9371 3.15.488358 -.415 3.15 1.4584 -.537 3.15.3961 -.476 3.15 -.8356 -.44911 6 3.5 3.984698 -.611 3.5.487838.743 3.5 1.418788 -.1971 3.5.11 -.55 3.5 -.8914 -.433 7 3.375 3.97698355 -.944 3.375.49119.578 3.375 1.41634.1849 3.375 -.581 -.19736 3.375 -.94553 -.4174 8 3.5 3.973348.871 3.5.498446.87696 3.5 1.417681.4839 3.5 -.548 -.175 3.5 -.99766 -.4165 9 3.65 3.9734193.936 3.65.5948.116196 3.65 1.4785.69784 3.65 -.74 -.1484 3.65.4787 -.3879 3 3.75 3.9778317.55619 3.75.53933.14867 3.75 1.43158.973 3.75 -.959 -.1556 3.75.968 -.37444 31 3.875 3.9843559.7957 3.875.541791.16781 3.875 1.443661.1717 3.875 -.169 -.1436 3.875.1438 -.36353 3 4 3.9939468.9996 4.5671.18943 4 1.4591.145853 4 -.11933 -.854 4.1885 -.3551 33 4.15 4.643343.11734 4.15.58633.78 4.15 1.47733.16656 4.15 -.1999 -.6858 4.15.39 -.3498 34 4.5 4.19893.131198 4.5.6137.343 4.5 1.49815.183588 4.5 -.13856 -.5473 4.5.7656 -.34675 35 4.375 4.3749864.141349 4.375.641.35119 4.375 1.5963.197565 4.375 -.1454 -.44 4.375.3199 -.34776 36 4.5 4.551673.147641 4.5.6696.433 4.5 1.545659.795 4.5 -.159 -.3671 4.5.36337 -.3564 37 4.65 4.73641.149995 4.65.76.47446 4.65 1.571653.14573 4.65 -.15549 -.33 4.65.4745 -.36168 38 4.75 4.93717.148397 4.75.73937.4769 4.75 1.598474.17313 4.75 -.1596 -.331 4.75.4566 -.3751 39 4.875 4.119138.14899 4.875.761898.43958 4.875 1.65639.161 4.875 -.16376 -.371 4.875.49955 -.39316 4 5 4.1878379.133614 5.79393.36311 5 1.65654.1111 5 -.1684 -.458 5.5487 -.41595 41 5.15 4.14548549.1714 5.15.8193.4874 5.15 1.6793.6 5.15 -.1743 -.57 5.15.669 -.44361 4 5.5 4.1657478.14433 5.5.8541.983 5.5 1.7488.18947 5.5 -.18116 -.781 5.5.65614 -.4761 43 5.375 4.1736894.8556 5.375.8767.191434 5.375 1.77964.17356 5.375 -.196 -.94 5.375.71567 -.51378 44 5.5 4.18469.6916 5.5.9199.169979 5.5 1.74961.15386 5.5 -.181 -.11561 5.5.77989 -.55631 45 5.65 4.191538.38389 5.65.91446.145816 5.65 1.768853.131538 5.65 -.166 -.14 5.65.84943 -.6378 46 5.75 4.196941.11888 5.75.939673.119337 5.75 1.78595.16638 5.75 -.344 -.17193 5.75.949 -.6561 47 5.875 4.1984998 -.1615 5.875.954591.997 5.875 1.79865.79559 5.875 -.5553 -.448 5.875 -.69 -.7136 48 6 4.19639151 -.456 6.96596.61169 6 1.8857.5731 6 -.819 -.395 6 -.967 -.7751 49 6.15 4.197345 -.7496 6.15.97368.349 6.15 1.814911.68 6.15 -.3113 -.7669 6.15 -.1996 -.84157 5 6.5 4.1813643 -.1478 6.5.97749 -.8 6.5 1.817487 -.134 6.5 -.3456 -.31561 6.5 -.9816 -.916 51 6.375 4.168663 -.1345 6.375.97736 -.33 6.375 1.816195 -.4163 6.375 -.3857 -.35587 6.375 -.413 -.9881 5 6.5 4.15148548 -.1688 6.5.9733 -.674 6.5 1.8199 -.777 6.5 -.4955 -.3977 6.5 -.53575.6813 3