<Text-field style="_pstyle257" layout="_pstyle257">How Does Heat Dissipate?</Text-field>

Note: You may notice differences between this Maple worksheet and the equivalent Mathematica notebook. These differences were introduced to preserve the content of these modules and were necessary because of major functional differences between Maple and Mathematica.

<Text-field style="Heading 1" layout="Heading 1">Introduction</Text-field>

OBJECTIVE: Observe the physical interpretations of the contour and level curve plots and see an application of the Fourier Series. Do you want to see physical interpretations of the contours and level curves you have been drawing? This equation of the heat will enable you to do so. The heat equation is a partial differential equation and its solution employs the Fourier series in a meaningful way. Note: Part II requires knowledge of the Fourier series.

<Text-field style="_cstyle256" layout="Heading 1"><Font size="18">Technology Guidelines</Font></Text-field>

NOTE: If you have just finished a worksheet, restart Maple before executing a new worksheet. TO OPEN SECTIONS, Click on the Arrow sign at the left hand side of the screen or select Expand All Sections from the View drop down menu.

TO STOP AN EXECUTION Click on STOP button from the toolbar.

ORDER OF EXECUTION Execute commands in the order given. Do not skip any Maple Input lines within a given worksheet

Alternatively, you can execute the entire worksheet by selecting the Execute Worksheet command from the Edit drop down menu.

SAVING WORKSHEETS.

You can save anytime to any directory you choose, and it is wise to save often. EXPERIENCING MAJOR PROBLEMS Save if appropriate, and then shut down Maple and start it up again.

<Text-field style="Heading 1" layout="Heading 1">Background</Text-field>

The following problems represent solutions to the partial differential heat equation: LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2Ji1JJW1zdWJHRiQ2JS1GLDYkUSJ3RicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiQtRiw2JFEjeHhGJ0Y3RjcvJS9zdWJzY3JpcHRzaGlmdEdRIjBGJy1JI21vR0YkNi5RIj1GJy8lJ2l0YWxpY0dRJmZhbHNlRicvRjhRJ25vcm1hbEYnLyUmZmVuY2VHRkgvJSpzZXBhcmF0b3JHRkgvJSlzdHJldGNoeUdGSC8lKnN5bW1ldHJpY0dGSC8lKGxhcmdlb3BHRkgvJS5tb3ZhYmxlbGltaXRzR0ZILyUnYWNjZW50R0ZILyUnbHNwYWNlR1EsMC4yNzc3Nzc4ZW1GJy8lJ3JzcGFjZUdGZW4tSSZtZnJhY0dGJDYoLUYjNiYtSSNtbkdGJDYlUSIxRidGRkZJLUZDNi1RMSZJbnZpc2libGVUaW1lcztGJ0ZJRktGTUZPRlFGU0ZVRlcvRlpRJjAuMGVtRicvRmduRmVvLUYyNiVGNC1GIzYkLUYsNiRRInRGJ0Y3RjdGP0Y3LUYjNiZGKy1GIzYkLUklbXN1cEdGJDYlLUYsNiRRImNGJ0Y3LUZebzYlUSIyRidGRkZJLyUxc3VwZXJzY3JpcHRzaGlmdEdGQUY3RitGNy8lLmxpbmV0aGlja25lc3NHRmBvLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmFxLyUpYmV2ZWxsZWRHRkhGN0YrRjc=, where w(x,t) represents the temperature at position x at time t in a thin rod of length L with perfectly insulated sides. The value of the positive constant LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2JC1JJW1zdXBHRiQ2JS1GLDYlUSJjRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW5HRiQ2JFEiMkYnL0Y7USdub3JtYWxGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGQUYrRkE= is determined by the material from which the rod is made.

<Text-field style="Heading 1" layout="Heading 1">Part 1: Fixed End Temperatures</Text-field>

We first examine a very simple solution to the problem where we assume that the ends of the rod are immersed in ice and are maintained at a temperature of zero degrees centigrade. If the rod is of length 4, the 0 degrees centigrade. at both ends is represented by the boundary conditions w(0,t) = w(4,t) = 0. Assume also that the initial temperature distribution within the rod is given by LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2KEYrLUYjNiYtRiw2JVEid0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Y6USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGRC8lKXN0cmV0Y2h5R0ZELyUqc3ltbWV0cmljR0ZELyUobGFyZ2VvcEdGRC8lLm1vdmFibGVsaW1pdHNHRkQvJSdhY2NlbnRHRkQvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZTLUkobWZlbmNlZEdGJDYkLUYjNiYtRiw2JVEieEYnRjZGOS1GPTYtUSIsRidGQEZCL0ZGRjhGR0ZJRktGTUZPRlEvRlVRLDAuMzMzMzMzM2VtRictSSNtbkdGJDYkUSIwRidGQEZARkBGQC1GPTYtUSI9RidGQEZCRkVGR0ZJRktGTUZPL0ZSUSwwLjI3Nzc3NzhlbUYnL0ZVRmZvLUYjNiYtRiw2JVEkc2luRicvRjdGREZARjwtRlc2JC1GIzYmRistRiM2KC1GX282JFEiMkYnRkAtRj02LVExJkludmlzaWJsZVRpbWVzO0YnRkBGQkZFRkdGSUZLRk1GT0ZRRlQtRiw2JVElJnBpO0YnRl1wRkBGZ3BGZW5GQEYrRkBGQEZARitGQEYrRkA=. We first show that LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2KEYrLUYjNiYtRiw2JFEid0YnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Y3USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGQS8lKXN0cmV0Y2h5R0ZBLyUqc3ltbWV0cmljR0ZBLyUobGFyZ2VvcEdGQS8lLm1vdmFibGVsaW1pdHNHRkEvJSdhY2NlbnRHRkEvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZQLUkobWZlbmNlZEdGJDYlLUYjNiYtRiw2JFEieEYnRjYtRjo2LlEiLEYnLyUnaXRhbGljR0ZBRj1GPy9GQ1EldHJ1ZUYnRkRGRkZIRkpGTEZOL0ZSUSwwLjMzMzMzMzNlbUYnLUYsNiRRInRGJ0Y2RjZGaG5GPUY2LUY6Ni5RIj1GJ0ZobkY9Rj9GQkZERkZGSEZKRkwvRk9RLDAuMjc3Nzc3OGVtRicvRlJGZW8tRiM2Jy1JJW1zdXBHRiQ2JS1GOjYuUS8mRXhwb25lbnRpYWxFO0YnRmhuRj1GP0ZCRkRGRkZIRkpGTEZOL0ZSUSwwLjExMTExMTFlbUYnLUYjNiYtRjo2LlEqJnVtaW51czA7RidGaG5GPUY/RkJGREZGRkhGSkZML0ZPUSwwLjIyMjIyMjJlbUYnL0ZSRmdwLUYjNiYtSSNtbkdGJDYlUSIzRidGaG5GPS1GOjYtUTEmSW52aXNpYmxlVGltZXM7RidGPUY/RkJGREZGRkhGSkZMRk5GUUZeb0Y2RitGNi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGX3EtRiM2Ji1GLDYlUSRzaW5GJ0ZobkY9RjktRlQ2JS1GIzYmRistRiM2KC1GXHE2JVEiMkYnRmhuRj1GX3EtRiw2JVElJnBpO0YnRmhuRj1GX3FGWEY2RitGNkZobkY9RjZGK0Y2RitGNkYrRjY= is a solution to the differential equation with these boundary and initial conditions. Also determine the value of the constant c for this problem.

restart; w:=(x,t)->exp(-3*t)*sin(2*Pi*x): print(cat(`The temperature at the left end of the rod is `, w(0,t))); print(cat(`The temperature at the right end of the rod is `, w(4,t))); print(cat(`The initial temperature distribution is `, w(x,0))); print(cat(`The second derivative of the temperature with respect to position is `, diff(diff(w(x,t), x), x))); wxx:= diff(diff(w(x,t), x), x): print(cat(`The first derivative of the temperature with respect to time is `, diff(w(x,t), t))); wt:=diff(w(x,t), t): c:=sqrt(wt/wxx): print(cat(`c =`, c, ` or `, evalf(c))); if wxx = 1/c^2*wt then print(cat(w(x,t), ` satisifies the differential equation.`)) else print(`The equation is not satisfies the differential equation.`); fi;

Plot the temperature function at various positions along the rod as a function of time.

plot([w(0, t), w(0.8, t), w(1.6, t), w(2.4, t), w(3.2, t), w(4, t)], t=0..1, color=[COLOR(RGB,1,0,0),COLOR(RGB,1,1,0),COLOR(RGB,0,1,0),COLOR(RGB,0,1,1),COLOR(RGB,1,0,1),COLOR(RGB,0,0,1)], labels=[t, w], title=`X-Level Plots`,axes=boxed);

Interpret these plots. What is happening to the temperature along different parts of the rod? Plot the temperature function at various times as a function of the position on the rod.

plot([w(x, 0), w(x, 0.2), w(x, 0.4), w(x, 0.6), w(x, 0.8), w(x, 1)], x=0..4, color=[COLOR(RGB,1,0,0),COLOR(RGB,1,1,0),COLOR(RGB,0,1,0),COLOR(RGB,0,1,1),COLOR(RGB,1,0,1),COLOR(RGB,0,0,1)], labels=[t, w], title=`X-Level Plots`);

Interpret the results. What is happening to the temperature function as time goes on? Plot the contours for the temperature function and the temperature function as a surface.

with(plots):

contourplot(w(x, t), x=0..4, t=0..1, labels=[x, t], contours=13, axes=boxed, filled=true, numpoints=1000);

plot3d(w(x, t), x=0..4, t=0..1, labels=[x, t, w], axes=boxed, orientation=[-50, 45]);

What is happening as LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2Ji1GLDYlUSJ0RicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2LVEnJnJhcnI7RicvRjhRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZCLyUpc3RyZXRjaHlHRkIvJSpzeW1tZXRyaWNHRkIvJShsYXJnZW9wR0ZCLyUubW92YWJsZWxpbWl0c0dGQi8lJ2FjY2VudEdGQi8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRlEtRiw2JVEoJmluZmluO0YnRjRGN0Y+RitGPg== ? Is this what you would expect? Why?

<Text-field style="Heading 1" layout="Heading 1">You Try It: Part I - Insulated Ends </Text-field>

We will next examine a very simple solution to the problem in which it is assumed that the ends of the rod are insulated so that no heat can escape. If the rod is of length 4, this is represented by the boundary conditions LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2KEYrLUYjNiYtSSVtc3ViR0YkNiUtRiw2JVEid0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYkLUYsNiVRInhGJ0Y5RjwvRj1RJ25vcm1hbEYnLyUvc3Vic2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJ0ZELyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZPLyUpc3RyZXRjaHlHRk8vJSpzeW1tZXRyaWNHRk8vJShsYXJnZW9wR0ZPLyUubW92YWJsZWxpbWl0c0dGTy8lJ2FjY2VudEdGTy8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRmhuLUkobWZlbmNlZEdGJDYkLUYjNiYtSSNtbkdGJDYkRkhGRC1GSjYtUSIsRidGREZNL0ZRRjtGUkZURlZGWEZaRmZuL0ZqblEsMC4zMzMzMzMzZW1GJy1GLDYlUSJ0RidGOUY8RkRGREZELUZKNi1RIj1GJ0ZERk1GUEZSRlRGVkZYRlovRmduUSwwLjI3Nzc3NzhlbUYnL0ZqbkZgcC1GIzYmRjNGSS1GXG82JC1GIzYmLUZhbzYkUSI0RidGREZjb0Zpb0ZERkRGREYrRkRGK0ZE = 0. Why? Can you explain this in your own words? Assume also that the initial temperature distribution within the rod is given by LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2KEYrLUYjNiYtRiw2JFEid0YnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Y3USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGQS8lKXN0cmV0Y2h5R0ZBLyUqc3ltbWV0cmljR0ZBLyUobGFyZ2VvcEdGQS8lLm1vdmFibGVsaW1pdHNHRkEvJSdhY2NlbnRHRkEvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZQLUkobWZlbmNlZEdGJDYlLUYjNiYtRiw2JFEieEYnRjYtRjo2LlEiLEYnLyUnaXRhbGljR0ZBRj1GPy9GQ1EldHJ1ZUYnRkRGRkZIRkpGTEZOL0ZSUSwwLjMzMzMzMzNlbUYnLUkjbW5HRiQ2JVEiMEYnRmhuRj1GNkZobkY9RjYtRjo2LlEiPUYnRmhuRj1GP0ZCRkRGRkZIRkpGTC9GT1EsMC4yNzc3Nzc4ZW1GJy9GUkZmby1GIzYmLUYsNiVRJGNvc0YnRmhuRj1GOS1GVDYlLUYjNiZGKy1GIzYoLUZfbzYlUSIyRidGaG5GPS1GOjYtUTEmSW52aXNpYmxlVGltZXM7RidGPUY/RkJGREZGRkhGSkZMRk5GUS1GLDYlUSUmcGk7RidGaG5GPUZmcEZYRjZGK0Y2RmhuRj1GNkYrRjZGK0Y2. First, show that LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2KEYrLUYjNiYtRiw2JFEid0YnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Y3USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGQS8lKXN0cmV0Y2h5R0ZBLyUqc3ltbWV0cmljR0ZBLyUobGFyZ2VvcEdGQS8lLm1vdmFibGVsaW1pdHNHRkEvJSdhY2NlbnRHRkEvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZQLUkobWZlbmNlZEdGJDYlLUYjNiYtRiw2JFEieEYnRjYtRjo2LlEiLEYnLyUnaXRhbGljR0ZBRj1GPy9GQ1EldHJ1ZUYnRkRGRkZIRkpGTEZOL0ZSUSwwLjMzMzMzMzNlbUYnLUYsNiRRInRGJ0Y2RjZGaG5GPUY2LUY6Ni5RIj1GJ0ZobkY9Rj9GQkZERkZGSEZKRkwvRk9RLDAuMjc3Nzc3OGVtRicvRlJGZW8tRiM2Jy1JJW1zdXBHRiQ2JS1GOjYuUS8mRXhwb25lbnRpYWxFO0YnRmhuRj1GP0ZCRkRGRkZIRkpGTEZOL0ZSUSwwLjExMTExMTFlbUYnLUYjNiYtRjo2LlEqJnVtaW51czA7RidGaG5GPUY/RkJGREZGRkhGSkZML0ZPUSwwLjIyMjIyMjJlbUYnL0ZSRmdwLUYjNiYtSSNtbkdGJDYlUSIzRidGaG5GPS1GOjYtUTEmSW52aXNpYmxlVGltZXM7RidGPUY/RkJGREZGRkhGSkZMRk5GUUZeb0Y2RitGNi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRidGX3EtRiM2Ji1GLDYlUSRjb3NGJ0ZobkY9RjktRlQ2JS1GIzYmRistRiM2KC1GXHE2JVEiMkYnRmhuRj1GX3EtRiw2JVElJnBpO0YnRmhuRj1GX3FGWEY2RitGNkZobkY9RjZGK0Y2RitGNkYrRjY= is a solution to the differential equation with these boundary and initial conditions. Also identify the constant c.

restart;

w:=(x, t)-> exp(-3*t)*cos(2*Pi*x): print(cat(`The change in temperature at the left end of the rod is `, diff(w(0, t), x))): print(cat(`The change in temperature at the right end of the rod is `, diff(w(4, t), x))): print(cat(`The initial temperature distribution is `, w(x, 0))); print(cat(`The second derivative of the temperature with respect to position is `, diff(diff(w(x, t), x),x))): wxx:=diff(diff(w(x, t), x),x): print(cat(`The first derivative of the temperarure with respect to time is `, diff(w(x, t), t))): wt:=diff(w(x, t), t): c:=sqrt(wt/wxx): print(`c= `, c, ` or `, evalf(c)); if wxx=1/c^2*wt then print(cat(w(x, t), ` satisifes the differential equation.`)); else print(cat(`The equation is not satisfied.`)); fi;

Plot the temperature function at various positions along the rod as a function of time.

Interpret these plots. What is happening to the temperature along different parts of the rod?

Plot the temperature function at various times as a function of the position on the rod.

plot([w(x, 0), w(x, 0.2), w(x, 0.4), w(x, 0.6), w(x, 0.8), w(x, 1)], x=0..4, color=[COLOR(RGB,1,0,0),COLOR(RGB,0,1,1),COLOR(RGB,0,1,0),COLOR(RGB,0,1,1),COLOR(RGB,1,0,1),COLOR(RGB,0,0,1)], labels=[t, w], title=`X-Level Plots`);

Interpret the results. What is happening to the temperature function as time goes on? Plot the contours for the temperature function and the temperature function as a surface.

with(plots): contourplot(w(x, t), x=0..4, t=0..1, labels=[x, t], contours=13, axes=boxed, filled=true,numpoints=1000); plot3d(w(x, t), x=0..4, t=0..1, labels=[x, t, w], axes=boxed, orientation=[-50, 45], labels=["x","y","z"]);

How does this differ from the temperature function in the rod with fixed ends?

<Text-field style="Heading 1" layout="Heading 1">Part II: Insulated Ends and Use of the Fourier Series (demonstration only)</Text-field>

If you have not learned the Fourier series, you should not do this part of the module. Suppose that the rod of length 4 has insulated ends, that is, LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2KEYrLUYjNiYtSSVtc3ViR0YkNiUtRiw2JVEid0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYkLUYsNiVRInhGJ0Y5RjwvRj1RJ25vcm1hbEYnLyUvc3Vic2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYtUTAmQXBwbHlGdW5jdGlvbjtGJ0ZELyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZPLyUpc3RyZXRjaHlHRk8vJSpzeW1tZXRyaWNHRk8vJShsYXJnZW9wR0ZPLyUubW92YWJsZWxpbWl0c0dGTy8lJ2FjY2VudEdGTy8lJ2xzcGFjZUdRJjAuMGVtRicvJSdyc3BhY2VHRmhuLUkobWZlbmNlZEdGJDYkLUYjNiYtSSNtbkdGJDYkRkhGRC1GSjYtUSIsRidGREZNL0ZRRjtGUkZURlZGWEZaRmZuL0ZqblEsMC4zMzMzMzMzZW1GJy1GLDYlUSJ0RidGOUY8RkRGREZELUZKNi1RIj1GJ0ZERk1GUEZSRlRGVkZYRlovRmduUSwwLjI3Nzc3NzhlbUYnL0ZqbkZgcC1GIzYmRjNGSS1GXG82JC1GIzYmLUZhbzYkUSI0RidGREZjb0Zpb0ZERkRGREYrRkRGK0ZE = 0. Assume a discontinuous initial temperature distribution that is 12 for the left LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkmbWZyYWNHRiQ2KC1GIzYkLUkjbW5HRiQ2JFEiMUYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJ0Y0LUYjNiQtRjE2JFEiNEYnRjRGNC8lLmxpbmV0aGlja25lc3NHRjMvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGQC8lKWJldmVsbGVkR1EmZmFsc2VGJ0Y0 of the rod and 0 elsewhere. This discontinuous initial temperature function necessitates the use of the Fourier series in order to get a continuous approximation to this discontinuous initial state. = 0. First, we determine the Fourier Series corresponding to the initial temperature state for the rod.

restart;

a:=n->(1/2)*int(12*cos(n*Pi*x/4), x=0..1): fourcos:=evalf((a(0)/2+sum(a(j)*cos(j*Pi*x/4), j=1..30))): print(cat(`The Fourier cosine series out to n = 30 is `, fourcos)); plot(fourcos, x=0..4);

The solution to the heat equation can now be written as LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2I1EhRictRiM2Jy1JJm1mcmFjR0YkNigtRiM2JkYrLUYjNiYtRiw2JFEiYUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RMCZBcHBseUZ1bmN0aW9uO0YnL0Y8USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGRi8lKXN0cmV0Y2h5R0ZGLyUqc3ltbWV0cmljR0ZGLyUobGFyZ2VvcEdGRi8lLm1vdmFibGVsaW1pdHNHRkYvJSdhY2NlbnRHRkYvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZVLUkobWZlbmNlZEdGJDYlLUYjNiQtSSNtbkdGJDYlUSIwRicvJSdpdGFsaWNHRkZGQkY7RltvRkJGO0YrRjstRiM2JC1GaG42JVEiMkYnRltvRkJGOy8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGZ28vJSliZXZlbGxlZEdGRi1GPzYuUSIrRidGW29GQkZERkdGSUZLRk1GT0ZRL0ZUUSwwLjIyMjIyMjJlbUYnL0ZXRmBwLUYjNigtSSttdW5kZXJvdmVyR0YkNictRj82LlEmJlN1bTtGJ0Zbb0ZCRkRGRy9GSlEldHJ1ZUYnRksvRk5GW3EvRlBGW3FGUUZTL0ZXUSwwLjE2NjY2NjdlbUYnLUYjNiYtRiw2JFEibkYnRjstRj82LlEiPUYnRltvRkJGREZHRklGS0ZNRk9GUS9GVFEsMC4yNzc3Nzc4ZW1GJy9GV0ZpcS1GaG42JUZkb0Zbb0ZCRjstRiw2JFEoJmluZmluO0YnRjtGUS8lLGFjY2VudHVuZGVyR0ZGRistSSdtc3BhY2VHRiQ2Ji8lJ2hlaWdodEdRJjAuMGV4RicvJSZ3aWR0aEdRJDUuMEYnLyUmZGVwdGhHRmdyLyUqbGluZWJyZWFrR1ElYXV0b0YnLUYjNipGKy1GIzYmRjhGPi1GWTYlLUYjNiRGYnFGO0Zbb0ZCRjstRj82LVExJkludmlzaWJsZVRpbWVzO0YnRkJGREZHRklGS0ZNRk9GUUZTRlYtSSVtc3VwR0YkNiUtRj82LlEvJkV4cG9uZW50aWFsRTtGJ0Zbb0ZCRkRGR0ZJRktGTUZPRlFGUy9GV1EsMC4xMTExMTExZW1GJy1GIzYmLUY/Ni5RKiZ1bWludXMwO0YnRltvRkJGREZHRklGS0ZNRk9GUUZfcEZhcC1GIzYmLUZcdDYlLUZZNiUtRjI2KC1GIzYmRmJxRmhzLUYsNiVRJSZwaTtGJ0Zbb0ZCRjstRiM2JC1GaG42JVEiNEYnRltvRkJGO0Zib0Zlb0Zob0Zqb0Zbb0ZCRl9vLyUxc3VwZXJzY3JpcHRzaGlmdEdGam5GaHMtRiw2JFEidEYnRjtGO0YrRjtGanVGaHMtRiM2Ji1GLDYlUSRjb3NGJ0Zbb0ZCRj4tRlk2JS1GIzYkLUYyNigtRiM2KEZicUZoc0ZidUZocy1GLDYkUSJ4RidGO0Y7RmV1RmJvRmVvRmhvRmpvRjtGW29GQkY7RitGO0YrRjtGK0Y7RitGOw== and we will write write it out to n = 30.

w1(x,t):=a(0)/2+add(a(j)*exp(-(j*Pi/4)^2*t)*cos(j*Pi*x/4),j = 1..30): evalf(%);

w:=unapply(w1(x,t),(x,t));

Plot the temperature function at various positions along the rod as a function of time.

xlevelplots:=plot([w(0, t), w(0.8, t), w(1.6, t), w(2.4, t), w(3.2, t), w(4, t)], t=0..5, color=[COLOR(RGB,1,0,0),COLOR(RGB,1,1,0),COLOR(RGB,0,1,0),COLOR(RGB,0,1,1),COLOR(RGB,1,0,1),COLOR(RGB,0,0,1)], labels=[t, w], title=`X-Level Plots`): print(xlevelplots);

Interpret these plots. What is happening to the temperature along different parts of the rod?

Plot the temperature function at various times as a function of the position on the rod. Interpret the results.

tlevelplots:=plot([w(x, 0), w(x, 0.3), w(x, .6), w(x, 1), w(x, 1.5), w(x, 2)], x=0..4, color=[COLOR(RGB,1,0,0),COLOR(RGB,1,1,0),COLOR(RGB,0,1,0),COLOR(RGB,0,1,1),COLOR(RGB,1,0,1),COLOR(RGB,0,0,1)], labels=[x, w], title=`T-Level Plots`): print(tlevelplots);

Interpret the results. What is happening to the temperature function as time goes on?

Plot the contours for the temperature function and the temperature function as a surface.

with(plots): contourplot(w(x, t), x=0..4, t=0..5, labels=[x, t], contours=13, axes=boxed, filled=true,numpoints=1000); plot3d(w(x, t), x=0..4, t=0..5, labels=[x, t, "w"], axes=boxed, orientation=[-50, 45], labels=["x","y","z"]);