1 00:00:00 --> 00:00:06 2 00:00:09 --> 00:00:15 This time, we started solving differential equations. 3 00:00:13 --> 00:00:19 This is the third lecture of the term, and I have yet to 4 00:00:18 --> 00:00:24 solve a single differential equation in this class. 5 00:00:22 --> 00:00:28 Well, that will be rectified from now until the end of the 6 00:00:27 --> 00:00:33 term. So, once you learn separation 7 00:00:30 --> 00:00:36 of variables, which is the most elementary 8 00:00:34 --> 00:00:40 method there is, the single, I think the single 9 00:00:38 --> 00:00:44 most important equation is the one that's called the first 10 00:00:43 --> 00:00:49 order linear equation, both because it occurs 11 00:00:47 --> 00:00:53 frequently in models because it's solvable, 12 00:00:50 --> 00:00:56 and-- I think that's enough. If you drop the course after 13 00:00:56 --> 00:01:02 today you will still have learned those two important 14 00:00:59 --> 00:01:05 methods: separation of variables, and first order 15 00:01:02 --> 00:01:08 linear equations. So, what does such an equation 16 00:01:06 --> 00:01:12 look like? Well, I'll write it in there. 17 00:01:09 --> 00:01:15 There are several ways of writing it, but I think the most 18 00:01:13 --> 00:01:19 basic is this. I'm going to use x as the 19 00:01:16 --> 00:01:22 independent variable because that's what your book does. 20 00:01:20 --> 00:01:26 But in the applications, it's often t, 21 00:01:23 --> 00:01:29 time, that is the independent variable. 22 00:01:25 --> 00:01:31 And, I'll try to give you examples which show that. 23 00:01:29 --> 00:01:35 So, the equation looks like this. 24 00:01:31 --> 00:01:37 I'll find some function of x times y prime plus some other 25 00:01:36 --> 00:01:42 function of x times y is equal to yet another function of x. 26 00:01:40 --> 00:01:46 Obviously, the x doesn't have 27 00:01:46 --> 00:01:52 the same status here that y does, so y is extremely limited 28 00:01:50 --> 00:01:56 in how it can appear in the equation. 29 00:01:53 --> 00:01:59 But, x can be pretty much arbitrary in those places. 30 00:01:57 --> 00:02:03 So, that's the equation we are talking about, 31 00:02:00 --> 00:02:06 and I'll put it up. This is the first version of 32 00:02:05 --> 00:02:11 it, and we'll call them purple. Now, why is that called the 33 00:02:10 --> 00:02:16 linear equation? The word linear is a very 34 00:02:14 --> 00:02:20 heavily used word in mathematics, science, 35 00:02:18 --> 00:02:24 and engineering. For the moment, 36 00:02:20 --> 00:02:26 the best simple answer is because it's linear in y and y 37 00:02:25 --> 00:02:31 prime, the variables y and y prime. 38 00:02:31 --> 00:02:37 Well, y prime is not a variable. 39 00:02:33 --> 00:02:39 Well, you will learn, in a certain sense, 40 00:02:36 --> 00:02:42 it helps to think of it as one, not right now perhaps, 41 00:02:41 --> 00:02:47 but think of it as linear. The most closely analogous 42 00:02:45 --> 00:02:51 thing would be a linear equation, a real linear 43 00:02:49 --> 00:02:55 equation, the kind you studied in high school, 44 00:02:53 --> 00:02:59 which would look like this. It would have two variables, 45 00:02:58 --> 00:03:04 and, I guess, constant coefficients, 46 00:03:00 --> 00:03:06 equal c. Now, that's a linear equation. 47 00:03:05 --> 00:03:11 And that's the sense in which this is linear. 48 00:03:09 --> 00:03:15 It's linear in y prime and y, which are the analogs of the 49 00:03:14 --> 00:03:20 variables y1 and y2. A little bit of terminology, 50 00:03:18 --> 00:03:24 if c is equal to zero, it's called homogeneous, 51 00:03:22 --> 00:03:28 the same way this equation is called homogeneous, 52 00:03:26 --> 00:03:32 as you know from 18.02, if the right-hand side is zero. 53 00:03:32 --> 00:03:38 So, c of x I should write here, but I won't. 54 00:03:37 --> 00:03:43 That's called homogeneous. Now, this is a common form for 55 00:03:42 --> 00:03:48 the equation, but it's not what it's called 56 00:03:46 --> 00:03:52 standard form. The standard form for the 57 00:03:50 --> 00:03:56 equation, and since this is going to be a prime course of 58 00:03:55 --> 00:04:01 confusion, which is probably completely correct, 59 00:04:00 --> 00:04:06 but a prime source of confusion is what I meant. 60 00:04:06 --> 00:04:12 The standard linear form, and I'll underline linear is 61 00:04:11 --> 00:04:17 the first co-efficient of y prime is taken to be one. 62 00:04:16 --> 00:04:22 So, you can always convert that to a standard form by simply 63 00:04:21 --> 00:04:27 dividing through by it. And if I do that, 64 00:04:25 --> 00:04:31 the equation will look like y prime plus, now, 65 00:04:29 --> 00:04:35 it's common to not call it b anymore, the coefficient, 66 00:04:35 --> 00:04:41 because it's really b over a. And, therefore, 67 00:04:40 --> 00:04:46 it's common to adopt, yet, a new letter for it. 68 00:04:44 --> 00:04:50 And, the standard one that many people use is p. 69 00:04:48 --> 00:04:54 How about the right-hand side? We needed a letter for that, 70 00:04:53 --> 00:04:59 too. It's c over a, 71 00:04:54 --> 00:05:00 but we'll call it q. So, when I talk about the 72 00:04:58 --> 00:05:04 standard linear form for a linear first order equation, 73 00:05:03 --> 00:05:09 it's absolutely that that I'm talking about. 74 00:05:08 --> 00:05:14 Now, you immediately see that there is a potential for 75 00:05:13 --> 00:05:19 confusion here because what did I call the standard form for a 76 00:05:20 --> 00:05:26 first-order equation? So, I'm going to say, 77 00:05:25 --> 00:05:31 not this. The standard first order form, 78 00:05:29 --> 00:05:35 what would that be? Well, it would be y prime 79 00:05:35 --> 00:05:41 equals, and everything else on the left-hand side. 80 00:05:39 --> 00:05:45 So, it would be y prime. And now, if I turn this into 81 00:05:44 --> 00:05:50 the standard first-order form, it would be negative p of x y 82 00:05:50 --> 00:05:56 plus q of x. 83 00:05:53 --> 00:05:59 But, of course, nobody would write negative p 84 00:05:57 --> 00:06:03 of x. So, now, I explicitly want to 85 00:06:00 --> 00:06:06 say that this is a form which I will never use for this 86 00:06:05 --> 00:06:11 equation, although half the books of the world do. 87 00:06:12 --> 00:06:18 In short, this poor little first-order equation belongs to 88 00:06:15 --> 00:06:21 two ethnic groups. It's both a first order 89 00:06:18 --> 00:06:24 equation, and therefore, its standard form should be 90 00:06:21 --> 00:06:27 written this way, but it's also a linear 91 00:06:24 --> 00:06:30 equation, and therefore its standard form should be used 92 00:06:28 --> 00:06:34 this way. Well, it has to decide, 93 00:06:30 --> 00:06:36 and I have decided for it. It is, above all, 94 00:06:34 --> 00:06:40 a linear equation, not just a first-order 95 00:06:37 --> 00:06:43 equation. And, in this course, 96 00:06:40 --> 00:06:46 this will always be the standard form. 97 00:06:43 --> 00:06:49 Now, well, what on earth is the difference? 98 00:06:46 --> 00:06:52 If you don't do it that way, the difference is entirely in 99 00:06:51 --> 00:06:57 the sin(p). But, if you get the sign of p 100 00:06:55 --> 00:07:01 wrong in the answers, it is just a disaster from that 101 00:07:00 --> 00:07:06 point on. A trivial little change of sign 102 00:07:04 --> 00:07:10 in the answer produces solutions and functions which have totally 103 00:07:09 --> 00:07:15 different behavior. And, you are going to be really 104 00:07:13 --> 00:07:19 lost in this course. So, maybe I should draw a line 105 00:07:17 --> 00:07:23 through it to indicate, please don't pay any attention 106 00:07:21 --> 00:07:27 to this whatsoever, except that we are not going to 107 00:07:25 --> 00:07:31 do that. Okay, well, what's so important 108 00:07:28 --> 00:07:34 about this equation? Well, number one, 109 00:07:32 --> 00:07:38 it can always be solved. That's a very, 110 00:07:36 --> 00:07:42 very big thing in differential equations. 111 00:07:39 --> 00:07:45 But, it's also the equation which arises in a variety of 112 00:07:44 --> 00:07:50 models. Now, I'm just going to list a 113 00:07:48 --> 00:07:54 few of them. All of them I think you will 114 00:07:51 --> 00:07:57 need either in part one or part two of problem sets over these 115 00:07:57 --> 00:08:03 first couple of problem sets, or second and third maybe. 116 00:08:04 --> 00:08:10 But, of them, I'm going to put at the very 117 00:08:09 --> 00:08:15 top of the list of what I'll call here, I'll give it two 118 00:08:16 --> 00:08:22 names: the temperature diffusion model, well, it would be better 119 00:08:24 --> 00:08:30 to call it temperature concentration by analogy, 120 00:08:29 --> 00:08:35 temperature concentration model. 121 00:08:35 --> 00:08:41 There's the mixing model, which is hardly less important. 122 00:08:40 --> 00:08:46 In other words, it's almost as important. 123 00:08:43 --> 00:08:49 You have that in your problem set. 124 00:08:46 --> 00:08:52 And then, there are other, slightly less important models. 125 00:08:52 --> 00:08:58 There is the model of radioactive decay. 126 00:08:55 --> 00:09:01 There's the model of a bank interest, bank account, 127 00:09:00 --> 00:09:06 various motion models, you know, Newton's Law type 128 00:09:04 --> 00:09:10 problems if you can figure out a way of getting rid of the second 129 00:09:10 --> 00:09:16 derivative, some motion problems. 130 00:09:15 --> 00:09:21 A classic example is the motion of a rocket being fired off, 131 00:09:23 --> 00:09:29 etc., etc., etc. Now, today I have to pick a 132 00:09:28 --> 00:09:34 model. And, the one I'm going to pick 133 00:09:33 --> 00:09:39 is this temperature concentration model. 134 00:09:40 --> 00:09:46 So, this is going to be today's model. 135 00:09:43 --> 00:09:49 Tomorrow's model in the recitation, I'm asking the 136 00:09:47 --> 00:09:53 recitations to, among other things, 137 00:09:50 --> 00:09:56 make sure they do a mixing problem, A) to show you how to 138 00:09:55 --> 00:10:01 do it, and B) because it's on the problem sets. 139 00:10:00 --> 00:10:06 That's not a good reason, but it's not a bad one. 140 00:10:04 --> 00:10:10 The others are either in part one or we will take them up 141 00:10:09 --> 00:10:15 later in the term. This is not going to be the 142 00:10:13 --> 00:10:19 only lecture on the linear equation. 143 00:10:16 --> 00:10:22 There will be another one next week of equal importance. 144 00:10:21 --> 00:10:27 But, we can't do everything today. 145 00:10:24 --> 00:10:30 So, let's talk about the temperature concentration model, 146 00:10:28 --> 00:10:34 except I'm going to change its name. 147 00:10:32 --> 00:10:38 I'm going to change its name to the conduction diffusion model. 148 00:10:39 --> 00:10:45 I'll put conduction over there, and diffusion over here, 149 00:10:44 --> 00:10:50 let's say, since, as you will see, 150 00:10:47 --> 00:10:53 the similarities, they are practically the same 151 00:10:52 --> 00:10:58 model. All that's changed from one to 152 00:10:55 --> 00:11:01 the other is the name of the ideas. 153 00:11:00 --> 00:11:06 In one case, you call it temperature, 154 00:11:03 --> 00:11:09 and the other, you should call it 155 00:11:06 --> 00:11:12 concentration. But, the actual mathematics 156 00:11:10 --> 00:11:16 isn't identical. So, let's begin with 157 00:11:14 --> 00:11:20 conduction. All right, so, 158 00:11:16 --> 00:11:22 I need a simple physical situation that I'm modeling. 159 00:11:21 --> 00:11:27 So, imagine a tank of some liquid. 160 00:11:24 --> 00:11:30 Water will do as well as anything. 161 00:11:27 --> 00:11:33 And, in the inside is a suspended, somehow, 162 00:11:31 --> 00:11:37 is a chamber. A metal cube will do, 163 00:11:37 --> 00:11:43 and let's suppose that its walls are partly insulated, 164 00:11:42 --> 00:11:48 not so much that no heat can get through. 165 00:11:46 --> 00:11:52 There is no such thing as perfect insulation anyway, 166 00:11:51 --> 00:11:57 except maybe an absolute perfect vacuum. 167 00:11:56 --> 00:12:02 Now, inside, so here on the outside is 168 00:11:59 --> 00:12:05 liquid. Okay, on the inside is, 169 00:12:02 --> 00:12:08 what I'm interested in is the temperature of this thing. 170 00:12:10 --> 00:12:16 I'll call that T. Now, that's different from the 171 00:12:13 --> 00:12:19 temperature of the external water bath. 172 00:12:17 --> 00:12:23 So, I'll call that T sub e, T for temperature 173 00:12:21 --> 00:12:27 measured in Celsius, let's say, for the sake of 174 00:12:25 --> 00:12:31 definiteness. But, this is the external 175 00:12:28 --> 00:12:34 temperature. So, I'll indicate it with an e. 176 00:12:33 --> 00:12:39 Now, what is the model? Well, in other words, 177 00:12:38 --> 00:12:44 how do I set up a differential equation to model the situation? 178 00:12:46 --> 00:12:52 Well, it's based on a physical law, which I think you know, 179 00:12:53 --> 00:12:59 you've had simple examples like this, the so-called Newton's Law 180 00:13:01 --> 00:13:07 of cooling, -- -- which says that the rate of 181 00:13:06 --> 00:13:12 change, the temperature of the heat goes from the outside to 182 00:13:12 --> 00:13:18 the inside by conduction only. Heat, of course, 183 00:13:16 --> 00:13:22 can travel in various ways, by convection, 184 00:13:19 --> 00:13:25 by conduction, as here, or by radiation, 185 00:13:23 --> 00:13:29 are the three most common. Of these, I only want one, 186 00:13:28 --> 00:13:34 namely transmission of heat by conduction. 187 00:13:33 --> 00:13:39 And, that's the way it's probably a little better to call 188 00:13:37 --> 00:13:43 it the conduction model, rather than the temperature 189 00:13:41 --> 00:13:47 model, which might involve other ways for the heat to be 190 00:13:46 --> 00:13:52 traveling. So, dt, the independent 191 00:13:49 --> 00:13:55 variable, is not going to be x, as it was over there. 192 00:13:53 --> 00:13:59 It's going to be t for time. So, maybe I should write that 193 00:13:58 --> 00:14:04 down. t equals time. 194 00:13:59 --> 00:14:05 Capital T equals temperature in degrees Celsius. 195 00:14:05 --> 00:14:11 So, you can put in the degrees Celsius if you want. 196 00:14:08 --> 00:14:14 So, it's proportional to the temperature difference between 197 00:14:13 --> 00:14:19 these two. Now, how shall I write the 198 00:14:16 --> 00:14:22 difference? Write it this way because if 199 00:14:19 --> 00:14:25 you don't you will be in trouble. 200 00:14:21 --> 00:14:27 Now, why do I write it that way? 201 00:14:24 --> 00:14:30 Well, I write it that way because I want this constant to 202 00:14:28 --> 00:14:34 be positive, a positive constant. 203 00:14:32 --> 00:14:38 In general, any constant, so, parameters which are 204 00:14:35 --> 00:14:41 physical, have some physical significance, 205 00:14:39 --> 00:14:45 one always wants to arrange the equation so that they are 206 00:14:43 --> 00:14:49 positive numbers, the way people normally think 207 00:14:47 --> 00:14:53 of these things. This is called the 208 00:14:49 --> 00:14:55 conductivity. The conductivity of what? 209 00:14:52 --> 00:14:58 Well, I don't know, of the system of the situation, 210 00:14:56 --> 00:15:02 the conductivity of the wall, or the wall if the metal were 211 00:15:01 --> 00:15:07 just by itself. At any rate, 212 00:15:03 --> 00:15:09 it's a constant. It's thought of as a constant. 213 00:15:07 --> 00:15:13 And, y positive, well, because if the external 214 00:15:10 --> 00:15:16 temperature is bigger than the internal temperature, 215 00:15:14 --> 00:15:20 I expect T to rise, the internal temperature to 216 00:15:18 --> 00:15:24 rise. That means dT / dt, 217 00:15:22 --> 00:15:28 its slope, should be positive. So, in other words, 218 00:15:26 --> 00:15:32 if Te is bigger than T, I expect this number to be 219 00:15:29 --> 00:15:35 positive. And, that tells you that k must 220 00:15:32 --> 00:15:38 be a positive constant. If I had turned it the other 221 00:15:36 --> 00:15:42 way, expressed the difference in the reverse order, 222 00:15:40 --> 00:15:46 K would then be negative, have to be negative in order 223 00:15:44 --> 00:15:50 that this turn out to be positive in that situation I 224 00:15:47 --> 00:15:53 described. And, since nobody wants 225 00:15:50 --> 00:15:56 negative values of k, you have to write the equation 226 00:15:53 --> 00:15:59 in this form rather than the other way around. 227 00:15:56 --> 00:16:02 So, there's our differential equation. 228 00:15:59 --> 00:16:05 It will probably have an initial condition. 229 00:16:02 --> 00:16:08 So, it could be the temperature at the starting time should be 230 00:16:06 --> 00:16:12 some given number, T zero. 231 00:16:10 --> 00:16:16 But, the condition could be given in other ways. 232 00:16:13 --> 00:16:19 One can ask, what's the temperature as time 233 00:16:15 --> 00:16:21 goes to infinity, for example? 234 00:16:17 --> 00:16:23 There are different ways of getting that initial condition. 235 00:16:21 --> 00:16:27 Okay, that's the conduction model. 236 00:16:23 --> 00:16:29 What would the diffusion model be? 237 00:16:25 --> 00:16:31 The diffusion model, mathematically, 238 00:16:27 --> 00:16:33 would be, word for word, the same. 239 00:16:31 --> 00:16:37 The only difference is that now, what I imagine is I'll draw 240 00:16:36 --> 00:16:42 the picture the same way, except now I'm going to put, 241 00:16:41 --> 00:16:47 label the inside not with a T but with a C, 242 00:16:45 --> 00:16:51 C for concentration. It's in an external water bath, 243 00:16:50 --> 00:16:56 let's say. So, there is an external 244 00:16:53 --> 00:16:59 concentration. And, what I'm talking about is 245 00:16:57 --> 00:17:03 some chemical, let's say salt will do as well 246 00:17:01 --> 00:17:07 as anything. So, C is equal to salt 247 00:17:06 --> 00:17:12 concentration inside, and Ce would be the salt 248 00:17:12 --> 00:17:18 concentration outside, outside in the water bath. 249 00:17:18 --> 00:17:24 Now, I imagine some mechanism, so this is a salt solution. 250 00:17:26 --> 00:17:32 That's a salt solution. And, I imagine some mechanism 251 00:17:32 --> 00:17:38 by which the salt can diffuse, it's a diffusion model now, 252 00:17:36 --> 00:17:42 diffuse from here into the air or possibly out the other way. 253 00:17:40 --> 00:17:46 And that's usually done by vaguely referring to the outside 254 00:17:45 --> 00:17:51 as a semi-permeable membrane, semi-permeable, 255 00:17:48 --> 00:17:54 so that the salt will have a little hard time getting through 256 00:17:53 --> 00:17:59 but permeable, so that it won't be blocked 257 00:17:56 --> 00:18:02 completely. So, there's a membrane. 258 00:18:00 --> 00:18:06 You write the semi-permeable membrane outside, 259 00:18:06 --> 00:18:12 outside the inside. Well, I give up. 260 00:18:10 --> 00:18:16 You know, membrane somewhere. Sorry, membrane wall. 261 00:18:17 --> 00:18:23 How's that? Now, what's the equation? 262 00:18:21 --> 00:18:27 Well, the equation is the same, except it's called the 263 00:18:28 --> 00:18:34 diffusion equation. I don't think Newton got his 264 00:18:34 --> 00:18:40 name on this. The diffusion equation says 265 00:18:37 --> 00:18:43 that the rate at which the salt diffuses across the membrane, 266 00:18:42 --> 00:18:48 which is the same up to a constant as the rate at which 267 00:18:47 --> 00:18:53 the concentration inside changes, is some constant, 268 00:18:51 --> 00:18:57 usually called k still, okay. 269 00:18:54 --> 00:19:00 Do I contradict? Okay, let's keep calling it k1. 270 00:18:58 --> 00:19:04 Now it's different, times Ce minus C. 271 00:19:01 --> 00:19:07 And, for the same reason as before, if the external 272 00:19:05 --> 00:19:11 concentration is bigger than the internal concentration, 273 00:19:10 --> 00:19:16 we expect salt to flow in. That will make C rise. 274 00:19:15 --> 00:19:21 It will make this positive, and therefore, 275 00:19:18 --> 00:19:24 we want k to be positive, just k1 to be positive for the 276 00:19:22 --> 00:19:28 same reason it had to be positive before. 277 00:19:25 --> 00:19:31 So, in each case, the model that I'm talking 278 00:19:28 --> 00:19:34 about is the differential equation. 279 00:19:31 --> 00:19:37 So, maybe I should, let's put that, 280 00:19:33 --> 00:19:39 make that clear. Or, I would say that this first 281 00:19:38 --> 00:19:44 order differential equation models this physical situation, 282 00:19:42 --> 00:19:48 and the same thing is true on the other side over here. 283 00:19:46 --> 00:19:52 This is the diffusion equation, and this is the conduction 284 00:19:50 --> 00:19:56 equation. Now, if you are in any doubt 285 00:19:53 --> 00:19:59 about the power of differential equations, the point is, 286 00:19:57 --> 00:20:03 when I talk about this thing, I don't have to say which of 287 00:20:01 --> 00:20:07 these I'm following. I'll use neutral variables like 288 00:20:06 --> 00:20:12 Y and X to solve these equations. 289 00:20:09 --> 00:20:15 But, with a single stroke, I will be handling those 290 00:20:13 --> 00:20:19 situations together. And, that's the power of the 291 00:20:16 --> 00:20:22 method. Now, you obviously must be 292 00:20:19 --> 00:20:25 wondering, look, these look very, 293 00:20:22 --> 00:20:28 very special. He said he was going to talk 294 00:20:25 --> 00:20:31 about the first, general first-order equation. 295 00:20:28 --> 00:20:34 But, these look rather special to me. 296 00:20:31 --> 00:20:37 Well, not too special. How should we write it? 297 00:20:36 --> 00:20:42 Suppose I write, let's take the temperature 298 00:20:40 --> 00:20:46 equation just to have something definite. 299 00:20:44 --> 00:20:50 Notice that it's in a form corresponding to Newton's Law. 300 00:20:48 --> 00:20:54 But it is not in the standard linear form. 301 00:20:52 --> 00:20:58 Let's put it in standard linear form, so at least you could see 302 00:20:57 --> 00:21:03 that it's a linear equation. So, if I put it in standard 303 00:21:02 --> 00:21:08 form, it's going to look like DTDTD little t plus KT is equal 304 00:21:07 --> 00:21:13 to K times TE. Now, compare that with the 305 00:21:12 --> 00:21:18 general, the way the general equation is supposed to look, 306 00:21:16 --> 00:21:22 the yellow box over there, the standard linear form. 307 00:21:20 --> 00:21:26 How are they going to compare? Well, this is a pretty general 308 00:21:24 --> 00:21:30 function. This is general. 309 00:21:26 --> 00:21:32 This is a general function of T because I can make the external 310 00:21:31 --> 00:21:37 temperature. I could suppose it behaves in 311 00:21:34 --> 00:21:40 anyway I like, steadily rising, 312 00:21:37 --> 00:21:43 decaying exponentially, maybe oscillating back and 313 00:21:40 --> 00:21:46 forth for some reason. The only way in which it's not 314 00:21:46 --> 00:21:52 general is that this K is a constant. 315 00:21:50 --> 00:21:56 So, I will ask you to be generous. 316 00:21:53 --> 00:21:59 Let's imagine the conductivity is changing over time. 317 00:21:57 --> 00:22:03 So, this is usually constant, but there's no law which says 318 00:22:03 --> 00:22:09 it has to be. How could a conductivity change 319 00:22:08 --> 00:22:14 over time? Well, we could suppose that 320 00:22:12 --> 00:22:18 this wall was made of slowly congealing Jell-O, 321 00:22:16 --> 00:22:22 for instance. It starts out as liquid, 322 00:22:20 --> 00:22:26 and then it gets solid. And, Jell-O doesn't transmit 323 00:22:26 --> 00:22:32 heat, I believe, quite as well as liquid does, 324 00:22:30 --> 00:22:36 as a liquid would. Is Jell-O a solid or liquid? 325 00:22:36 --> 00:22:42 I don't know. Let's forget about that. 326 00:22:40 --> 00:22:46 So, with this understanding, so let's say not necessarily 327 00:22:46 --> 00:22:52 here, but not necessarily, I can think of this, 328 00:22:51 --> 00:22:57 therefore, by allowing K to vary with time. 329 00:22:55 --> 00:23:01 And the external temperature to vary with time. 330 00:23:00 --> 00:23:06 I can think of it as a general, linear equation. 331 00:23:07 --> 00:23:13 So, these models are not special. 332 00:23:09 --> 00:23:15 They are fairly general. Well, I did promise you I would 333 00:23:13 --> 00:23:19 solve an equation, and that this lecture, 334 00:23:16 --> 00:23:22 I still have not solved any equations. 335 00:23:18 --> 00:23:24 OK, time to stop temporizing and solve. 336 00:23:21 --> 00:23:27 So, I'm going to, in order not to play favorites 337 00:23:24 --> 00:23:30 with these two models, I'll go back to, 338 00:23:27 --> 00:23:33 and to get you used to thinking of the variables all the time, 339 00:23:31 --> 00:23:37 that is, you know, be eclectic switching from one 340 00:23:35 --> 00:23:41 variable to another according to which particular lecture you 341 00:23:39 --> 00:23:45 happened to be sitting in. So, let's take our equation in 342 00:23:47 --> 00:23:53 the form, Y prime plus P of XY, the general form using the old 343 00:23:57 --> 00:24:03 variables equals Q of X. Solve me. 344 00:24:04 --> 00:24:10 Well, there are different ways of describing the solution 345 00:24:07 --> 00:24:13 process. No matter how you do it, 346 00:24:10 --> 00:24:16 it amounts to the same amount of work and there is always a 347 00:24:14 --> 00:24:20 trick involved at each one of them since you can't suppress a 348 00:24:18 --> 00:24:24 trick by doing the problem some other way. 349 00:24:21 --> 00:24:27 The way I'm going to do it, I think, is the best. 350 00:24:24 --> 00:24:30 That's why I'm giving it to you. 351 00:24:26 --> 00:24:32 It's the easiest to remember. It leads to the least work, 352 00:24:30 --> 00:24:36 but I have colleagues who would fight with me about that point. 353 00:24:36 --> 00:24:42 So, since they are not here to fight with me I am free to do 354 00:24:41 --> 00:24:47 whatever I like. One of the main reasons for 355 00:24:45 --> 00:24:51 doing it the way I'm going to do is because I want you to get 356 00:24:51 --> 00:24:57 what our word into your consciousness, 357 00:24:55 --> 00:25:01 two words, integrating factor. I'm going to solve this 358 00:25:00 --> 00:25:06 equation by finding and integrating factor of the form U 359 00:25:05 --> 00:25:11 of X. What's an integrating factor? 360 00:25:09 --> 00:25:15 Well, I'll show you not by writing an elaborate definition 361 00:25:13 --> 00:25:19 on the board, but showing you what its 362 00:25:16 --> 00:25:22 function is. It's a certain function, 363 00:25:18 --> 00:25:24 U of X, I don't know what it is, but here's what I wanted to 364 00:25:23 --> 00:25:29 do. I want to multiply, 365 00:25:24 --> 00:25:30 I'm going to drop the X's a just so that the thing looks 366 00:25:28 --> 00:25:34 less complicated. So, what I want to do is 367 00:25:33 --> 00:25:39 multiply this equation through by U of X. 368 00:25:36 --> 00:25:42 That's why it's called a factor because you're going to multiply 369 00:25:42 --> 00:25:48 everything through by it. So, it's going to look like UY 370 00:25:47 --> 00:25:53 prime plus PUY equals QU, and now, so far, 371 00:25:51 --> 00:25:57 it's just a factor. What makes it an integrating 372 00:25:55 --> 00:26:01 factor is that this, after I do that, 373 00:25:58 --> 00:26:04 I want this to turn out to be the derivative of something with 374 00:26:04 --> 00:26:10 respect to X. You see the motivation for 375 00:26:08 --> 00:26:14 that. If this turns out to be the 376 00:26:10 --> 00:26:16 derivative of something, because I've chosen U so 377 00:26:13 --> 00:26:19 cleverly, then I will be able to solve the equation immediately 378 00:26:17 --> 00:26:23 just by integrating this with respect to X, 379 00:26:19 --> 00:26:25 and integrating that with respect to X. 380 00:26:21 --> 00:26:27 You just, then, integrate both sides with 381 00:26:24 --> 00:26:30 respect to X, and the equation is solved. 382 00:26:26 --> 00:26:32 Now, the only question is, what should I choose for U? 383 00:26:31 --> 00:26:37 Well, if you think of the product formula, 384 00:26:34 --> 00:26:40 there might be many things to try here. 385 00:26:37 --> 00:26:43 But there's only one reasonable thing to try. 386 00:26:40 --> 00:26:46 Try to pick U so that it's the derivative of U times Y. 387 00:26:45 --> 00:26:51 See how reasonable that is? If I use the product rule on 388 00:26:49 --> 00:26:55 this, the first term is U times Y prime. 389 00:26:52 --> 00:26:58 The second term would be U prime times Y. 390 00:26:56 --> 00:27:02 Well, I've got the Y there. So, this will work. 391 00:27:01 --> 00:27:07 It works if, what's the condition that you 392 00:27:05 --> 00:27:11 must satisfy in order for that to be true? 393 00:27:09 --> 00:27:15 Well, it must be that after it to the differentiation, 394 00:27:15 --> 00:27:21 U prime turns out to be P times U. 395 00:27:19 --> 00:27:25 So, is it clear? This is something we want to be 396 00:27:24 --> 00:27:30 equal to, and the thing I will try to do it is by choosing U in 397 00:27:30 --> 00:27:36 such a way that this equality will take place. 398 00:27:37 --> 00:27:43 And then I will be able to solve the equation. 399 00:27:40 --> 00:27:46 And so, here's what my U prime must satisfy. 400 00:27:43 --> 00:27:49 Hey, we can solve that. But please don't forget that P 401 00:27:47 --> 00:27:53 is P of X. It's a function of X. 402 00:27:49 --> 00:27:55 So, if you separate variables, I'm going to do this. 403 00:27:53 --> 00:27:59 So, what is it, DU over U equals P of X times 404 00:27:56 --> 00:28:02 DX. If I integrate that, 405 00:27:59 --> 00:28:05 so, separate variables, integrate, and you're going to 406 00:28:03 --> 00:28:09 get DU over U integrates to the be the log of U, 407 00:28:07 --> 00:28:13 and the other side integrates to be the integral of P of X DX. 408 00:28:12 --> 00:28:18 Now, you can put an arbitrary constant there, 409 00:28:16 --> 00:28:22 or you can think of it as already implied by the 410 00:28:19 --> 00:28:25 indefinite integral. Well, that doesn't tell us, 411 00:28:23 --> 00:28:29 yet, what U is. What should U be? 412 00:28:26 --> 00:28:32 Notice, I don't have to find every possible U, 413 00:28:29 --> 00:28:35 which works. All I'm looking for is one. 414 00:28:34 --> 00:28:40 All I want is a single view which satisfies that equation. 415 00:28:38 --> 00:28:44 Well, U equals the integral, E to the integral of PDX. 416 00:28:42 --> 00:28:48 That's not too beautiful looking, but by differential 417 00:28:46 --> 00:28:52 equations, things can get so complicated that in a week or 418 00:28:50 --> 00:28:56 two, you will think of this as an extremely simple formula. 419 00:28:55 --> 00:29:01 So, there is a formula for our integrating factor. 420 00:29:00 --> 00:29:06 We found it. We will always be able to write 421 00:29:05 --> 00:29:11 an integrating factor. Don't worry about the arbitrary 422 00:29:11 --> 00:29:17 constant because you only need one such U. 423 00:29:17 --> 00:29:23 So: no arbitrary constant since only one U needed. 424 00:29:23 --> 00:29:29 And, that's the solution, the way we solve the linear 425 00:29:29 --> 00:29:35 equation. OK, let's take over, 426 00:29:35 --> 00:29:41 and actually do it. I think it would be better to 427 00:29:42 --> 00:29:48 summarize it as a clear-cut method. 428 00:29:48 --> 00:29:54 So, let's do that. So, what's our method? 429 00:29:54 --> 00:30:00 It's the method for solving Y prime plus PY equals Q. 430 00:30:04 --> 00:30:10 Well, the first place, make sure it's in standard 431 00:30:07 --> 00:30:13 linear form. If it isn't, 432 00:30:09 --> 00:30:15 you must put it in that form. Notice, the formula for the 433 00:30:13 --> 00:30:19 integrating factor, the formula for the integrating 434 00:30:16 --> 00:30:22 factor involves P, the integral of PDX. 435 00:30:19 --> 00:30:25 So, you'd better get the right P. 436 00:30:21 --> 00:30:27 Otherwise, you are sunk. OK, so put it in standard 437 00:30:25 --> 00:30:31 linear form. That way, you will have the 438 00:30:28 --> 00:30:34 right P. Notice that if you wrote it in 439 00:30:32 --> 00:30:38 that form, and all you remembered was E to the integral 440 00:30:35 --> 00:30:41 PDX, the P would have the wrong sign. 441 00:30:38 --> 00:30:44 If you're going to write, that P should have a negative 442 00:30:41 --> 00:30:47 sign there. So, do it this way, 443 00:30:43 --> 00:30:49 and no other way. Otherwise, you will get 444 00:30:46 --> 00:30:52 confused and get wrong signs. And, as I say, 445 00:30:49 --> 00:30:55 that will produce wrong answers, and not just slightly 446 00:30:52 --> 00:30:58 wrong answers, but disastrously wrong answers 447 00:30:55 --> 00:31:01 from the point of view of the modeling if you really want 448 00:30:59 --> 00:31:05 answers to physical problems. So, here's a standard linear 449 00:31:06 --> 00:31:12 form. Then, find the integrating 450 00:31:09 --> 00:31:15 factor. So, calculate E to the 451 00:31:13 --> 00:31:19 integral, PDX, the integrating factor, 452 00:31:17 --> 00:31:23 and that multiply both, I'm putting this as both, 453 00:31:22 --> 00:31:28 underlined that as many times as you have room in your notes. 454 00:31:31 --> 00:31:37 Multiply both sides by this integrating factor by E to the 455 00:31:40 --> 00:31:46 integral PDX. And then, integrate. 456 00:31:46 --> 00:31:52 OK, let's take a simple example. 457 00:31:51 --> 00:31:57 Suppose we started with the equation XY prime minus Y 458 00:31:59 --> 00:32:05 equals, I had X2, X3, something like that, 459 00:32:06 --> 00:32:12 X3, I think, yeah, X2. 460 00:32:12 --> 00:32:18 OK, what's the first thing to do? 461 00:32:16 --> 00:32:22 Put it in standard form. So, step zero will be to write 462 00:32:23 --> 00:32:29 it as Y prime minus one over X times Y equals X2. 463 00:32:30 --> 00:32:36 Let's do the work first, and then I'll talk about 464 00:32:34 --> 00:32:40 mistakes. Well, we now calculate the 465 00:32:39 --> 00:32:45 integrating factor. So, I would do it in steps. 466 00:32:43 --> 00:32:49 You can integrate negative one over X, right? 467 00:32:48 --> 00:32:54 That integrates to minus log X. So, the integrating factor is E 468 00:32:54 --> 00:33:00 to the integral of this, DX. 469 00:32:57 --> 00:33:03 So, it's E to the negative log X. 470 00:33:02 --> 00:33:08 Now, in real life, that's not the way to leave 471 00:33:06 --> 00:33:12 that. What is E to the negative log 472 00:33:10 --> 00:33:16 X? Well, think of it as E to the 473 00:33:13 --> 00:33:19 log X to the minus one. Or, in other words, 474 00:33:18 --> 00:33:24 it is E to the log X is X. So, it's one over X. 475 00:33:23 --> 00:33:29 So, the integrating factor is one over X. 476 00:33:27 --> 00:33:33 OK, multiply both sides by the integrating factor. 477 00:33:34 --> 00:33:40 Both sides of what? Both sides of this: 478 00:33:37 --> 00:33:43 the equation written in standard form, 479 00:33:41 --> 00:33:47 and both sides. So, it's going to be one over 480 00:33:45 --> 00:33:51 XY prime minus one over X2 Y is equal to X2 times one over X, 481 00:33:51 --> 00:33:57 which is simply X. Now, if you have done the work 482 00:33:55 --> 00:34:01 correctly, you should be able, now, to integrate the left-hand 483 00:34:01 --> 00:34:07 side directly. So, I'm going to write it this 484 00:34:06 --> 00:34:12 way. I always recommend that you put 485 00:34:09 --> 00:34:15 it as extra step, well, put it as an extra step 486 00:34:13 --> 00:34:19 the reason for using that integrating factor, 487 00:34:17 --> 00:34:23 in other words, that the left-hand side is 488 00:34:20 --> 00:34:26 supposed to be, now, one over X times Y prime. 489 00:34:24 --> 00:34:30 I always put it that because there's always a chance you made 490 00:34:29 --> 00:34:35 a mistake or forgot something. Look at it, mentally 491 00:34:34 --> 00:34:40 differentiated using the product rule just to check that, 492 00:34:39 --> 00:34:45 in fact, it turns out to be the same as the left-hand side. 493 00:34:43 --> 00:34:49 So, what do we get? One over X times Y prime plus Y 494 00:34:47 --> 00:34:53 times the derivative of one over X, which indeed is negative one 495 00:34:53 --> 00:34:59 over X2. And now, finally, 496 00:34:55 --> 00:35:01 that's 3A, continue, do the integration. 497 00:34:58 --> 00:35:04 So, you're going to get, let's see if we can do it all 498 00:35:02 --> 00:35:08 on one board, one over X times Y is equal to 499 00:35:06 --> 00:35:12 X plus a constant, X, sorry, X2 over two plus a 500 00:35:09 --> 00:35:15 constant. And, the final step will be, 501 00:35:15 --> 00:35:21 therefore, now I want to isolate Y by itself. 502 00:35:21 --> 00:35:27 So, Y will be equal to multiply through by X. 503 00:35:26 --> 00:35:32 X3 over two plus C times X. And, that's the solution. 504 00:35:34 --> 00:35:40 OK, let's do one a little slightly more complicated. 505 00:35:40 --> 00:35:46 Let's try this one. Now, my equation is going to be 506 00:35:45 --> 00:35:51 one, I'll still keep two, Y and X, as the variables. 507 00:35:51 --> 00:35:57 I'll use T and F for a minute or two. 508 00:35:57 --> 00:36:03 One plus cosine X, so, I'm not going to give you 509 00:36:03 --> 00:36:09 this one in standard form either. 510 00:36:07 --> 00:36:13 It's a trick question. Y prime minus sine X times Y is 511 00:36:14 --> 00:36:20 equal to anything reasonable, I guess. 512 00:36:19 --> 00:36:25 I think X, 2X, make it more exciting. 513 00:36:24 --> 00:36:30 OK, now, I think I should warn you where the mistakes are just 514 00:36:32 --> 00:36:38 so that you can make all of them. 515 00:36:38 --> 00:36:44 So, this is mistake number one. You don't put it in standard 516 00:36:44 --> 00:36:50 form. Mistake number two: 517 00:36:47 --> 00:36:53 generally people can do step one fine. 518 00:36:51 --> 00:36:57 Mistake number two is, this is my most common mistake, 519 00:36:57 --> 00:37:03 so I'm very sensitive to it. But that doesn't mean if you 520 00:37:03 --> 00:37:09 make it, you'll get any sympathy from me. 521 00:37:06 --> 00:37:12 I don't give sympathy to myself. 522 00:37:08 --> 00:37:14 You are so intense, so happy at having found the 523 00:37:11 --> 00:37:17 integrating factor, you forget to multiply Q by the 524 00:37:15 --> 00:37:21 integrating factor also. You just handle the left-hand 525 00:37:19 --> 00:37:25 side of the equation, if you forget about the 526 00:37:22 --> 00:37:28 right-hand side. So, the emphasis on the both 527 00:37:25 --> 00:37:31 here is the right-hand, please include the Q. 528 00:37:29 --> 00:37:35 Please include the right-hand side. 529 00:37:33 --> 00:37:39 Any other mistakes? Well, nothing that I can think 530 00:37:37 --> 00:37:43 of. Well, maybe only, 531 00:37:38 --> 00:37:44 anyway, we are not going to make any mistakes the rest of 532 00:37:43 --> 00:37:49 this lecture. So, what do we do? 533 00:37:45 --> 00:37:51 We write this in standard form. So, it's going to look like Y 534 00:37:50 --> 00:37:56 prime minus sine X, sine X divided by one plus 535 00:37:54 --> 00:38:00 cosine X times Y equals, my heart sinks because I know 536 00:37:59 --> 00:38:05 I'm supposed to integrate something like this. 537 00:38:04 --> 00:38:10 And, boy, that's going to give me problems. 538 00:38:07 --> 00:38:13 Well, not yet. With the integrating factor? 539 00:38:11 --> 00:38:17 The integrating factor is, well, we want to calculate the 540 00:38:16 --> 00:38:22 integral of negative sine X over one plus cosine. 541 00:38:20 --> 00:38:26 That's the integral of PDX. And, after that, 542 00:38:23 --> 00:38:29 we have to exponentiate it. Well, can you do this? 543 00:38:28 --> 00:38:34 Yeah, but if you stare at it a little while, 544 00:38:31 --> 00:38:37 you can see that the top is the derivative of the bottom. 545 00:38:38 --> 00:38:44 That is great. That means it integrates to be 546 00:38:42 --> 00:38:48 the log of one plus cosine X. Is that right, 547 00:38:46 --> 00:38:52 one over one plus cosine X times the derivative of this, 548 00:38:51 --> 00:38:57 which is negative cosine X. Therefore, the integrating 549 00:38:56 --> 00:39:02 factor is E to that. In other words, 550 00:38:59 --> 00:39:05 it is one plus cosine X. Therefore, so this was step 551 00:39:05 --> 00:39:11 zero. Step one, we found the 552 00:39:08 --> 00:39:14 integrating factor. And now, step two, 553 00:39:12 --> 00:39:18 we multiply through the integrating factor. 554 00:39:17 --> 00:39:23 And what do we get? We multiply through the 555 00:39:21 --> 00:39:27 standard for equation by the integrating factor, 556 00:39:26 --> 00:39:32 if you do that, what you get is, 557 00:39:29 --> 00:39:35 well, Y prime gets the coefficient one plus cosine X, 558 00:39:35 --> 00:39:41 Y prime minus sign X equals 2X. Oh, dear. 559 00:39:40 --> 00:39:46 Well, I hope somebody would giggle at this point. 560 00:39:45 --> 00:39:51 What's giggle-able about it? Well, that all this was totally 561 00:39:50 --> 00:39:56 wasted work. It's called spinning your 562 00:39:53 --> 00:39:59 wheels. No, it's not spinning your 563 00:39:56 --> 00:40:02 wheels. It's doing what you're supposed 564 00:39:59 --> 00:40:05 to do, and finding out that you wasted the entire time doing 565 00:40:05 --> 00:40:11 what you were supposed to do. Well, in other words, 566 00:40:10 --> 00:40:16 that net effect of this is to end up with the same equation we 567 00:40:16 --> 00:40:22 started with. But, what is the point? 568 00:40:19 --> 00:40:25 The point of having done all this was because now the 569 00:40:24 --> 00:40:30 left-hand side is exactly the derivative of something, 570 00:40:29 --> 00:40:35 and the left-hand side should be the derivative of what? 571 00:40:35 --> 00:40:41 Well, it should be the derivative of one plus cosine X 572 00:40:39 --> 00:40:45 times Y, all prime. Now, you can check that that's 573 00:40:43 --> 00:40:49 in fact the case. It's one plus cosine X, 574 00:40:47 --> 00:40:53 Y prime, plus minus sine X, the derivative of this side 575 00:40:51 --> 00:40:57 times Y. So, if you had thought, 576 00:40:54 --> 00:41:00 in looking at the equation, to say to yourself, 577 00:40:58 --> 00:41:04 this is a derivative of that, maybe I'll just check right 578 00:41:03 --> 00:41:09 away to see if it's the derivative of one plus cosine X 579 00:41:08 --> 00:41:14 sine. You would have saved that work. 580 00:41:12 --> 00:41:18 Well, you don't have to be brilliant or clever, 581 00:41:16 --> 00:41:22 or anything like that. You can follow your nose, 582 00:41:19 --> 00:41:25 and it's just, I want to give you a positive 583 00:41:23 --> 00:41:29 experience in solving linear equations, not too negative. 584 00:41:28 --> 00:41:34 Anyway, so we got to this point. 585 00:41:31 --> 00:41:37 So, now this is 2X, and now we are ready to solve 586 00:41:37 --> 00:41:43 the equation, which is the solution now will 587 00:41:42 --> 00:41:48 be one plus cosine X times Y is equal to X2 plus a constant, 588 00:41:49 --> 00:41:55 and so Y is equal to X2 divided by X2 plus a constant divided by 589 00:41:56 --> 00:42:02 one plus cosine X. Suppose I have given you an 590 00:42:01 --> 00:42:07 initial condition, which I didn't. 591 00:42:03 --> 00:42:09 But, suppose the initial condition said that Y of zero 592 00:42:07 --> 00:42:13 were one, for instance. Then, the solution would be, 593 00:42:10 --> 00:42:16 so, this is an if, I'm throwing in at the end just 594 00:42:14 --> 00:42:20 to make it a little bit more of a problem, how would I put, 595 00:42:17 --> 00:42:23 then I could evaluate the constant by using the initial 596 00:42:21 --> 00:42:27 condition. What would it be? 597 00:42:23 --> 00:42:29 This would be, on the left-hand side, 598 00:42:25 --> 00:42:31 one, on the right-hand side would be C over two. 599 00:42:30 --> 00:42:36 So, I would get one equals C over two. 600 00:42:34 --> 00:42:40 Is that correct? Cosine of zero is one, 601 00:42:39 --> 00:42:45 so that's two down below. Therefore, C is equal to two, 602 00:42:46 --> 00:42:52 and that would then complete the solution. 603 00:42:51 --> 00:42:57 We would be X2 plus two over one plus cosine X. 604 00:42:57 --> 00:43:03 Now, you can do this in general, of course, 605 00:43:03 --> 00:43:09 and get a general formula. And, we will have occasion to 606 00:43:10 --> 00:43:16 use that next week. But for now, 607 00:43:13 --> 00:43:19 why don't we concentrate on the most interesting case, 608 00:43:18 --> 00:43:24 namely that of the most linear equation, with constant 609 00:43:23 --> 00:43:29 coefficient, that is, so let's look at the linear 610 00:43:27 --> 00:43:33 equation with constant coefficient, because that's the 611 00:43:32 --> 00:43:38 one that most closely models the conduction and diffusion 612 00:43:37 --> 00:43:43 equations. So, what I'm interested in, 613 00:43:41 --> 00:43:47 is since this is the, of them all, 614 00:43:43 --> 00:43:49 probably it's the most important case is the one where 615 00:43:47 --> 00:43:53 P is a constant because of its application to that. 616 00:43:50 --> 00:43:56 And, many of the other, the bank account, 617 00:43:53 --> 00:43:59 for example, all of those will use a 618 00:43:55 --> 00:44:01 constant coefficient. So, how is the thing going to 619 00:43:58 --> 00:44:04 look? Well, I will use the cooling. 620 00:44:01 --> 00:44:07 Let's use the temperature model, for example. 621 00:44:05 --> 00:44:11 The temperature model, the equation will be DTDT plus 622 00:44:09 --> 00:44:15 KT is equal to. Now, notice on the right-hand 623 00:44:13 --> 00:44:19 side, this is a common error. You don't put TE. 624 00:44:16 --> 00:44:22 You have to put KTE because that's what the equation says. 625 00:44:21 --> 00:44:27 If you think units, you won't have any trouble. 626 00:44:25 --> 00:44:31 Units have to be compatible on both sides of a differential 627 00:44:30 --> 00:44:36 equation. And therefore, 628 00:44:32 --> 00:44:38 whatever the units were for capital KT, I'd have to have the 629 00:44:36 --> 00:44:42 same units on the right-hand side, which indicates I cannot 630 00:44:40 --> 00:44:46 have KT on the left of the differential equation, 631 00:44:43 --> 00:44:49 and just T on the right, and expect the units to be 632 00:44:47 --> 00:44:53 compatible. That's not possible. 633 00:44:49 --> 00:44:55 So, that's a good way of remembering that if you're 634 00:44:52 --> 00:44:58 modeling temperature or concentration, 635 00:44:54 --> 00:45:00 you have to have the K on both sides. 636 00:44:57 --> 00:45:03 OK, let's do, now, a lot of this we are going 637 00:45:00 --> 00:45:06 to do in our head now because this is really too easy. 638 00:45:05 --> 00:45:11 What's the integrating factor? Well, the integrating factor is 639 00:45:10 --> 00:45:16 going to be the integral of K, the coefficient now is just K. 640 00:45:16 --> 00:45:22 P is a constant, K, and if I integrate KDT, 641 00:45:20 --> 00:45:26 I get KT, and I exponentiate that. 642 00:45:23 --> 00:45:29 So, the integrating factor is E to the KT. 643 00:45:28 --> 00:45:34 I multiply through both sides, multiply by E to the KT, 644 00:45:34 --> 00:45:40 and what's the resulting equation? 645 00:45:38 --> 00:45:44 Well, it's going to be , I'll write it in the compact 646 00:45:44 --> 00:45:50 form. It's going to be E to the KT 647 00:45:47 --> 00:45:53 times T, all prime. The differentiation is now, 648 00:45:53 --> 00:45:59 of course, with respect to the time. 649 00:45:57 --> 00:46:03 And, that's equal to KTE, whatever that is, 650 00:46:02 --> 00:46:08 times E to the KT. This is a function of T, 651 00:46:08 --> 00:46:14 of course, the function of little time, sorry, 652 00:46:12 --> 00:46:18 little T time. OK, and now, 653 00:46:15 --> 00:46:21 finally, we are going to integrate. 654 00:46:18 --> 00:46:24 What's the answer? Well, it is E to the, 655 00:46:22 --> 00:46:28 so, are we going to get E to the KT times T is, 656 00:46:27 --> 00:46:33 sorry, K little t, K times time times the 657 00:46:31 --> 00:46:37 temperature is equal to the integral of KTE. 658 00:46:37 --> 00:46:43 I'll put the fact that it's a function of T inside just to 659 00:46:42 --> 00:46:48 remind you, E to the KT, and now I'll put the arbitrary 660 00:46:46 --> 00:46:52 constant. Let's put in the arbitrary 661 00:46:49 --> 00:46:55 constant explicitly. So, what will T be? 662 00:46:53 --> 00:46:59 OK, T will look like this, finally. 663 00:46:56 --> 00:47:02 It will be E to the negative KT. 664 00:46:59 --> 00:47:05 That's on the outside. Then, you will integrate. 665 00:47:04 --> 00:47:10 Of course, the difficulty of doing this integral depends 666 00:47:08 --> 00:47:14 entirely upon how this external temperature varies. 667 00:47:12 --> 00:47:18 But anyways, it's going to be K times that 668 00:47:16 --> 00:47:22 function, which I haven't specified, E to the KT plus C 669 00:47:20 --> 00:47:26 times E to the negative KT. Now, some people, 670 00:47:24 --> 00:47:30 many, in fact, that almost always, 671 00:47:27 --> 00:47:33 in the engineering literature, almost never write indefinite 672 00:47:32 --> 00:47:38 integrals because an indefinite integral is indefinite. 673 00:47:38 --> 00:47:44 In other words, this covers not just one 674 00:47:40 --> 00:47:46 function, but a whole multitude of functions which differ from 675 00:47:44 --> 00:47:50 each other by an arbitrary constant. 676 00:47:46 --> 00:47:52 So, in a formula like this, there's a certain vagueness, 677 00:47:49 --> 00:47:55 and it's further compounded by the fact that I don't know 678 00:47:53 --> 00:47:59 whether the arbitrary constant is here. 679 00:47:55 --> 00:48:01 I seem to have put it explicitly on the outside the 680 00:47:58 --> 00:48:04 way you're used to doing from calculus. 681 00:48:02 --> 00:48:08 Many people, therefore, prefer, 682 00:48:04 --> 00:48:10 and I think you should learn this, to do what is done in the 683 00:48:08 --> 00:48:14 very first section of the notes called definite integral 684 00:48:13 --> 00:48:19 solutions. If there's an initial condition 685 00:48:16 --> 00:48:22 saying that the internal temperature at time zero is some 686 00:48:20 --> 00:48:26 given value, what they like to do is make this thing definite 687 00:48:25 --> 00:48:31 by integrating here from zero to T, and making this a dummy 688 00:48:30 --> 00:48:36 variable. You see, what that does is it 689 00:48:35 --> 00:48:41 gives you a particular function, whereas, I'm sorry I didn't put 690 00:48:42 --> 00:48:48 in the DT one minus two. What it does is that when time 691 00:48:48 --> 00:48:54 is zero, all this automatically disappears, and the arbitrary 692 00:48:55 --> 00:49:01 constant will then be, it's T. 693 00:49:00 --> 00:49:06 So, in other words, C times this, 694 00:49:02 --> 00:49:08 which is one, is that equal to [T?]. 695 00:49:05 --> 00:49:11 In other words, if I make this zero, 696 00:49:07 --> 00:49:13 that I can write C as equal to this arbitrary starting value. 697 00:49:12 --> 00:49:18 Now, when you do this, the essential thing, 698 00:49:15 --> 00:49:21 and we're going to come back to this next week, 699 00:49:18 --> 00:49:24 but right away, because K is positive, 700 00:49:21 --> 00:49:27 I want to emphasize that so much at the beginning of the 701 00:49:25 --> 00:49:31 period, I want to conclude by showing you what its 702 00:49:29 --> 00:49:35 significance is. This part disappears because K 703 00:49:34 --> 00:49:40 is positive. The conductivity is positive. 704 00:49:38 --> 00:49:44 This part disappears as T goes to zero. 705 00:49:41 --> 00:49:47 This goes to zero as T goes to infinity. 706 00:49:45 --> 00:49:51 So, this is a solution that remains. 707 00:49:48 --> 00:49:54 This, therefore, is called the steady state 708 00:49:52 --> 00:49:58 solution, the thing which the temperature behaves like, 709 00:49:57 --> 00:50:03 as T goes to infinity. This is called the trangent. 710 00:50:01 --> 00:50:07 because it disappears as T goes to infinity. 711 00:50:07 --> 00:50:13 It depends on the initial condition, but it disappears, 712 00:50:11 --> 00:50:17 which shows you, then, in the long run for this 713 00:50:15 --> 00:50:21 type of problem the initial condition makes no difference. 714 00:50:20 --> 00:50:26 The function behaves always the same way as T goes to infinity.