Question:
For the Laplace equation
The separation of variables
solution for this problem is
Fourier coefficientsmust chosen so that they satisfy
Plot this separation of variables solution for 0 and for a few values greater than zero like 1, 0.5, 0.25. Then comment on whether a solution exists at 0 and for 0.
This example should illustrate that typical improperly posed problems might have solutions if the data are perfectly smooth and their Taylor series have finite radii of convergence. But if there is a singularity, like the kink in the triangular profile, all bets are off.
You might know that if you talk about instability of ordinary differential equations, you wonder about what happens to the solution for infinite time. But in this problem you do not let the time” coordinate go to infinity. The problem is not large , but large “wave number
. The large wave number problem is really unique to partial differential equations. (If you had a system of infinitely many ordinary differential equations, you might also run into it.)
Include your code, if any.
Answer:
Skim through the Fourier series tables of a table book for an appropriate trangular wave . You will find that the Fourier coefficients are:
Given that, you can plot the solution
Of course, on a computer you must stop summing at some finite value of . Try plotting increasing values of to figure out what the infinite sum will look like.
First plot 0. You should find that the triangular profile is produced just fine by the infinite sum.
Now plot a nonzero value of , like 0.5. Try different values 3, then 5, 10, 20, 40, and 100. Keep your eyes on the numbers on the axes. Use the results to argue that the exact solution at 0.5 does not exist.
When using matlab or octave, you could use a program like badlap_run.m. You should first look at what is in there. Then you need to save it in your matlab working directory as a file called badlap_run.m
. You will also need the function that performs the sum, badlap_sum.m. You need to save this in your matlab working directory as a file called badlap_sum.m
.
Then inside matlab or octave issue the command badlap_run
to plot. Change the values of n_max
and y
as needed and run the program again. Use help print
for information on how to make hardcopies of some of the plots you need to make your case.