Controlling Error and Conditional Statements:
Measuring error and the Residual:
If we are attempting to find a numerical solution of an equation f(x) = 0 then there are a few different ways we can measure the error of our approximation. The most direct way to calculate the error would be as:
{Error at step n} = en= xn− x
Where xn is the n-th approximation alsoxis the true value. Nevertheless we usually do not know thevalue of x, or we wouldn’t be trying to approximate it. This makes it not possible to know the errordirectly as well as so we must be more clever.
One possible strategy that habitually works is to run a program until the approximation xnstopschanging. The trouble with this is that sometimes doesn’t work. Just for the reason that the program stopchanging does not necessarily mean that xn is close to the true solution.
For Newton’s method we have the subsequent principle- At every step the number of significantdigits roughly doubles. While this is asignificant statement about the error (ever since it signifiesNewton’s method converges really quickly) it is to some extent hard to use in a program.
Rather than calculate how close xn is to x* in this and many other situations it is much more practicalto measure how close the equation is to being satisfied in other words how close yn= f(xn) is to0. We will use the quantity rn= f(xn) − 0, called the residual, in many different situations. Mostof the time we merely care about the size of rn, so we look at |rn| = |f(xn)|.
The if ... end statement
If we have a definite tolerance for |rn| = |f(xn)| then we are able to incorporate that into our Newtonmethod program using an if ... end statement-
function x = mynewton(f,f1,x0,n,tol)% Solves f(x) = 0 by doing n steps of Newton’s method starting at x0.% Inputs: f -- the function, input as an inline% f1 -- it’s derivative, input as an inline% x0 -- starting guess, a number% tol -- desired tolerance, prints a warning if |f(x)|>tol% Output: x -- the approximate solutionx = x0; % set x equal to the initial guess x0for i = 1:n % Do n timesx = x - f(x)/f1(x) % Newton’s formulaendr = abs(f(x))if r > tolwarning (’The desired correctness was not attained’)end
In this program if checks if abs(y) > tol is true or else else not. If it is true then it does the whole thingamong there and end. If not accurate then it omits ahead to end.
In the command window explains the function
> f = inline(’x^3-5’,’x’)
And its derivative
> f1 = inline(’3*x^2’,’x’).
Then utilize the program with n = 5, tol = .01, and x0 = 2. Next, change tol to 10−10 and repeat.The loop: while ... end
While the preceding program will tell us if it worked or not we still have to input n the number ofsteps to take. Therefore far for a well-behaved problem if we make n too small then the tolerance will notbe attained and we will have to go back as well as increase it or if we make n too big then the programwill take more steps than necessary.
One way to manage the quantity of steps taken is to iterate until the residual |rn| = |f(x)| = |y| issmall enough. In Mat lab this is merely accomplished with a while ... end loop.
function x = mynewtontol(f, f1, x0, tol)
x = x0; % set x equal to the initial guess x0y = f(x);while abs(y) > tol % Do until the tolerence is reached.x = x - y/f1(x) % Newton’s formulay = f(x)end
The statement while end is a loopalike to for end however instead of going through theloop a fixed number of times it keeps going so long as the statement abs(y) > tol is true.
One clear drawback of the program is that abs(y) might never get slighter than tol. If thishappens the program would carry onrunning over as well as over until we stop it. Try this by setting theforbearance to a really small number:
>tol = 10^(-100)
after that run the program again for f(x) = x3 − 5.
You are able to use Ctrl-c to stop the program when it’s stuck.
One way to evade an infinite loop is add a counter variable say i in addition to a maximum number ofiterations to the programs.
Utilizing the while statement this can be accomplished as:
function x = mynewtontol(f,f1,x0,tol)x = x0; i=0; % set x equal to the initial guess x0. set counter to zeroy = f(x);while abs(y) > tol & i < 1000 % Do until the tolerence is reached or max iter.x = x - y/f1(x) % Newton’s formulay = f(x)i = i+1;end
Latest technology based Matlab Programming Online Tutoring Assistance
Tutors, at the www.tutorsglobe.com, take pledge to provide full satisfaction and assurance in Matlab Programming help via online tutoring. Students are getting 100% satisfaction by online tutors across the globe. Here you can get homework help for Matlab Programming, project ideas and tutorials. We provide email based Matlab Programming help. You can join us to ask queries 24x7 with live, experienced and qualified online tutors specialized in Matlab Programming. Through Online Tutoring, you would be able to complete your homework or assignments at your home. Tutors at the TutorsGlobe are committed to provide the best quality online tutoring assistance for Matlab Programming Homework help and assignment help services. They use their experience, as they have solved thousands of the Matlab Programming assignments, which may help you to solve your complex issues of Matlab Programming. TutorsGlobe assure for the best quality compliance to your homework. Compromise with quality is not in our dictionary. If we feel that we are not able to provide the homework help as per the deadline or given instruction by the student, we refund the money of the student without any delay.
Theory and lecture notes of The QR Method all along with the key concepts of the qr method, linear algebra, assignment help, homework help. Tutorsglobe offers homework help, assignment help and tutor’s assistance on The QR Method.
Vibrational Spectroscopy tutorial all along with the key concepts of Vibrational coordinates, Newtonian mechanics, Quantum mechanics, Intensities
tutorsglobe.com chemical co-ordination assignment help-homework help by online co-ordination systems tutors
Limitation of Refraction Methods tutorial all along with the key concepts of Direct Waves, Vertical Velocities, Hidden Layers, Blind Zones, best-fit arrival, Cross-over distance
Photometry tutorial all along with the key concepts of Photometry and eye, Photometric quantities, Photometric versus radiometric quantities, Photometric measurement methods
www.tutorsglobe.com offers oxidation homework help, oxidation assignment help, online tutoring assistance, organic chemistry solutions by online qualified tutor's help.
tutorsglobe.com structure of sarcomere assignment help-homework help by online muscles tutors
tutorsglobe.com clinical manifestations of lyme disease assignment help-homework help by online lyme disease tutors
the neutral atmosphere tutorial all along with the key concepts of equation of state, equation of hydrostatic equilibrium, snell's law, water vapor, relative humidity, mixing ratio, partial pressure of saturated air, propagation delay, refractivity
Get top-rated General Ecology Assignment Help service from PhD experts and get authentic paper at viable prices to secure A++
Winding wire is solid wire, that to permit closer winding while making electromagnetic coils, is insulated just with varnish, than the thicker plastic or another insulation generally employed on electrical wire.
tutorsglobe.com physiological effects of sleep assignment help-homework help by online sleep tutors
Fate of Organic Matter in Sedimentary Basins tutorial all along with the key concepts of Accumulation of Organic Matter, Diagenesis, Catagenesis, Metagenesis, Transformation of Organic Matter and Kerogen to Petroleum
Dyeing Mechanisms tutorial all along with the key concepts of Textiles, Nontextile materials, Process of Dyeing, natural dye, Cellulose fibers
Respiration tutorial all along with the key concepts of Retrieving Glucose from other Molecules, Retrieval from Sucrose, Glycolysis, Potential Energy of Glucose, Breakdown of Glucose to Pyruvic Acid and The Krebs cycle
1963776
Questions Asked
3689
Tutors
1463723
Questions Answered
Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!