{"id":76,"date":"2018-07-11T11:49:44","date_gmt":"2018-07-11T11:49:44","guid":{"rendered":"http:\/\/itp15.epgpbooks.inflibnet.ac.in\/?post_type=chapter&p=76"},"modified":"2019-05-15T08:40:56","modified_gmt":"2019-05-15T08:40:56","slug":"iterative-methods-to-solve-equation-f-x-0-method-of-false-position","status":"publish","type":"chapter","link":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/chapter\/iterative-methods-to-solve-equation-f-x-0-method-of-false-position\/","title":{"rendered":"Iterative Methods to solve equation f (x) = 0: Method of False Position"},"content":{"raw":"
\"epgp<\/a>\r\n<\/span><\/div>\r\n
\r\n\r\n \r\n

3.1\u00a0 Introduction\u00a0<\/strong><\/p>\r\n \r\n

In \u00a0our \u00a0previous \u00a0lecture \u00a0(module), \u00a0we \u00a0learned \u00a0about \u00a0nature \u00a0of \u00a0iterative \u00a0methods, classification \u00a0of \u00a0iterative \u00a0methods \u00a0(bracketing \u00a0or \u00a0open) \u00a0to \u00a0find \u00a0roots \u00a0of \u00a0an \u00a0equation\u00a0f <\/em>(<\/span>x<\/em>) = 0<\/span>. We began with Bisection method to find roots of <\/span>f <\/em>(<\/span>x<\/em>) = 0 . It was observed\u00a0<\/span>that through Bisection method is a very simple, steady and reliable method, it suffers from drawback of being slow. It would be advisable to seek its modification to make it faster. The regular falsi method or method of false position is an attempt in that direction. This also is closed method. As in, Bisection method, at the beginning of iterative process\u00a0<\/span>two end points <\/span>a <\/em>and <\/span>b <\/em>are found so that\u00a0f <\/em>(a<\/em>)\u00b7 f <\/em>(b<\/em>) < 0 to ensure that at least one root lies\u00a0<\/span>between <\/span>a <\/em>and <\/span>b<\/em>.<\/span><\/p>\r\n \r\n

It was also pointed out in disadvantages of Bisection method that it does not take into account of magnitude of function values at end points to generate estimates for the next iteration. It just considers sign of function values. This makes it slow.<\/span><\/p>\r\n \r\n

Thus, in method of false position, the function values at the end points <\/span>a <\/em>and <\/span>b <\/em>are taken into account to generate the estimates for next iteration.<\/span><\/p>\r\n \r\n

Method of false position differs with Bisection method only in the \u00a0methodology \u00a0of finding <\/span>c <\/em>. In Bisection Method, interval was subdivided into equal parts, \u00a0<\/span>c <\/em>was mid point given by (<\/span>a <\/em>+ <\/span>b<\/em>) 2.<\/span><\/p>\r\n \r\n

No account was taken of function values at <\/span>a <\/em>and <\/span>b <\/em>. As at the root, function value is zero, the point, where function value is smaller is expected to be nearer to the root. It can be visualized also from the graph below.<\/span><\/p>\r\n\r\n<\/div>\r\n

\r\n\r\n\"\"\r\n\r\n \r\n\r\n3.2 Iterative Process\u00a0<\/strong>\r\n\r\n \r\n\r\nSo, in regular falsi method, to take into consideration the function values at \u00a0a <\/em>and b <\/em>,\u00a0straight line is drawn, joining\u00a0<\/span>(<\/span>a<\/em>, <\/span>f <\/em>(<\/span>a<\/em>)) and (<\/span>b<\/em>, <\/span>f <\/em>(<\/span>b<\/em>)) .<\/span>\r\n\r\n<\/div>\r\n
\r\n\r\n \r\n\r\nThe point, where it cuts \u00a0X \u00a0<\/em>axis is, the new estimate of the root. Mathematically, estimate of formula for c <\/em>can be derived as follows:\r\n\r\n<\/div>\r\n

Equation of straight line joining\u00a0<\/span>(<\/span>a<\/em>, <\/span>f <\/em>(<\/span>a<\/em>)) and (<\/span>b<\/em>, <\/span>f <\/em>(<\/span>b<\/em>))\u00a0<\/span>is given by :<\/span><\/p>\r\n\r\n

\r\n\r\n\"\"\r\n\r\nIf \u00a0f <\/em>(c<\/em>) = 0,c <\/em>is the root otherwise if f <\/em>(a<\/em>)\u00b7 f <\/em>(c<\/em>) > 0 , \u00a0c <\/em>takes the role of \u00a0a <\/em>otherwise \u00a0c \u00a0<\/em>takes the role of b<\/em>.\r\n\r\n \r\n\r\n3.3\u00a0 Graphical illustration of the Iterative Process\u00a0<\/strong>\r\n\r\n \r\n\r\nCase \u2013 1:f(<\/em><\/strong>x)<\/em><\/strong>\r\n\r\n<\/div>\r\n
\r\n\r\n\"\"\r\n\r\n<\/div>\r\n\"\"\r\n\r\n3.4\u00a0 Different Stopping Criterion\u00a0<\/strong>\r\n\r\n \r\n

Now, once we have observed the iterative process, let us focus on different stopping criterions. The sequence of <\/span>c <\/em>k\u2019s approaches the root. That is why, in case of the Bisection section algorithm, iterative process was stopped, when<\/span><\/p>\r\n \r\n

| <\/strong>c <\/em>k <\/strong>\u2013 <\/strong>c <\/em>k-1<\/strong>| <<\/strong>\u00ce <\/span>\u2026 (3)<\/strong><\/p>\r\n\r\n

\r\n\r\n \r\n

It is called x \u2013 tolerance criterion (X-TOL<\/strong>). The recent most c <\/em>k is the estimate of the root. There is another stopping criterion called function tolerance (F<\/strong>\u2013TOL)<\/strong>. At root, function value is zero, so if estimate is quite near the root then function value would be small. Hence, many times, one would like to stop when<\/p>\r\n \r\n\r\n|f ( <\/strong>c <\/em>k<\/strong>)| < <\/strong>d \u00a0\u2026 (4)<\/strong>\r\n\r\n \r\n\r\nwhere d is some small positive constant like \u20ac.\r\n\r\n \r\n

You would like to ask: <\/span>which one is better or what basis decision should be made to choose the particular criterion<\/em>? Well, it depends on the method, nature of function near the root, objective of finding the root etc. Many times, a combination of different criterions are also applied.<\/span><\/p>\r\n \r\n

If slope of the curve is small near the root, curve is almost horizontal, and then function tolerance may not be appropriate stopping criterion, because curve is rising slowly, function values in neighborhood of the root are going to be small, so even if estimate <\/span>c <\/em>is not sufficiently near the root, one may stop. This is visible from the following graph:<\/span><\/p>\r\n\r\n<\/div>\r\n

\r\n\r\n\"\"\r\n

On the other hand, if it is expected that slope is high near the root or sequence\u00a0 of approximations c <\/em>k\u2019s may converge to the root, like in case of almost vertical graph, then function tolerance ensures that estimate is good approximation to the root.<\/p>\r\n\r\n<\/div>\r\n

\r\n\r\n\"\"\r\n\r\nIn simple words, for <\/span>f <\/em>(<\/span>x<\/em>) , if\u00a0<\/span>f <\/em>\u00a2(<\/span>x<\/em>)\u00a0\u00a0<\/span>is high near the root, FTOL would be better to use as\u00a0<\/span>stopping criterion.\u00a0<\/span>\r\n\r\n<\/div>\r\n\"\"\r\n

\"\"<\/p>\r\n

3.5 Algorithm for Regular False Method\u00a0<\/span><\/strong><\/p>\r\n \r\n

\"\"<\/div>\r\n \r\n

Let us solve examples to find root of an equation using method of false position. For comparison purpose, we shall find the roots of the same equations, used for illustrating Bisection Method.<\/span><\/p>\r\n\r\n

\r\n\r\n \r\n\r\nExample 1:\r\n\r\n \r\n\r\nf <\/em>(x<\/em>) = 3x <\/em>- cos x <\/em>-1\r\n\r\n \r\n\r\nTake a <\/em>= 0, f <\/em>(a<\/em>) = -2.0, c<\/em>0 = a<\/em>\r\n\r\n \r\n\r\nb <\/em>= 1, f <\/em>(b<\/em>) = 1.459698\r\n\r\nf <\/em>(a<\/em>). f <\/em>(b<\/em>) < 0\r\n\r\n \r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
k<\/strong><\/td>\r\na<\/strong><\/td>\r\nf(a)<\/strong><\/td>\r\nb<\/strong><\/td>\r\nf(b)<\/strong><\/td>\r\nc<\/strong>k<\/strong><\/td>\r\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\r\nc<\/strong>k<\/strong>- c<\/strong>k-1<\/strong><\/td>\r\n<\/tr>\r\n
1<\/td>\r\n0.000000<\/td>\r\n-2.000000<\/td>\r\n1.000000<\/td>\r\n1.459698<\/td>\r\n0.578085<\/td>\r\n-0.103255<\/td>\r\n0.578085<\/td>\r\n<\/tr>\r\n
2<\/td>\r\n0.578085<\/td>\r\n-0.103255<\/td>\r\n1.000000<\/td>\r\n1.459698<\/td>\r\n0.605959<\/td>\r\n-0.004081<\/td>\r\n0.027873<\/td>\r\n<\/tr>\r\n
3<\/td>\r\n0.605959<\/td>\r\n-0.004081<\/td>\r\n1.000000<\/td>\r\n1.459698<\/td>\r\n0.607057<\/td>\r\n-0.000159<\/td>\r\n0.001099<\/td>\r\n<\/tr>\r\n
4<\/td>\r\n0.607057<\/td>\r\n-0.000159<\/td>\r\n1.000000<\/td>\r\n1.459698<\/td>\r\n0.607100<\/td>\r\n-0.000006<\/td>\r\n0.000043<\/td>\r\n<\/tr>\r\n
5<\/td>\r\n0.607100<\/td>\r\n-0.000006<\/td>\r\n1.000000<\/td>\r\n1.459698<\/td>\r\n0.607102<\/td>\r\n-0.000000<\/td>\r\n0.000002<\/td>\r\n<\/tr>\r\n
6<\/td>\r\n0.607102<\/td>\r\n-0.000000<\/td>\r\n1.000000<\/td>\r\n1.459698<\/td>\r\n0.607102<\/td>\r\n-0.000000<\/td>\r\n0.000000<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n \r\n\r\nExample 2:\r\n\r\n \r\n\r\nf <\/em>(x<\/em>) = x<\/em>3 + x <\/em>-1\r\n\r\n \r\n\r\nTake a <\/em>= 0, f <\/em>(a<\/em>) = -1,c0 = a <\/em>= 0\r\n\r\n \r\n\r\nb <\/em>=1, f <\/em>(b<\/em>) =1\r\n\r\n \r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
k<\/strong><\/td>\r\na<\/strong><\/td>\r\nf(a)<\/strong><\/td>\r\nb<\/strong><\/td>\r\nf(b)<\/strong><\/td>\r\nc<\/strong>k<\/strong><\/td>\r\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\r\nc<\/strong>k<\/strong>- c<\/strong>k-1<\/strong><\/td>\r\n<\/tr>\r\n
1<\/td>\r\n0.000000<\/td>\r\n-1.000000<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.500000<\/td>\r\n-0.375000<\/td>\r\n0.500000<\/td>\r\n<\/tr>\r\n
2<\/td>\r\n0.500000<\/td>\r\n-0.375000<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.636364<\/td>\r\n-0.105935<\/td>\r\n0.136364<\/td>\r\n<\/tr>\r\n
3<\/td>\r\n0.636364<\/td>\r\n-0.105935<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.671196<\/td>\r\n-0.026428<\/td>\r\n0.034832<\/td>\r\n<\/tr>\r\n
4<\/td>\r\n0.671196<\/td>\r\n-0.026428<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.679662<\/td>\r\n-0.006376<\/td>\r\n0.008466<\/td>\r\n<\/tr>\r\n
5<\/td>\r\n0.679662<\/td>\r\n-0.006376<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.681691<\/td>\r\n-0.001525<\/td>\r\n0.002029<\/td>\r\n<\/tr>\r\n
6<\/td>\r\n0.681691<\/td>\r\n-0.001525<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.682176<\/td>\r\n-0.000364<\/td>\r\n0.000485<\/td>\r\n<\/tr>\r\n
7<\/td>\r\n0.682176<\/td>\r\n-0.000364<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.682292<\/td>\r\n-0.000087<\/td>\r\n0.000116<\/td>\r\n<\/tr>\r\n
8<\/td>\r\n0.682292<\/td>\r\n-0.000087<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.682319<\/td>\r\n-0.000021<\/td>\r\n0.000028<\/td>\r\n<\/tr>\r\n
9<\/td>\r\n0.682319<\/td>\r\n-0.000021<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.682326<\/td>\r\n-0.000005<\/td>\r\n0.000007<\/td>\r\n<\/tr>\r\n
10<\/td>\r\n0.682326<\/td>\r\n-0.000005<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.682327<\/td>\r\n-0.000001<\/td>\r\n0.000002<\/td>\r\n<\/tr>\r\n
11<\/td>\r\n0.682327<\/td>\r\n-0.000001<\/td>\r\n1.000000<\/td>\r\n1.000000<\/td>\r\n0.682328<\/td>\r\n-0.000000<\/td>\r\n0.000000<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n \r\n\r\nProgramming Code in C:\r\n\r\nThe above output was obtained by running the following C program.\r\n\r\n \r\n\r\n\/\/Regular falsi method\r\n\r\n \r\n\r\n#include<stdio.h>\r\n\r\n#include<conio.h>\r\n\r\n#include<math.h>\r\n\r\n#define fnx(x) pow (X, 3)-X+1\r\n\r\n \r\n\r\nvoid main()\r\n\r\n{\r\n\r\nfloat a, b, c,fa, fb, fc, error,co; int i=0;\r\n\r\nint count=0;\r\n\r\nfloat err = pow(10.0,-6.0)\/2;\r\n\r\n \r\n\r\nprintf(\"\\nEnter initial value of a = \"); scanf(\"%f\", &a);\r\n\r\nfa = fnx(a);\r\n\r\nprintf(\"\\n value of f(a) = %f\", fa);\r\n\r\n \r\n\r\nprintf(\"\\nEnter initial value of b = \"); scanf(\"%f\", &b);\r\n\r\nfb = fnx(b);\r\n\r\nprintf(\"\\n value of f(b) = %f\", fb); co=a;\r\n\r\nprintf(\"\\nNumber\\t A\\t F(A)\\t B\\t F(B)\\t C\\t F(C)\\t Error\");\r\n\r\n \r\n\r\ndo\r\n\r\n<\/div>\r\n
\r\n\r\n{\r\n\r\ni++;\r\n\r\nc = (a*fb - b*fa) \/ (fb - fa);\r\n\r\nfc = fnx(c); if(fc==0)\r\n\r\n{\r\n\r\nprintf(\"\\ndesired root is %f\", c);\r\n\r\n\/\/exit(0);\r\n\r\n}\r\n\r\nerror=fabs(co-c); co=c;\r\n\r\nprintf(\"\\n\\n%d\\t %f\\t %f\\t %f\\t %f\\t %f\\t %f\\t %f\", i, a, fa, b, fb, c, fc, error); if(fa*fc >0)\r\n\r\n{\r\n\r\n<\/div>\r\n

}<\/span><\/p>\r\n\r\n

\r\n\r\nelse\r\n\r\n{\r\n\r\n \r\n\r\n}\r\n\r\na = c; fa= fc;\r\n\r\n \r\n\r\nb = c; fb = fc;\r\n\r\n<\/div>\r\n
\r\n\r\n}while(error >err);\r\n\r\n \r\n\r\nprintf(\"\\n\\nThe root of the equation = %f\", c);\r\n\r\n \r\n\r\ngetch();\r\n\r\n}\r\n\r\n \r\n\r\n3.7 Observations\u00a0<\/strong>\r\n\r\n \r\n

1)\u00a0 From the tables, it is clear that same desired accuracy is obtained in less number of iterations. Question is, is this the\u00a0 case for these two specific equations, or in general we can expect False Position Method to perform better as compared to Bisection Method. Well, in general method of false position is expected to converge faster to the root than the Bisection Method. That should be also, as to generate next approximation, it takes account of function values by the formula.<\/p>\r\n \r\n

c\u00a0 \u00a0<\/em>= <\/span>a <\/em>\u00d7f <\/em>(b<\/em>) -b <\/em>\u00d7f <\/em>(a<\/em>)<\/span><\/p>\r\n

k\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0\u00a0<\/em>f <\/em>(<\/span>b<\/em>) - <\/span>f <\/em>(<\/span>a<\/em>)<\/span><\/p>\r\n\r\n<\/div>\r\n

\r\n

and c<\/em>k\u00a0<\/em>is just not the mid point of the interval.\u00a0 c<\/em>k\u00a0<\/em>is the point, where straight line\u00a0joining the two points (a, f(a)) and (b, f(b)) on the curve intersects with X axis. One very important measure of speed of convergence is \u201cOrder of Convergence\u201d. We shall discuss over \u201cOrder of Convergence\u201d in one of next modules. Method of False Position is regarded as having superlinear convergence.<\/span><\/p>\r\n \r\n

It can be seen, from the Figures [2,3,4,5,6 ], as well as computational tables above, one endpoint has become fixed. The moment curve does not change its shape from concave to convex or vice versa in the bracketing interval [a, b], always in this method, <\/span>c<\/em>k <\/em>\u2019s will approach the root from one side. As a result <\/span>c<\/em>k <\/em>\u2019s successively replace one of the two <\/span>a<\/em>\u2019s or <\/span>b<\/em>\u2019s, as the case may be, and the other remains fixed. For example, in example 1 as well as example 2, from the very beginning <\/span>b <\/em>is remaining constant. Other shape figures are also presented below for completeness.<\/span><\/p>\r\n\r\n<\/div>\r\n

Case\u20133:\u00a0<\/span><\/strong><\/p>\r\n\r\n

\r\n\r\n\"\"\r\n\r\n<\/div>\r\n
\r\n\r\n[a <\/em>remains fixed throughoutthe process. Only b <\/em>is changing, sequence ofc'sthat is Its approach the root)\r\n\r\n<\/div>\r\nCase\u20134:\u00a0<\/span><\/strong>\r\n
\r\n\r\n\"\"\r\n\r\n<\/div>\r\n
\r\n\r\n \r\n

[ b <\/em>remains fixed, a <\/em>keeps on changing, sequence of c<\/em>k <\/em>\u2019s <\/strong>(a<\/em>\u2019s) approach the root ]<\/p>\r\n \r\n

Thus, as\u00a0<\/span>c<\/em>k <\/em>\u2019s approach the root only from one side of the root, interval length no more\u00a0<\/span>approaches zero, it will always be greater than either |<\/span>a <\/em>\u2013 root| or |<\/span>b <\/em>\u2013 root| depending upon which of the endpoint is fixed.<\/span><\/p>\r\n

e.g. If b point isfixed, bracketing subinterval length shall always be >I<\/span>b- <\/em>root 1.<\/span><\/p>\r\n\r\n<\/div>\r\n

\r\n\r\n\"\"\r\n

Thus, no more number of intervals can be predetermined to achieve desired accuracy.2) There could be cases, though rare, where Method False Position may get slower as compared to Bisection Method. For example, look at the following graph of the equation f(x) = x10 \u2013 1 = 0.<\/p>\r\n\r\n<\/div>\r\n

\r\n\r\n\"\"\r\n\r\n \r\n

[ False position method is slower than Bisection method ]<\/p>\r\n \r\n

In \u00a0this \u00a0example, \u00a0sequence \u00a0of \u00a0approximations \u00a0is \u00a0converging \u00a0slowly. \u00a0<\/span>c<\/em>k<\/em>\u2013 \u00a0<\/strong>c<\/em>k<\/em>-1\u00a0\u00a0\u00a0 is\u00a0<\/span>becoming small, though not sufficiently near the root. The reason being, slope of the curve is large, at and in neighborhood of the root. That is why; many solvers apply F \u2013 TOL in conjunction with X\u2013TOL stopping criterion to find roots of an equation using method of False Position Sequence Following solution tables give comparison of Bisection Method and method of False Position in this specific example, which confirms that for this example, Bisection Method is better.<\/span><\/p>\r\n

Table:<\/span><\/p>\r\n

f <\/em>(<\/span>x<\/em>) = <\/span>x<\/em>10 -1<\/span><\/p>\r\n

a <\/em>= 0, <\/span>f <\/em>(<\/span>a<\/em>) = 1.00000<\/span><\/p>\r\n

b <\/em>= 1.3, <\/span>f <\/em>(<\/span>b<\/em>) = 12.785842<\/span><\/p>\r\n\r\n<\/div>\r\n

\r\n

Comparison Table\u00a0<\/strong><\/p>\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
k<\/strong><\/td>\r\nBisection Method<\/strong><\/td>\r\nFalse Position Method<\/strong><\/td>\r\n<\/td>\r\n<\/td>\r\n<\/tr>\r\n
<\/td>\r\nc<\/strong>k<\/strong><\/td>\r\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\r\nc<\/strong>k<\/strong><\/td>\r\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\r\n<\/tr>\r\n
1<\/td>\r\n0.650000<\/td>\r\n-0.986537<\/td>\r\n0.094300<\/td>\r\n-1.000000<\/td>\r\n<\/tr>\r\n
2<\/td>\r\n0.975000<\/td>\r\n-0.223671<\/td>\r\n0.181759<\/td>\r\n-1.000000<\/td>\r\n<\/tr>\r\n
3<\/td>\r\n1.137500<\/td>\r\n2.626718<\/td>\r\n0.262874<\/td>\r\n-0.999998<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
18<\/td>\r\n0.999999<\/td>\r\n-0.000013<\/td>\r\n0.928885<\/td>\r\n-0.521790<\/td>\r\n<\/tr>\r\n
19<\/td>\r\n1.000001<\/td>\r\n0.000012<\/td>\r\n0.943436<\/td>\r\n-0.441369<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n.<\/td>\r\n<\/tr>\r\n
59<\/td>\r\n---<\/td>\r\n---<\/td>\r\n0.999999<\/td>\r\n-0.000014<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/div>\r\n \r\n\r\nAdvantages:<\/strong>\r\n\r\n \r\n\r\n1) It is faster than Bisection Method.\r\n\r\n2) It is also simpl\r\n\r\n3) It guarantees convergenc\r\n\r\n4) Only one function evaluation per iteration is required.\r\n\r\n \r\n

Disadvantages:<\/strong><\/p>\r\n \r\n\r\n1) Not self starting. One needs two initial guesses a and b such that f(a)\u2219<\/strong>f(b) < 0.\r\n\r\n2) Though, faster than Bisection Method, still regarded as slow.\r\n\r\n3) In rare cases, it may become slower than Bisection Method.\r\n\r\n \r\n

-----------------o-------------o----------------o------------------o-----------------o------------o<\/p>\r\n\r\n\r\n\r\n\r\n
you can view video on Iterative Methods to solve equation f (x) = 0: Method of False Position<\/strong><\/td>\r\n\"\"<\/a><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
Suggested Reading:<\/strong><\/div>\r\n
<\/div>\r\n
1. Numerical Methods for Engineers by Steven C Chapra & Raymond P Canale, Fifth Edition, Tata McGraw Hill Publication, Special Indian Edition.<\/div>\r\n
2.\u00a0A Friendly Introduction to Numerical Analysis by Brian Bradie, Pearson Education.<\/div>\r\n
3.\u00a0Numerical Mathematics and Computing by Ward Cheney & David Kincaid, fifth Edition, Cengage Learning.<\/div>\r\n
4.\u00a0Computer Oriented Numerical Methods by Dr. N Datta, Vikas Publication.<\/div>\r\n
5. Numerical Methods with Programs in C by T Veerarajan & T Ramachandran, Second Edition, Tata McGraw<\/div>\r\n
Hill Publication.<\/div>\r\n
6.\u00a0Numerical Methods by V. Rajaraman, Third Edition, Prentice - Hall India Pvt. Ltd.<\/div>\r\n
7.\u00a0Numerical Methods with C++ Programming by RM Somasundaram & RM Chandrasekaran, Prentice -Hall India Pvt. Ltd.<\/div>\r\n
8. Applied Numerical Analysis by C F Gerald & P O Wheatley, Seventh Edition,\u00a0Pearson Education Asia, New Delhi.<\/div>\r\n
9. Numerical Methods by Dr. V. N. Vedamurthy & Dr. N.Ch. S.N. Iyengar, Vikas Publication.<\/div>\r\n
10. Numerical Analysis by Richard L. Burden, J. Douglas Faires, Cengage Publishcation.<\/div>\r\n
11. Numerical Methods with programs in BASIC, FORTRAN, Pascal and C++ by S. Balachandra Rao & C. K. Shantha, Revised Edition, Universities Press.<\/div>\r\n
12.\u00a0A Textbook of Computer Based Numerical and Statistical Techniques by A. K. Jaiswal & Anju Khandelwal, New Age International (P) Ltd, Publishers.<\/div>\r\n
13.\u00a0http:\/\/www.math.niu.edu\/~dattab\/MATH435.2013\/<\/div>\r\n
14.https:\/\/ece.uwaterloo.ca\/~dwharder\/NumericalAnalysis\/10RootFinding\/falseposition\/<\/div>\r\n ","rendered":"
\"epgp<\/a>
\n<\/span><\/div>\n
\n

 <\/p>\n

3.1\u00a0 Introduction\u00a0<\/strong><\/p>\n

 <\/p>\n

In \u00a0our \u00a0previous \u00a0lecture \u00a0(module), \u00a0we \u00a0learned \u00a0about \u00a0nature \u00a0of \u00a0iterative \u00a0methods, classification \u00a0of \u00a0iterative \u00a0methods \u00a0(bracketing \u00a0or \u00a0open) \u00a0to \u00a0find \u00a0roots \u00a0of \u00a0an \u00a0equation\u00a0f <\/em>(<\/span>x<\/em>) = 0<\/span>. We began with Bisection method to find roots of <\/span>f <\/em>(<\/span>x<\/em>) = 0 . It was observed\u00a0<\/span>that through Bisection method is a very simple, steady and reliable method, it suffers from drawback of being slow. It would be advisable to seek its modification to make it faster. The regular falsi method or method of false position is an attempt in that direction. This also is closed method. As in, Bisection method, at the beginning of iterative process\u00a0<\/span>two end points <\/span>a <\/em>and <\/span>b <\/em>are found so that\u00a0f <\/em>(a<\/em>)\u00b7 f <\/em>(b<\/em>) < 0 to ensure that at least one root lies\u00a0<\/span>between <\/span>a <\/em>and <\/span>b<\/em>.<\/span><\/p>\n

 <\/p>\n

It was also pointed out in disadvantages of Bisection method that it does not take into account of magnitude of function values at end points to generate estimates for the next iteration. It just considers sign of function values. This makes it slow.<\/span><\/p>\n

 <\/p>\n

Thus, in method of false position, the function values at the end points <\/span>a <\/em>and <\/span>b <\/em>are taken into account to generate the estimates for next iteration.<\/span><\/p>\n

 <\/p>\n

Method of false position differs with Bisection method only in the \u00a0methodology \u00a0of finding <\/span>c <\/em>. In Bisection Method, interval was subdivided into equal parts, \u00a0<\/span>c <\/em>was mid point given by (<\/span>a <\/em>+ <\/span>b<\/em>) 2.<\/span><\/p>\n

 <\/p>\n

No account was taken of function values at <\/span>a <\/em>and <\/span>b <\/em>. As at the root, function value is zero, the point, where function value is smaller is expected to be nearer to the root. It can be visualized also from the graph below.<\/span><\/p>\n<\/div>\n

\n

\"\"<\/p>\n

 <\/p>\n

3.2 Iterative Process\u00a0<\/strong><\/p>\n

 <\/p>\n

So, in regular falsi method, to take into consideration the function values at \u00a0a <\/em>and b <\/em>,\u00a0straight line is drawn, joining\u00a0<\/span>(<\/span>a<\/em>, <\/span>f <\/em>(<\/span>a<\/em>)) and (<\/span>b<\/em>, <\/span>f <\/em>(<\/span>b<\/em>)) .<\/span><\/p>\n<\/div>\n

\n

 <\/p>\n

The point, where it cuts \u00a0X \u00a0<\/em>axis is, the new estimate of the root. Mathematically, estimate of formula for c <\/em>can be derived as follows:<\/p>\n<\/div>\n

Equation of straight line joining\u00a0<\/span>(<\/span>a<\/em>, <\/span>f <\/em>(<\/span>a<\/em>)) and (<\/span>b<\/em>, <\/span>f <\/em>(<\/span>b<\/em>))\u00a0<\/span>is given by :<\/span><\/p>\n

\n

\"\"<\/p>\n

If \u00a0f <\/em>(c<\/em>) = 0,c <\/em>is the root otherwise if f <\/em>(a<\/em>)\u00b7 f <\/em>(c<\/em>) > 0 , \u00a0c <\/em>takes the role of \u00a0a <\/em>otherwise \u00a0c \u00a0<\/em>takes the role of b<\/em>.<\/p>\n

 <\/p>\n

3.3\u00a0 Graphical illustration of the Iterative Process\u00a0<\/strong><\/p>\n

 <\/p>\n

Case \u2013 1:f(<\/em><\/strong>x)<\/em><\/strong><\/p>\n<\/div>\n

\n

\"\"<\/p>\n<\/div>\n

\"\"<\/p>\n

3.4\u00a0 Different Stopping Criterion\u00a0<\/strong><\/p>\n

 <\/p>\n

Now, once we have observed the iterative process, let us focus on different stopping criterions. The sequence of <\/span>c <\/em>k\u2019s approaches the root. That is why, in case of the Bisection section algorithm, iterative process was stopped, when<\/span><\/p>\n

 <\/p>\n

| <\/strong>c <\/em>k <\/strong>\u2013 <\/strong>c <\/em>k-1<\/strong>| <<\/strong>\u00ce <\/span>\u2026 (3)<\/strong><\/p>\n

\n

 <\/p>\n

It is called x \u2013 tolerance criterion (X-TOL<\/strong>). The recent most c <\/em>k is the estimate of the root. There is another stopping criterion called function tolerance (F<\/strong>\u2013TOL)<\/strong>. At root, function value is zero, so if estimate is quite near the root then function value would be small. Hence, many times, one would like to stop when<\/p>\n

 <\/p>\n

|f ( <\/strong>c <\/em>k<\/strong>)| < <\/strong>d \u00a0\u2026 (4)<\/strong><\/p>\n

 <\/p>\n

where d is some small positive constant like \u20ac.<\/p>\n

 <\/p>\n

You would like to ask: <\/span>which one is better or what basis decision should be made to choose the particular criterion<\/em>? Well, it depends on the method, nature of function near the root, objective of finding the root etc. Many times, a combination of different criterions are also applied.<\/span><\/p>\n

 <\/p>\n

If slope of the curve is small near the root, curve is almost horizontal, and then function tolerance may not be appropriate stopping criterion, because curve is rising slowly, function values in neighborhood of the root are going to be small, so even if estimate <\/span>c <\/em>is not sufficiently near the root, one may stop. This is visible from the following graph:<\/span><\/p>\n<\/div>\n

\n

\"\"<\/p>\n

On the other hand, if it is expected that slope is high near the root or sequence\u00a0 of approximations c <\/em>k\u2019s may converge to the root, like in case of almost vertical graph, then function tolerance ensures that estimate is good approximation to the root.<\/p>\n<\/div>\n

\n

\"\"<\/p>\n

In simple words, for <\/span>f <\/em>(<\/span>x<\/em>) , if\u00a0<\/span>f <\/em>\u00a2(<\/span>x<\/em>)\u00a0\u00a0<\/span>is high near the root, FTOL would be better to use as\u00a0<\/span>stopping criterion.\u00a0<\/span><\/p>\n<\/div>\n

\"\"<\/p>\n

\"\"<\/p>\n

3.5 Algorithm for Regular False Method\u00a0<\/span><\/strong><\/p>\n

 <\/p>\n

\"\"<\/div>\n

 <\/p>\n

Let us solve examples to find root of an equation using method of false position. For comparison purpose, we shall find the roots of the same equations, used for illustrating Bisection Method.<\/span><\/p>\n

\n

 <\/p>\n

Example 1:<\/p>\n

 <\/p>\n

f <\/em>(x<\/em>) = 3x <\/em>– cos x <\/em>-1<\/p>\n

 <\/p>\n

Take a <\/em>= 0, f <\/em>(a<\/em>) = -2.0, c<\/em>0 = a<\/em><\/p>\n

 <\/p>\n

b <\/em>= 1, f <\/em>(b<\/em>) = 1.459698<\/p>\n

f <\/em>(a<\/em>). f <\/em>(b<\/em>) < 0<\/p>\n

 <\/p>\n\n\n\n\n\n\n\n\n\n
k<\/strong><\/td>\na<\/strong><\/td>\nf(a)<\/strong><\/td>\nb<\/strong><\/td>\nf(b)<\/strong><\/td>\nc<\/strong>k<\/strong><\/td>\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\nc<\/strong>k<\/strong>– c<\/strong>k-1<\/strong><\/td>\n<\/tr>\n
1<\/td>\n0.000000<\/td>\n-2.000000<\/td>\n1.000000<\/td>\n1.459698<\/td>\n0.578085<\/td>\n-0.103255<\/td>\n0.578085<\/td>\n<\/tr>\n
2<\/td>\n0.578085<\/td>\n-0.103255<\/td>\n1.000000<\/td>\n1.459698<\/td>\n0.605959<\/td>\n-0.004081<\/td>\n0.027873<\/td>\n<\/tr>\n
3<\/td>\n0.605959<\/td>\n-0.004081<\/td>\n1.000000<\/td>\n1.459698<\/td>\n0.607057<\/td>\n-0.000159<\/td>\n0.001099<\/td>\n<\/tr>\n
4<\/td>\n0.607057<\/td>\n-0.000159<\/td>\n1.000000<\/td>\n1.459698<\/td>\n0.607100<\/td>\n-0.000006<\/td>\n0.000043<\/td>\n<\/tr>\n
5<\/td>\n0.607100<\/td>\n-0.000006<\/td>\n1.000000<\/td>\n1.459698<\/td>\n0.607102<\/td>\n-0.000000<\/td>\n0.000002<\/td>\n<\/tr>\n
6<\/td>\n0.607102<\/td>\n-0.000000<\/td>\n1.000000<\/td>\n1.459698<\/td>\n0.607102<\/td>\n-0.000000<\/td>\n0.000000<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n

Example 2:<\/p>\n

 <\/p>\n

f <\/em>(x<\/em>) = x<\/em>3 + x <\/em>-1<\/p>\n

 <\/p>\n

Take a <\/em>= 0, f <\/em>(a<\/em>) = -1,c0 = a <\/em>= 0<\/p>\n

 <\/p>\n

b <\/em>=1, f <\/em>(b<\/em>) =1<\/p>\n

 <\/p>\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
k<\/strong><\/td>\na<\/strong><\/td>\nf(a)<\/strong><\/td>\nb<\/strong><\/td>\nf(b)<\/strong><\/td>\nc<\/strong>k<\/strong><\/td>\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\nc<\/strong>k<\/strong>– c<\/strong>k-1<\/strong><\/td>\n<\/tr>\n
1<\/td>\n0.000000<\/td>\n-1.000000<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.500000<\/td>\n-0.375000<\/td>\n0.500000<\/td>\n<\/tr>\n
2<\/td>\n0.500000<\/td>\n-0.375000<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.636364<\/td>\n-0.105935<\/td>\n0.136364<\/td>\n<\/tr>\n
3<\/td>\n0.636364<\/td>\n-0.105935<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.671196<\/td>\n-0.026428<\/td>\n0.034832<\/td>\n<\/tr>\n
4<\/td>\n0.671196<\/td>\n-0.026428<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.679662<\/td>\n-0.006376<\/td>\n0.008466<\/td>\n<\/tr>\n
5<\/td>\n0.679662<\/td>\n-0.006376<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.681691<\/td>\n-0.001525<\/td>\n0.002029<\/td>\n<\/tr>\n
6<\/td>\n0.681691<\/td>\n-0.001525<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.682176<\/td>\n-0.000364<\/td>\n0.000485<\/td>\n<\/tr>\n
7<\/td>\n0.682176<\/td>\n-0.000364<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.682292<\/td>\n-0.000087<\/td>\n0.000116<\/td>\n<\/tr>\n
8<\/td>\n0.682292<\/td>\n-0.000087<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.682319<\/td>\n-0.000021<\/td>\n0.000028<\/td>\n<\/tr>\n
9<\/td>\n0.682319<\/td>\n-0.000021<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.682326<\/td>\n-0.000005<\/td>\n0.000007<\/td>\n<\/tr>\n
10<\/td>\n0.682326<\/td>\n-0.000005<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.682327<\/td>\n-0.000001<\/td>\n0.000002<\/td>\n<\/tr>\n
11<\/td>\n0.682327<\/td>\n-0.000001<\/td>\n1.000000<\/td>\n1.000000<\/td>\n0.682328<\/td>\n-0.000000<\/td>\n0.000000<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n

Programming Code in C:<\/p>\n

The above output was obtained by running the following C program.<\/p>\n

 <\/p>\n

\/\/Regular falsi method<\/p>\n

 <\/p>\n

#include<stdio.h><\/p>\n

#include<conio.h><\/p>\n

#include<math.h><\/p>\n

#define fnx(x) pow (X, 3)-X+1<\/p>\n

 <\/p>\n

void main()<\/p>\n

{<\/p>\n

float a, b, c,fa, fb, fc, error,co; int i=0;<\/p>\n

int count=0;<\/p>\n

float err = pow(10.0,-6.0)\/2;<\/p>\n

 <\/p>\n

printf(“\\nEnter initial value of a = “); scanf(“%f”, &a);<\/p>\n

fa = fnx(a);<\/p>\n

printf(“\\n value of f(a) = %f”, fa);<\/p>\n

 <\/p>\n

printf(“\\nEnter initial value of b = “); scanf(“%f”, &b);<\/p>\n

fb = fnx(b);<\/p>\n

printf(“\\n value of f(b) = %f”, fb); co=a;<\/p>\n

printf(“\\nNumber\\t A\\t F(A)\\t B\\t F(B)\\t C\\t F(C)\\t Error”);<\/p>\n

 <\/p>\n

do<\/p>\n<\/div>\n

\n

{<\/p>\n

i++;<\/p>\n

c = (a*fb – b*fa) \/ (fb – fa);<\/p>\n

fc = fnx(c); if(fc==0)<\/p>\n

{<\/p>\n

printf(“\\ndesired root is %f”, c);<\/p>\n

\/\/exit(0);<\/p>\n

}<\/p>\n

error=fabs(co-c); co=c;<\/p>\n

printf(“\\n\\n%d\\t %f\\t %f\\t %f\\t %f\\t %f\\t %f\\t %f”, i, a, fa, b, fb, c, fc, error); if(fa*fc >0)<\/p>\n

{<\/p>\n<\/div>\n

}<\/span><\/p>\n

\n

else<\/p>\n

{<\/p>\n

 <\/p>\n

}<\/p>\n

a = c; fa= fc;<\/p>\n

 <\/p>\n

b = c; fb = fc;<\/p>\n<\/div>\n

\n

}while(error >err);<\/p>\n

 <\/p>\n

printf(“\\n\\nThe root of the equation = %f”, c);<\/p>\n

 <\/p>\n

getch();<\/p>\n

}<\/p>\n

 <\/p>\n

3.7 Observations\u00a0<\/strong><\/p>\n

 <\/p>\n

1)\u00a0 From the tables, it is clear that same desired accuracy is obtained in less number of iterations. Question is, is this the\u00a0 case for these two specific equations, or in general we can expect False Position Method to perform better as compared to Bisection Method. Well, in general method of false position is expected to converge faster to the root than the Bisection Method. That should be also, as to generate next approximation, it takes account of function values by the formula.<\/p>\n

 <\/p>\n

c\u00a0 \u00a0<\/em>= <\/span>a <\/em>\u00d7f <\/em>(b<\/em>) –b <\/em>\u00d7f <\/em>(a<\/em>)<\/span><\/p>\n

k\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0\u00a0<\/em>f <\/em>(<\/span>b<\/em>) – <\/span>f <\/em>(<\/span>a<\/em>)<\/span><\/p>\n<\/div>\n

\n

and c<\/em>k\u00a0<\/em>is just not the mid point of the interval.\u00a0 c<\/em>k\u00a0<\/em>is the point, where straight line\u00a0joining the two points (a, f(a)) and (b, f(b)) on the curve intersects with X axis. One very important measure of speed of convergence is \u201cOrder of Convergence\u201d. We shall discuss over \u201cOrder of Convergence\u201d in one of next modules. Method of False Position is regarded as having superlinear convergence.<\/span><\/p>\n

 <\/p>\n

It can be seen, from the Figures [2,3,4,5,6 ], as well as computational tables above, one endpoint has become fixed. The moment curve does not change its shape from concave to convex or vice versa in the bracketing interval [a, b], always in this method, <\/span>c<\/em>k <\/em>\u2019s will approach the root from one side. As a result <\/span>c<\/em>k <\/em>\u2019s successively replace one of the two <\/span>a<\/em>\u2019s or <\/span>b<\/em>\u2019s, as the case may be, and the other remains fixed. For example, in example 1 as well as example 2, from the very beginning <\/span>b <\/em>is remaining constant. Other shape figures are also presented below for completeness.<\/span><\/p>\n<\/div>\n

Case\u20133:\u00a0<\/span><\/strong><\/p>\n

\n

\"\"<\/p>\n<\/div>\n

\n

[a <\/em>remains fixed throughoutthe process. Only b <\/em>is changing, sequence ofc’sthat is Its approach the root)<\/p>\n<\/div>\n

Case\u20134:\u00a0<\/span><\/strong><\/p>\n

\n

\"\"<\/p>\n<\/div>\n

\n

 <\/p>\n

[ b <\/em>remains fixed, a <\/em>keeps on changing, sequence of c<\/em>k <\/em>\u2019s <\/strong>(a<\/em>\u2019s) approach the root ]<\/p>\n

 <\/p>\n

Thus, as\u00a0<\/span>c<\/em>k <\/em>\u2019s approach the root only from one side of the root, interval length no more\u00a0<\/span>approaches zero, it will always be greater than either |<\/span>a <\/em>\u2013 root| or |<\/span>b <\/em>\u2013 root| depending upon which of the endpoint is fixed.<\/span><\/p>\n

e.g. If b point isfixed, bracketing subinterval length shall always be >I<\/span>b- <\/em>root 1.<\/span><\/p>\n<\/div>\n

\n

\"\"<\/p>\n

Thus, no more number of intervals can be predetermined to achieve desired accuracy.2) There could be cases, though rare, where Method False Position may get slower as compared to Bisection Method. For example, look at the following graph of the equation f(x) = x10 \u2013 1 = 0.<\/p>\n<\/div>\n

\n

\"\"<\/p>\n

 <\/p>\n

[ False position method is slower than Bisection method ]<\/p>\n

 <\/p>\n

In \u00a0this \u00a0example, \u00a0sequence \u00a0of \u00a0approximations \u00a0is \u00a0converging \u00a0slowly. \u00a0<\/span>c<\/em>k<\/em>\u2013 \u00a0<\/strong>c<\/em>k<\/em>-1\u00a0\u00a0\u00a0 is\u00a0<\/span>becoming small, though not sufficiently near the root. The reason being, slope of the curve is large, at and in neighborhood of the root. That is why; many solvers apply F \u2013 TOL in conjunction with X\u2013TOL stopping criterion to find roots of an equation using method of False Position Sequence Following solution tables give comparison of Bisection Method and method of False Position in this specific example, which confirms that for this example, Bisection Method is better.<\/span><\/p>\n

Table:<\/span><\/p>\n

f <\/em>(<\/span>x<\/em>) = <\/span>x<\/em>10 -1<\/span><\/p>\n

a <\/em>= 0, <\/span>f <\/em>(<\/span>a<\/em>) = 1.00000<\/span><\/p>\n

b <\/em>= 1.3, <\/span>f <\/em>(<\/span>b<\/em>) = 12.785842<\/span><\/p>\n<\/div>\n

\n

Comparison Table\u00a0<\/strong><\/p>\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
k<\/strong><\/td>\nBisection Method<\/strong><\/td>\nFalse Position Method<\/strong><\/td>\n<\/td>\n<\/td>\n<\/tr>\n
<\/td>\nc<\/strong>k<\/strong><\/td>\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\nc<\/strong>k<\/strong><\/td>\nf(c<\/strong>k<\/strong>)<\/strong><\/td>\n<\/tr>\n
1<\/td>\n0.650000<\/td>\n-0.986537<\/td>\n0.094300<\/td>\n-1.000000<\/td>\n<\/tr>\n
2<\/td>\n0.975000<\/td>\n-0.223671<\/td>\n0.181759<\/td>\n-1.000000<\/td>\n<\/tr>\n
3<\/td>\n1.137500<\/td>\n2.626718<\/td>\n0.262874<\/td>\n-0.999998<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
18<\/td>\n0.999999<\/td>\n-0.000013<\/td>\n0.928885<\/td>\n-0.521790<\/td>\n<\/tr>\n
19<\/td>\n1.000001<\/td>\n0.000012<\/td>\n0.943436<\/td>\n-0.441369<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n.<\/td>\n<\/tr>\n
59<\/td>\n—<\/td>\n—<\/td>\n0.999999<\/td>\n-0.000014<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n

 <\/p>\n

Advantages:<\/strong><\/p>\n

 <\/p>\n

1) It is faster than Bisection Method.<\/p>\n

2) It is also simpl<\/p>\n

3) It guarantees convergenc<\/p>\n

4) Only one function evaluation per iteration is required.<\/p>\n

 <\/p>\n

Disadvantages:<\/strong><\/p>\n

 <\/p>\n

1) Not self starting. One needs two initial guesses a and b such that f(a)\u2219<\/strong>f(b) < 0.<\/p>\n

2) Though, faster than Bisection Method, still regarded as slow.<\/p>\n

3) In rare cases, it may become slower than Bisection Method.<\/p>\n

 <\/p>\n

—————–o————-o—————-o——————o—————–o————o<\/p>\n\n\n\n
you can view video on Iterative Methods to solve equation f (x) = 0: Method of False Position<\/strong><\/td>\n\"\"<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n
Suggested Reading:<\/strong><\/div>\n
<\/div>\n
1. Numerical Methods for Engineers by Steven C Chapra & Raymond P Canale, Fifth Edition, Tata McGraw Hill Publication, Special Indian Edition.<\/div>\n
2.\u00a0A Friendly Introduction to Numerical Analysis by Brian Bradie, Pearson Education.<\/div>\n
3.\u00a0Numerical Mathematics and Computing by Ward Cheney & David Kincaid, fifth Edition, Cengage Learning.<\/div>\n
4.\u00a0Computer Oriented Numerical Methods by Dr. N Datta, Vikas Publication.<\/div>\n
5. Numerical Methods with Programs in C by T Veerarajan & T Ramachandran, Second Edition, Tata McGraw<\/div>\n
Hill Publication.<\/div>\n
6.\u00a0Numerical Methods by V. Rajaraman, Third Edition, Prentice – Hall India Pvt. Ltd.<\/div>\n
7.\u00a0Numerical Methods with C++ Programming by RM Somasundaram & RM Chandrasekaran, Prentice -Hall India Pvt. Ltd.<\/div>\n
8. Applied Numerical Analysis by C F Gerald & P O Wheatley, Seventh Edition,\u00a0Pearson Education Asia, New Delhi.<\/div>\n
9. Numerical Methods by Dr. V. N. Vedamurthy & Dr. N.Ch. S.N. Iyengar, Vikas Publication.<\/div>\n
10. Numerical Analysis by Richard L. Burden, J. Douglas Faires, Cengage Publishcation.<\/div>\n
11. Numerical Methods with programs in BASIC, FORTRAN, Pascal and C++ by S. Balachandra Rao & C. K. Shantha, Revised Edition, Universities Press.<\/div>\n
12.\u00a0A Textbook of Computer Based Numerical and Statistical Techniques by A. K. Jaiswal & Anju Khandelwal, New Age International (P) Ltd, Publishers.<\/div>\n
13.\u00a0http:\/\/www.math.niu.edu\/~dattab\/MATH435.2013\/<\/div>\n
14.https:\/\/ece.uwaterloo.ca\/~dwharder\/NumericalAnalysis\/10RootFinding\/falseposition\/<\/div>\n

 <\/p>\n","protected":false},"author":4,"menu_order":4,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":["prof-savita-r-gandhi"],"pb_section_license":""},"chapter-type":[],"contributor":[58],"license":[],"class_list":["post-76","chapter","type-chapter","status-publish","hentry","contributor-prof-savita-r-gandhi"],"part":3,"_links":{"self":[{"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/chapters\/76","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/wp\/v2\/users\/4"}],"version-history":[{"count":7,"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/chapters\/76\/revisions"}],"predecessor-version":[{"id":584,"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/chapters\/76\/revisions\/584"}],"part":[{"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/parts\/3"}],"metadata":[{"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/chapters\/76\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/wp\/v2\/media?parent=76"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/pressbooks\/v2\/chapter-type?post=76"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/wp\/v2\/contributor?post=76"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/ebooks.inflibnet.ac.in\/itp15\/wp-json\/wp\/v2\/license?post=76"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}