© Distribution of this video is restricted by its owner
Transcript ×
Auto highlight
Font-size
00:03 Yeah. Something people want to. , so continue that um inspirational

00:20 And I started at they're almost simplest . There's a simpler version of it's

00:30 . Mm hmm. So we'll talk little bit more about that then how

00:36 can do so chapters or go to something more accurate but using sensual

00:43 Sunday extrapolation. We're trapped. And I talked about another matter. Eventually

00:51 is same sort of there's always something . Alright. So here about Yes

01:00 remind her. But the tropical rule was 6 3 approximate the integral or

01:07 area under the curve area under or by chapter swords. Or there is

01:17 the straight plane piece by speed. approximation of the functions affect a bunch

01:24 notes. And you evaluate the function those notes and then rule as it

01:31 things correctly. There's the function between two or two neighboring points behave.

01:39 a straight cut and I'll come back that thing. So Sandy's going from

01:47 formulas were fairly simple in the end the best behalf of the function bodies

01:54 the and then just some of the . So that was 16 and a

02:07 version of it. And the reason that's the person with respect to the

02:14 of the interval. And the reason talking about it is I think we'll

02:19 it later essentially. It's an extrapolation is kind of based on mr cursed

02:28 of the tractors are so so in case the gun like when we talk

02:39 division they sort of refined somehow. or make it smaller and smaller as

02:44 function in this case and and has divided very attractive too forever stepped in

02:51 the final process. So H and former is in fact depending on them

02:58 on the given subject about being his and his father something. So this

03:08 that half of the 10 point values then in the summer of internal values

03:14 by the winds. Okay. Um as always. Well, that's that's

03:26 , yep. Mhm. So what now in this kind of notion of

03:32 intervals, manufacturers through and looking at finance, finance of the intervals?

03:40 you can see basically uh some of points, what they're saying us in

03:48 immediate course of of division. So only thing that happens is that in

03:55 you add some new functional evaluations and we change the weights of the functions

04:01 are actually and then we'll do Then we'll get this occurs the

04:08 So the first one was easy, the two endpoints and the width of

04:14 interval from the start finish from A B And just one track. Thanks

04:26 . 5 10 minutes. Some divided this can use the big point,

04:35 the intervals of april life simple. then they got the endpoint and the

04:42 points and that's the first step is The life does not demand as a

04:48 by two. And then there was of the influence with the board and

04:51 there's the other guy, the other once a month and rearrange the

04:59 So wow this thing I'll have the endpoints divided by two manipulating this survey

05:18 . This is half of this kind the new volume, that was the

05:26 point. That was interesting. Thank . So because again age is kind

05:35 one way of thinking about it. you split the interval between A and

05:41 . Into two intervals instead of Um The length just um as

05:46 so that's kind of the scaling factor that half a song. So and

05:56 there is the value for them so and then ah okay so do this

06:09 . It's like this formula and red the previous slide kind of related a

06:16 bit so manipulated from an exercise and what and also mostly invisible come

06:23 It's not worry too much about the endpoints porch and it's called it see

06:27 now and then we take a look what happens when some something in the

06:37 That is not true to the N to have Someplace to go from Say

06:41 from to this one. So too that age is kind of different normally

06:52 you use advice one because that would the interval divided by tweet today and

06:58 instead of the end. So to the same meaning of age and both

07:04 , they have to plug in the two here because Divided by 2 to

07:08 N -1 instead of two. So can so this is basically says that

07:14 and this former starts for the boats this form. So then we got

07:19 here that's also that replaced age but reach whether it's just a minute here

07:27 that she wasn't provided and I want play with that. That was from

07:33 previous slide 1st 1. So now you look at this difference that was

07:40 , is this one that comes from deal our end. Come on

07:47 And this is the one that comes this expression but we multiply this by

07:55 they have so they can see me see so that disappears an obvious to

07:59 about ah those three expressions now. since now it's just the same

08:12 right? So we can notice for the even um ice and this formula

08:19 corresponds best with 2, 1 Of function values and this one. So

08:27 it's only the odd ones. That's . So they aren't ones essentially the

08:35 But I knew that was not in state to to the to but only

08:39 the two For the three invigoration. space beside as well and this proportion

08:45 you're looking at it. So if have already evaluated this one you get

08:50 next one but there's a different function these options for the new points,

08:55 of the proper, so that's So, so this is based on

09:06 personal formula and wow. So summarizing , ah and the reason for this

09:23 we're going to use this formula You don't have to do it this

09:27 . Obviously you can start from scratch whatever function values Savant or on your

09:32 and the proceeds directly with us, you can also go through this the

09:41 of what you're doing stan The previous with a factor of two for a

09:49 divine gratitude. And then you have you function in it. Okay,

09:58 little bit of discourse before coming to book. So doing it for multi

10:10 is sticking straight forward and for that , but it's the first Look at

10:16 standard interval from 0 to 1 and it's just the scale and arbitrary intervals

10:23 standard interval And sometimes this community 0-1 sometimes it's convenient to use the symmetric

10:31 around zero. So from one is , the look of the scaling a

10:36 shots accounts later. So, so especially you can sort of recognize the

10:47 of rule because yeah, these are two endpoints, There are scale events

10:52 1/2 and then you have all the points on the two and the two

10:58 also this basically the sum of all internal points, everything may be

11:04 but it's expected and then make it simpler notation, just look at

11:11 Then the coefficients here, various It's basically the same for all the

11:15 points. It's best for them one here And then we got the

11:25 And so in order to do that thing, it can be real taking

11:32 dimension at the time which is done this thing. And the chapters are

11:36 with respect to X. And that's . In principle, I do it

11:42 all the wise. It's not done . Is to keep the integral and

11:48 two guys commute so you can sign and now the best to have waited

11:56 more while his wife is determined for chapters overall. In this case for

12:02 extraction. And then you can still the integral in the Y direction

12:09 you know something his own rules also the wine. And they get a

12:15 way to some respect to the white . So basically you get in this

12:22 intervention, you got some kind of great yeah xy notes. And then

12:31 function values those notes. And then have the weighting factors with respect to

12:40 um X integration. So that's all internal points added together. And then

12:46 have to be a little bit careful the boundaries will lead us to ask

12:50 for the border, depending upon where points are, but it's diseases are

13:00 the end, points in exile. the after boundaries for X. It's

13:05 heart and then in between it's wow, and then the same thing

13:11 the Y direction in the corner. what happens? Did you use that

13:18 of formula? But on something that satisfy the, well 12 seconds.

13:27 let's say integrating 1-1, one over or something. Okay. Yes.

13:37 if um so in that case it's signal point that zero. Um So

13:53 that's a good question. So I thank you. But that's a

14:08 question. Whether so, what's the of the rest? Um So that's

14:18 low content, right? Is the to everything or one overhead?

14:24 So the world function is not, simple zero. So then in that

14:31 it into $20,000 makes sense. So it, hmm. So if

14:42 remember figure out the being dependent on relative some water dependent in that

14:58 I guess the first term that was no one was second or third.

15:05 . But that's so what's the best relative of one over next against,

15:11 know, one over X. Squared stepped up with the science. So

15:15 affects the suspicion is more with the of stuff because for some reason,

15:25 that's why there's no bound up on . So correctness for several magnitudes.

15:33 I think, yeah. Okay, . So then it's just increasing one

15:47 two dimensions just by looking at the thing and so you invested in the

15:52 product from the approximation of respect. . Um So This Sunday's two points

16:09 , the Gander there were is proportional the number of points that have to

16:16 the function evaluation of so much and the most functions that people use or

16:22 end up using various contexts they work evaluation or computing function values. It's

16:30 lot more than they had multiplying. need to be careful. So,

16:40 now it has the same spacing in the different dimensions. So that's new

16:45 in the a number of points in dimension. The means that the work

16:56 the power of the number of So, and you know, functional

17:04 in one point and then you might of an N by N grid or

17:08 squared from evaluations they keep falling in to the number of dimensions. But

17:17 type of work. Whereas the accuracy determined by that edge is facing in

17:27 dimension so that it did not So that's what I'm trying to

17:31 So it's kind of a complicated This is straightforward. Um that so

17:37 this is basically, did you stab in as we will end to

17:43 Okay, powers. So this is one over n squared. That's been

17:50 to the square just to remember. that is in terms of the trip

18:01 over all. A broad separation in of what type of functions 15 cigarettes

18:08 well. It does not account for past that's doing more work. I

18:15 have more dimensions given age, grand and other things may actually get worse

18:22 this. That in respect to an position that this is the case and

18:29 independent. Okay. Dustin stop. Oh right. Um, what about

18:48 multidimensional, I'm not coming back to This recursive formula um, chapters on

18:56 . 20 questions. All right. this is um, um, it's

19:04 first of formula looks like uh, this is the way it looked at

19:11 sports, the person of yes, the chapter is older than the first

19:20 . This was no reversion. I just breaking it down and starting get

19:25 this point where the length of the because side revival but for each

19:41 that's it. Mm hmm. Dr martin, it's not okay for

19:49 , sorry. So now again, what I wanted to say. So

19:55 so called bromberg algorithm depends on The first formulation of the trump resort

20:04 , it's just supplying the switches and idea that they used for derivatives before

20:12 started, when we talked about they have reformed their percentage differences of

20:20 derivative and for that the error and approximation was um mr terms a bunch

20:30 terms of the docs expansion. There powers of X square. So I

20:36 a student. second power is the is 6 and Richardson was down good

20:44 that doing, Looking at H&H of two and producing Those two expressions to

20:53 1st the X squared term and you something that is accurate to the edge

21:00 the order and then it's complete cold they wanted more, what are your

21:05 ? So it turns out that for trap assad role the or terror also

21:12 in even powers of mhm. So this is kind of a scheme,

21:21 will talk about construction scheme levels actually for this um which is an extrapolation

21:27 the chances differences. You see the kind of triangle schema, you can

21:34 this in your first parliament for this of when it comes to the tractor

21:39 rule is for eight squared and then go to the next column here that

21:44 seen stage four and 86 etcetera that stopped. And I think that silence

21:52 for its moral is postulating the same as reconstructed. Well, it's kind

21:59 so it's the same format to get to go back and look at re

22:04 extrapolations conduct for derivatives. Now you from his Romberg integration for the one

22:13 , but you can see it's the structure and the same weight and it

22:18 from again that this is kind of formula for the track the songs.

22:29 so that gives rise to exactly the um schema as well, come for

22:36 and getting better approximations for integration for given age by doing India combinations of

22:47 that united in the first column that the source of the form. So

22:58 best is the same structures who did . So I think that's going to

23:04 or less the same exercise. But there's no difference except just using there

23:11 simple instead of the, the symbol the think about democratic expression because the

23:18 , the term looks exactly the same it comes to the trap is

23:23 And the whole thing is just to successive borders for powers or age they

23:30 the more active formula and the coefficients on the power. So um,

23:38 not on the coefficients in front so that's why we have the same

23:46 of the algorithm eliminate it where term this area ends up being so

24:01 yeah, continue as the in general looks just like the terror of being

24:21 powers so yes, right. That's and that's why it's best if we

24:40 say sure and I'm not going through , but Somebody did want to get

24:47 of the Air Formula four, the are low impact just has even hours

24:54 age. Stop it. So the thing you should remember now this particular

25:05 but that it's even for ourselves ancient system that would suck anyone interested.

25:17 it was pretty much and so I yes. So this is just one

25:21 thing and kind of think they have example of work. So they also

25:28 out, so how do you know sometimes the internet, Okay. Um

25:39 the table so one can look at particular expression and that production wait um

25:50 this formula here um and stick it um let's form it up so no

26:02 can show that this ratio should they behavior, it's this power is the

26:12 behaves in a way subject to the , that depends on the derivatives,

26:21 correct order depending upon um before you to go and jump on improving the

26:27 . Um Nothing welcome to interval. in quotes. The absolute value of

26:40 derivatives in the term that you are that. So arena or digital comics

26:48 too. The chapters on room. , so this is kind of the

26:55 . So as long as this comes more or less close to forward to

27:01 pharaoh was the number of columns. how many times supply the Richardson extrapolation

27:08 the Oh no ! So one can through and look at the expression,

27:16 fine. And then you can actually that uh can for m equals

27:24 The first thing you notice the first of the circulation and plug it into

27:31 ratio. But I said it's a test for whether that's circulation gives a

27:38 outcome. So on the budget in and then it was to manipulate this

27:44 like that. That's something that That's the four plus in the

27:50 Mhm. So it kind of Ah This is proved by example which

27:58 not been approved at least at And then this is just a reminder

28:13 what they've already said, remind you for and supporting the slide and then

28:21 the boat it shows the code for the Romberg type algorithm or extensions and

28:30 for them, not this one and only comment. They have you in

28:37 texas. Yes, you can go here and I think you can call

28:41 by columns. But the point of bromberg algorithm or it's an extrapolation and

28:51 our orders and things to get good before. So obviously not that small

28:58 . So basically, definitely this is point of care about. So you

29:03 want to do more functions because when go across the columns here, if

29:08 don't do any more function evaluations, the same number of function evaluations.

29:13 only thing that happens yes, you uh it doesn't show explicitly but it's

29:20 hola what's happening. And then the thing that happens that well, my

29:25 from one column to the next ah weights of the different function values you

29:31 . It's changed and I'll come back think so that's why it's here.

29:40 £7 kind of that can do trial error, you know, take an

29:48 and see what happens and that increased age. But yes, you said

29:53 , you know that they're so It's , 2, 4, 6,

29:59 here in this column. Then they estimate what they should need and they

30:03 need to go further down. And more functions about addiction and work in

30:07 column, they do function evaluations for person um to double the number of

30:18 more or less every time they moved increase. So so that's part of

30:25 show. Here's where the function of are and this is just a combination

30:30 what's your priorities completely. Yeah. , so this structure is called

30:40 It goes over That's all right. the most common something people ever

30:57 Yes. So the concrete for um for the book the integration of

31:10 , which one 15 points the next of the special again by age or

31:21 smaller. Aah and then what the result is straightforward application of the chapters

31:34 rule. And then in this case doing Romberg ah algorithm all the way

31:43 the right so that speaks to refuse the highest order accuracy you can get

31:48 the number of functions that are your you'll make. So that's the shows

31:53 that you get very good approximation in case for mr Uh huh is to

32:04 for 16 right function evaluations seven. I cannot versus to get the same

32:17 accuracy. You need to get some somewhere that's playing around many more functional

32:26 for straightforward. The first college aged and 6 and this causes the

32:56 Um So breaking the flow I would . But we'll talk about the general

33:08 on extrapolation in the general context of internship. Talking about the integration and

33:14 up with something. So this is shows like definitely it's just an idea

33:26 just fine and just a little bit messy. Is it's such a nice

33:36 of two comes of age. Uh . Some orbiter for south delta and

33:44 of age. But the same idea try to where There are three some

33:50 in that sense. So I was smaller than their smaller and governments that

33:56 want to grow the accuracy. They to get rid of the first things

34:00 and then they start to work on second time after you might eliminated the

34:05 . So just on the concept of you can still do the exercise on

34:11 state cutting H half. It's just now things are a bit last

34:16 And then so much thanks to you this format to space to the alpha

34:22 instead of the coefficient A They divided two or 4 but By 2

34:29 And so that means you have to things to to the answer. Um

34:34 out the same proficient for the first of H and L. This abuses

34:43 the village. But if you've done now they can do that. Take

34:48 difference between these two guys and they a sulfa turned his fierce and then

34:57 Mhm. The first area term that the next 100 power over age and

35:04 you can keep doing that. The level of that, that's the

35:10 leaving term that it's ignored. Um smallest. So this is thanks

35:23 then. Um you this was the step start to work on. There

35:33 go. Okay, so so that's best of the general, the idea

35:43 the general extrapolation nothing named to the but even powers are in this case

35:52 will interval that pitch but it's a bit last year. Nice and simple

36:01 of But conceptually it's not in fact I'm not so sure. Yes.

36:11 could also be another conference for this form. Now that's the age.

36:18 it's one over some variable X. of course in the Lord to access

36:24 small market is one over age. I went over it so it's very

36:30 . So but now they can work the one over X. Term.

36:34 that all H. And then instead getting smaller it should get larger because

36:42 want this can be diminished, they converted to their ideas, it's in

36:47 case it's exposed towards infinity that this here. That's too this So now

36:58 of looking for convergence for some getting smaller than this convergence question,

37:04 mechanical knowledge and then the and stuff extend half your double X. And

37:12 that's the same idea. Just expressions then you got something new and eventually

37:18 we've got this is kind of the term now and get rid of

37:24 And then you can start to work the next reading from. Yeah,

37:29 approximation. Same idea about doing linear or whatever. So, and this

37:41 a complete example the numbers or that double X for each case here and

37:49 what the function like this and it 1 28. It was about something

37:56 So one step forward. Yes. you. Extrapolation idea. Then it

38:03 something ah well you're welcome. That's not some story that Okay,

38:25 think so. I don't know how Well, I don't have it on

38:28 spins. It's not what I project the chair. All right.

38:37 Um, we'll see. Try but when you try to Okay.

38:50 , it stays think but I know can manually can like right on the

39:00 for the Was one of the Okay, alright. Um, so

39:12 was about the general circulation idea. don't want to come back to integration

39:27 . Mm hmm, conceptually what they're in the comics is slightly different the

39:34 and attractive. So now like taking different you on the integration to try

39:49 ah mr driven by kind of a model or accuracy about the incident.

40:02 and this is who said that that's happened when they got or the approximation

40:12 the Stromberg, the idea of this and extrapolation. But the rule is

40:18 didn't even more function about the ever for different functions evaluations, they just

40:24 the white aggressively by doing the combination in the previous columns. So,

40:40 going after going to the role then if you just have one chapters on

40:49 there was just there for the function of the U. N. And

40:54 true for every we can manage our through for each one of these resources

41:00 them to take their average height. that's actually what it was. That

41:08 kind of a linear combination of the endpoints values. And for the shepherds

41:13 drool, A and B was the now. And you also know that

41:22 state took place the curve of the F for the straight line, the

41:28 of. So that was the tractor also what it says if to integrate

41:35 straight function, that is a straight between the two endpoints, A and

41:41 . R. A and B and . Then there is no error or

41:50 . So it's another way of looking it. So if in fact I

41:55 the function was this great times, is no approximation error and maybe the

42:05 error during the um occupies and all that but there's no approximation or translation

42:11 . They are starting to send annual finally and it takes years especially.

42:18 this is kind of trying to look integration formulas using this idea of what

42:25 of polynomial? Yes integrated, So for the trappers or rule the

42:34 is what's been your function from A B. It just takes the tractor

42:38 already calling for the function then you correct with no approximation. So so

42:49 look at that. So that means that's the first order polynomial. So

42:55 the purpose of road integrates any first and over correctly. Um First of

43:05 polynomial is a concept that So that should be integrated clean.

43:12 now reduced instead of the arbitrary into 8:30 for simplicity of notation. And

43:18 we'll figure out how to make a In a bit but there's the tsunami

43:23 stay in speed from 0 to 1 sustained they integrate the constant bond and

43:29 value is one. So that means the sum of it because function value

43:35 influences one was a function. So this place they have this notion of

43:40 combination of two. 10.5. This a plastic then this is a sure

43:48 thing going straight down the street. are the one this X. So

43:54 integrates X from 0 to 1 is squared over two. Now the function

43:59 of Um X0 is zero. And it's only the function value 10

44:05 So it's quite that one as So that's now. So now that

44:11 that the effect the function value in that formula um for ex if zero

44:18 zero so eight disappears. So we get mm hmm. So now we

44:24 two equations ah two for the two , A&B. And obviously you will

44:31 see that that day is also a . So what we got is

44:38 to do the centerville correctly. Uh error than in the half of what

44:43 got this formula which is exactly the is over. Uh but now basically

44:50 by not just taking off as a is on the street plan and they

44:57 say, well fallen order gets integrated . And we designed the format,

45:04 the first degree polynomial, correct. then I can just change it to

45:12 an arbitrary interval basically. So it like this formula. So for X0

45:18 why? It's a. And for the by the cuts being. So

45:23 doing this variable kind of substitution Yet this formula from 0 to 1

45:29 supplied to norbit orient to both A tv. And then they get

45:35 scaling factor because of the length of interval. Thank goodness for this

45:42 So now the rest of the general or rule. But now I'm starting

45:50 the correct integration of a polynomial, this case? Of course.

45:58 So survival. So now. so the next thing is well okay

46:12 slowly. So it's now I have by the way is not the

46:21 the simplest for the integrate constant. this for me um this they started

46:32 very simplistic. That means integrating step , constant function, staircase them

46:46 Straight lines. And the next thing wow, maybe I want something that

46:52 on a mobile first. So then same way the that time they pursued

47:01 and now basically I want to attend bad idea. Sizing take now three

47:08 um the two endpoints and then That's . We need three points to define

47:17 quadratic formula. That's very simple. Kind of the simplest thing is just

47:23 take a deep .15. There are ways of picking up boxes now that's

47:34 . Um mm hmm. So now idea is that fish for the infrequent

47:41 this function. But if the approximate in our academic then they want this

47:49 of the normal to be Exactly. this is kind of the professional woman

47:58 . And that means we have assume because when this We needed the three

48:14 values or three notes and correspondent function to define the polynomial and they just

48:22 to have um the function values at end points and then the reports and

48:30 of a narrative with the traveler's sword with again the evaluation function you can

48:41 and I think that's one of them the next year is nice but you

48:44 also I think it's the polynomial is using the lagrange functions. Then there

48:55 this cardinal polynomial and they used function the nodes used them for collections.

49:04 in the spirit of the ground formulation wow different in the polynomial because the

49:13 polynomial has depended on X. The it only depends on the chosen noble

49:20 but it's not that. So that the functions are just used in the

49:28 they worked to plug in Lagrange formulation have to get the constants then those

49:36 have the function values and A. B. And C. Will he

49:43 to the integration of the cardinal Oh Mhm. So no I

49:53 For simplicity. Now we have three . So it's kind of convenience and

49:57 a good point Correspondent. Except for and the Q. and points.

50:03 person nine is 1 for the It kind of is symmetric standard in

50:08 morning. So now and this idea it's location. Now we should have

50:16 The three function values there are such if f. Of X happens to

50:21 a quadratic polynomial then the integral All right. Same thing again.

50:31 we integrate the constant X. And squared. So these are the Up

50:36 the 2nd 3 Polynomial. And that from -1-1. Then The concept the

50:49 is intermodal length is two and the is one of the two. And

50:54 we plug in the function value. is one of all the different arguments

50:59 extra asbestos function values of three points one. Such especially some of the

51:07 . Um text and integrate text that's negative side on the positive side of

51:14 equals zero. Equal areas such as builders, zero. Um The Function

51:22 of X -1 is -1. So . So again Xterra at the midpoint

51:29 from and similarly from values from the value. So now we have three

51:37 for three and also we can figure out what these professionals are and then

51:42 the A. N. C. this sign and they have all of

51:49 . So so this is not what but this centers off it is of

51:57 form and it was designed such that quadratic pulling almost identical except approximation.

52:05 the approximation ever happens between from the in the integration of the pulling

52:14 So then we're just gonna stand into arbitrary point. So the manager gets

52:22 is from A to B. That's expression that looks like this. So

52:28 not two different um the weights are of difference compared to taking other chapters

52:40 there was half of it to the points and then it was just a

52:46 by the live. So the whites is yes ma'am. So the interval

52:54 of length that could have this food being monitored. They provide about

52:59 So this is the right, And then the formula has three.

53:03 that's fun. The sixth. and this formula is known as same

53:13 alone and that. So it's still simple, but it turns off go

53:23 inside somewhere, but there's an expense me order to report, but it's

53:31 ah more accurate than the chapters and not very much more complex. It's

53:39 , but my phone. So then the very simple mitigation format. But

53:46 advice you seem because it's a lot than, and it's not for

53:53 I know. So that's an example where things. That's right. It's

54:03 very subtle symptoms, picks up a bit more about courage and about and

54:13 we have against the americans, some too. So, So we have

54:19 investigation only six square which arm with and and see what happens.

54:30 ah, so I just standard into 0 to 1, do the very

54:37 one per episode for the chapters or . And that's just basically two endpoints

54:47 provided by harsh and that's what you now on the steps of the soups

54:52 . We also have the net point the proper waiting years and I will

54:56 the value that is quite different from use this pair of function that turns

55:04 to be the function that this integral . And then you can see that

55:10 is actually okay, they're much better and the practice now that wasn't quite

55:21 comparison I would say because and the rule will use three function evaluations.

55:30 we did one more function evaluations so would say. So this to be

55:36 fair one should have used to track . They especially have the same number

55:43 function evaluations because that is more apparent me. So I think that you

55:50 do that in the book but so you actually use um the trap is

55:56 instead of one same number of function but trappers or rules just to get

56:01 better result from this, you get but this is still much so that

56:08 that Yeah bye. Using a different of the three function values. Ah

56:19 is in the Simpson waking of the that is compared to if you use

56:25 chapters always here then you still have of these two values and then to

56:30 not before but the function value of meatball there, they still don't get

56:38 good an approximation. So that would Yes for 1/2 in the mid

56:46 And whereas these days So the 16 one point and to thirds of 46

56:54 the leg weights the hit points and differently then twice chapter one or to

57:02 results. So it's the waiting that makes a difference. Yeah, so

57:12 guess it's a than looking at so it took this formulation and trying to

57:22 interpolation polynomial Now if you have three and I choose basically 40 and age

57:32 age 2 3 points. So it's quite plus minus H No one has

57:38 But from 0 to 2 age then got this cardinal Obando meals that has

57:46 length of the For the distance between three points age for a pair of

57:53 and there are two main towards so are the cardinal polynomial. So then

57:58 we integrate the function values as you , they don't change the X.

58:04 are given for the global points. it's only the cardinal polynomial surfers.

58:11 variable X and y integrate this place this for the gap. Okay.

58:19 okay. So yes, yes, the function values. And then to

58:24 integration to work it out and that's this is going to work it out

58:30 somewhere along the lines. Eventually you this form of your simplified so so

58:37 the third one can start in this . It was convenient to use from

58:42 ground formulation because yeah, um the polynomial has nothing to do with the

58:50 that so that this was just a separation and you can find the weights

58:56 this formula for the integration that we or a B and C. S

59:05 waste for the function graduates. They out directly from the integration or

59:11 So we got these tractors often. this also shows that the symptoms exactly

59:25 greats. The director. Mm S So this is a good let's

59:36 about any questions so far. So also different integration formulas out there.

59:46 , talk about it. Also expect and sometimes there's no less quality moves

59:52 it's all place around with the white , different function values in the

59:58 The intensity. Thanks. Mm So now I would like to show

60:04 the accuracy is of the same sense all. And I'm kind of on

60:10 stage four and we'll see if I'm . Um, so again, go

60:14 to this simple taylor series expansion for same role. We have The 300

60:23 from there. They point and they for us At the midpoint is four

60:28 the endpoint is one and divided. you take the steel suspension or functional

60:38 passage to age. Then they tried walk down the expression or the right

60:46 side of the topic. Creation. think. So now they want the

60:51 hand side that is up there, fire. But it's just having these

61:00 equations Plus this is more than proper . So understanding, not multiply this

61:08 and this 54 and the system So straightforward competition of They were serious

61:23 this expression and an expression of this . And what it looks like.

61:30 not done yet more. Now we've something like this. So now we

61:40 the right hand side. Now we look at what the left hand side

61:45 like for that. Ah, they this form on the color. I

61:54 the answer derivative, the capital And then they're going to make it

62:00 okay, two Plus H version. , so we still have a small

62:06 to roll then they can just stay same suspension of the capital and

62:15 this form of expression and then we that um, integrate from 8 to

62:19 0. So this is double zero derivative or capital F is lower case

62:27 . So, And then it's There's lower case for one. Go

62:33 . The insurance. Mm hmm. of this was coming. So

62:45 we have the So this is the A. two H. So this

62:51 kind of this expression And um, what putting in the fact that this

62:59 here and this is that's the workplace . So we best support to

63:05 So that's just plugging in the So that's how So this is the

63:10 thing of the left hand side. now what is it going to do

63:16 take the difference what they have this inside right here we have expression for

63:26 using table series expansion and then now have it on the left hand side

63:33 both left and right hand side. friends are this thing and the left

63:42 side. And so this is the . And put that here. But

63:53 warrants. Yes. So they have expression that um so anyway, you

64:03 compare terms here hopefully. So these the same and the right left hand

64:09 and that was the right hand This term is also the same and

64:13 is the same. Um so the term that is different is in fact

64:23 . So now I said that in of the integral this morning, Israel

64:29 exactly material Simpson role for the interval an error term is mm hmm.

64:42 , so that was from a small but in general they want to do

64:45 larger interval. And then we apply send some jewel to a bunch of

64:50 access. Oh, it's called. that's why you and am I will

64:57 him the next time. That's all to have a bunch of these.

65:01 one order it disappears because of the of things that are proportional to

65:13 So, so in other words, just passed it. No, something

65:24 is determines the that's a long Um, that doesn't have a

65:33 Sorry. Yes, that's one another five which is, it's one symptom

65:39 whole thing. But you want the , all the things you divide the

65:47 between A and B into some intervals then you apply Simpson's rule for each

65:52 of these sounds into all the same as So let's start from the beginning

66:00 year. Any questions so far on , on this today, I'm showing

66:05 electron for me it's it was the trick that was used plus before and

66:14 ah left hand side. And use anti derivative and the privacy of expansion

66:21 on the left hand side and also a series expansion on the right hand

66:26 , look at the difference and then downwards the first. How about

66:39 wow ! All right now you don't just want to stay in some form

66:46 the interval all the intervals. And the results of the symptoms smaller.

66:59 that divided the interval between D. . A. And to end

67:06 And for each one of these segments going so far simple something. So

67:13 million years with the kind of symmetric I'm now the midpoint for each of

67:19 Simpson's rule. So that's an And then some of you're insulted of

67:26 Simpsons. We're going supply for these in cigarettes. So this is exactly

67:35 the sense of to do a bunch intervals over small segments to age and

67:41 . And we just divide the integral into something else. For informative,

67:46 memory. Very soft for all the , smaller intervals, no approximation stuff

67:54 approximation happens in the interval for the interval. Where are we going to

68:00 this? So so now we're looking during the integration and replacing this uh

68:13 the interval of length to each for Simpson formula which is left standpoint report

68:24 right. So this is basically to -2 and this is why I'm on

68:32 1610 and some of the age the thing is right one of the two

68:43 points of this million to one and in the middle The Vitamin D three

68:48 after transforming it will some interval of function. So now that's why the

68:57 comes in it's in formula but it's now the truant developments the Simpsons for

69:03 smaller so what kind of player all this expression? Ah Now get the

69:17 form. So if you take one these sub intervals and he takes the

69:27 summer vegetables ah The two endpoints of subject of all is in both the

69:36 and right something. So so this is the and for the sake of

69:46 previous interval and this is also right point of the figures and this is

69:52 left hand from from the next thing . So that's why you live with

69:57 guests. Now the waist this four the mid points dr Alexis Simpson install

70:09 they each fill sub interval. The appear fights especially left hand point and

70:19 and right hand front in and So now instead of having the two

70:25 staying on the way The same but it's up to the same chapter Sawyer

70:33 divided by two The Total vol three instead of just something the internal function

70:43 I think it's weighted by 424242 instead 24. So to me it's a

70:49 simple performance also kind of so memorize Simpson worked as opposed to the tractor

70:56 its various simple parents speaking. So Right so because now yeah that's fire

71:18 of this book potentially the error accumulates all the different sub intervals and subject

71:27 and end times pages blank for intervals can produce one power here but they're

71:38 so that's what becomes age to the in there. So so that's solid

71:48 it's for a very simple scheme. on the very simple example of one

71:55 record for the for the virus expert . The difference in waiting for the

72:02 number of function evaluations Stadium that That's fine. Mhm Thanks. Evaluating

72:13 expression instead of just having them waiting or by to market time of 4

72:19 instead of best not to find one us having them up very little.

72:30 are you doing? Multiplying by two forces shifts amongst by in terms of

72:36 binary systems. Very simple. You that cell free. Right. So

72:46 can see not only is oh jihadis four instead of those two possible proficiency

72:55 the 15 smaller and smaller. So works out. Okay. Stop

73:08 So um too many great informer's Yeah, that's something else.

73:30 So So we talked about the tractor rule, the basic 1 um two

73:40 we got this for this kind of then Now two steps from to grab

73:48 swords. So the interval 9 15 divided by two. And that they

73:55 for the endpoints and and then we the new points and this is the

74:02 . And if we then your insides users hamburg algorithm based on the

74:12 So now what happened? two points H&H of the two With a step

74:19 one and 2 trapezoid And then we do one step the Richardson extrapolation of

74:28 step of the Romberg. They use to to compute an approach estimate.

74:36 if you manipulate these two expressions you discover the actual against things. So

74:42 noted its exact same. Waiting other values best to forget if you use

74:50 a fun mix population on the so this also shows up and they have

75:01 but but these weights find out something adequate to they should before the same

75:09 up as Bloomberg God it's all the . But now we're kind of awesome

75:21 for starting for deriving the 60 based for normal skin. That's actually that's

75:29 into the second order, There's no instead of just looking at Stromberg algorithm

75:41 expressions and yes, you know the term. But driving the waiting positions

75:50 for the function values based on what normal to want to integrate correctly.

75:57 gives me a better intuition for what . It does a good job on

76:07 same thing on the market. So not really. So that was that's

76:19 think for today so we'll continue to a little bit more different and progression

76:27 us. It's time. Yes, had to come and so on.

76:35 back to that. If we talk compliance in other words, you're looking

76:42 this the chapters old rule. Integrated voice senior functions correctly. The Simpson

76:51 basically into it functions correctly which consists peaceful, he's still in the

77:03 If you think of it as the is successful the next interval and have

77:09 totally different political enormous success. The of political enormous because that happens to

77:14 your function then it's integrative. It rain. Ah It's obviously it works

77:27 , it doesn't mean that it's a flexible in terms of what options from

77:31 integrated function for the normal instead of single and it still doesn't Okay.

77:50 questions. Oh yes.

-
+