Hi, everyone!

I need help.

There are a system of 2 pde's:

diff(Y(x, t), x$2) = exp(-2*x*b)*(A(x, t)-Y(x, t)), diff(A(x, t), t) = exp(-2*x*b)*(Y(x, t)-A(x, t))

and initial and boundary conditions:

A(x, 0) = 0, Y(0, t) = 0.1, (D[1](Y))(0, t) = 0.

Goal:

For each b = 0, 0.05, 0.1.

1)to plot 3-d Y(x,t): 0<=x<=20,0<=t<=7.

2)to plot Y(x,4).

Are there any methods with no finite-difference mesh?

I realized the methods such as pds1 := pdsolve(sys, ibc, numeric, time = t, range = 0 .. 7) can't help me:

*Error, (in pdsolve/numeric/match_PDEs_BCs) cannot handle systems with multiple PDE describing the time dependence of the same dependent variable, or having no time dependence *

I found something, that can solve my system analytically:

pds := pdsolve(sys), where sys - my system without initial and boundary conditions. At the end of the output: huge monster, consisted of symbols and numbers :) And I couldn't affiliate init-bound conditions to it.

I use Maple 13.