Solution of the heat equation

Please choose:

	Case 1	Analytic solution is : exp(-a*t)*sin(x)+x*t/tend
	Case 2	Homogeneous boundary values and nonsmooth initial value u0(x) = 2x for 0<=x<=1/2 , u0(x) = 2-2x for 1/2<=x<=1)
	a self defined case	Please specify the formula for f(x,t) here, (for example x**5+t) using Fortran- conventions: This may extend over up to three lines. You may use the constants pi=π, e=exp(1) and sqrt2=sqrt(2.0) but no other variables than x and t. f(x,t) = And here we need the two boundary functions and the initial value function, also as formulas with the same limitations as above. Of course the only variable is t resp. x here! u(0,t)=φ0(t): φ0(t) = u(1,t)=φ1(t): φ1(t) = And finally u(x,0)=u0(x): u0(x) =

 

Now we need the thermal conductivity a :
a = Important: a >= 0.0000001 !

Please psecify the number of internal grid points N in x direction.
This results in delta_x=1/(N+1).
N = Important: 1<=N<=200 !

Please specify the endtime:
tend = Important: tend > 0 !

Here you must specify the tolerance for the time integrator:
tol =

Please choose the maximal order for the time integration:
1 means that you want to see Euler implicit exclusively.
ord = Important: ord in {1,...,5} !

For the graphical output please specify the number of time lines you want to see.
0 means that every internal time step is plotted.
grid = Important: 0 <= grid <= 200 !

And finally, you specify your view point via two rotation angles
Important : 0 <= xang <= 180, 0 <= zang <= 360 !
For xang = 0 = zang you get the x-axis screen horizontal, the t axis on screen vertical and the u-axis is vertically to the screen.

xang=

zang=

Click on "evaluate", in order to submit your input.