While pasting the question, I made this mistake.

thx

While pasting the question, I made this mistake.

thx

thx 

thx 

@barefoot1980 

I am using Maple 16, and its working.

 solvemp.mw

@barefoot1980 

I am using Maple 16, and its working.

 solvemp.mw

@Markiyan Hirnyk , thats why I am asking about it.

@Markiyan Hirnyk , thats why I am asking about it.

There is an article about this,

Numerical Calculations Using Maple: Why & How?

Numerical_Calculatio.pdf

acer is right about its 3rd party origin.

For more info, 

http://arxiv.org/pdf/chao-dyn/9812029v1.pdf

Anyway thx.

There is an article about this,

Numerical Calculations Using Maple: Why & How?

Numerical_Calculatio.pdf

acer is right about its 3rd party origin.

For more info, 

http://arxiv.org/pdf/chao-dyn/9812029v1.pdf

Anyway thx.

As you have mentioned that the best way (with which I agree) is to have a source term

(heat generation source) in the heat equation.  For example,

(rate of energy accumulation)   =    (rate of energy in by conduction and (or) convection) 

+ (rate of heat generation)

In your case, you are only considering conduction (I guess it from the heat equation you

considered), then

(rate of energy accumulation)   =    (rate of energy in by conduction) + (rate of heat 

generation)

 which then in mathematical form will be

alpha*diff(T(x,t),t) =  diff(T(x,t),x$2) + h(x,t),

   where h(x,t) is the source term (heat energy per unit volume per unit time).

   As for as I know, you can think of

 1. h(x,t) can be assumed to be constant, say h=1,2,3,

 2. h(x,t) can assumed to be time dependent only, i.e, h(t),

 3.  h(x,t) can be assumed to be position dependent, i.e, h(x),

 4.  h(x,t) can be assumed to be both time and position dependent, i.e. h(x,t) = t*exp(x),

     or h(x,t) = t* sin(x)

Now it depends on you, which one is more suitable for your problem. 

Thanks

@Preben Alsholm 

I am getting an error that I must enter the system as a set/list.

dsolve(eval({eq2 union bc}));

thx

@Preben Alsholm 

I am getting an error that I must enter the system as a set/list.

dsolve(eval({eq2 union bc}));

thx

 @Preben Alsholm

It seems that, the author adopted the visual approach for presenting their results.

As you can see that when N is much smaller than the one they mentioned in the paper

(N=infinity=9), the results are visually more attractive.

Thanks

 @Preben Alsholm

It seems that, the author adopted the visual approach for presenting their results.

As you can see that when N is much smaller than the one they mentioned in the paper

(N=infinity=9), the results are visually more attractive.

Thanks