Is it possible to solve (numerically or symbolically) the system of PDEs
sys:={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)) }
under the conditions
ibc:={A(x, 0) = 0, Y(0, t) = 0.1, D[1](Y)(0, t) = 0},
 where the parameter b takes the values 0,0.05,0.1, in Maple? The ranges are t=0..7, x=0..20.

my codes list below.note the last line of the codes.evalf is not effective? the output is a long expression. how to force the maple to evaluate the long expression into a numeric value?

--------------------codes---------------------

restart:
 with(Tolerances):
 ALL := [`$`(1 .. 3)]:
 solution_k := simplify(solve(map(proc (i) options operator, arrow:
 1/k[i] = (cos(Phi)*(x[i]-Tx)+sin(Phi)*(y[i]-Ty))/(-sin(Phi)*(x[i]-Tx)+cos(Phi)*(y[i]-Ty)) end proc, ALL), [Tx, Ty, Phi]), size):
 assign(solution_k):
 k := map(proc (i) options operator, arrow:
 tan(phi[i]) end proc, ALL):
 deg2rad := Tolerances:-`*`(Pi, Tolerances:-`^`(180, Tolerances:-`-`(1))):
 phi0[1] := Tolerances:-`*`(Pi, Tolerances:-`^`(3, Tolerances:-`-`(1))):
 phi0[2] := Pi:
 phi0[3] := Tolerances:-`-`(Tolerances:-`*`(Pi, Tolerances:-`^`(3, Tolerances:-`-`(1)))):
 alpha := Tolerances:-`*`(Pi, Tolerances:-`^`(10, Tolerances:-`-`(1))):
 L := 1:
 phi := phi0:
 x := map(proc (i) options operator, arrow:
 cos(phi0[i]+alpha)*L end proc, ALL):
 y := map(proc (i) options operator, arrow:
 sin(phi0[i]+alpha)*L end proc, ALL):
 angleError := Tolerances:-`&+-`(0, Tolerances:-`^`(10, -4)):
 locError := Tolerances:-`&+-`(0, Tolerances:-`^`(10, -4)):
 phi := map(proc (i) options operator, arrow:
 phi[i]+angleError end proc, ALL):
 x := map(proc (i) options operator, arrow:
 x[i]+locError end proc, ALL):
 y := map(proc (i) options operator, arrow:
 y[i]+locError end proc, ALL):
 evalf(Phi)

-----------------------------outputs--------------------

1. How can I represent the following in Maple?

(1) B(t+1) = B(t)+X (t,t+1)–D(t,t+1)

(2) T(t+1) =(m +1)/m(X(t,t+1)+X(t–1,t)+…+X (t–m+1,t–m+2))

restart:assume(M>0);

Eq1 := diff(psi(y), y, y, y, y)-M*(diff(psi(y), y, y))-Gr*b*y = 0;

bcs1:=psi(0)=0,(D@@2)(psi)(0)=0,psi(h)=-F/2,D(psi)(h)=A;

res1:=(dsolve(Eq1));

res2:=(dsolve({Eq1,bcs1},psi(y)));

match(rhs(res2)=rhs(res1),y,s);

s:

C3:=eval(_C3, s);

I am unable to find the constants. Anyway around this?

i am solving 4 ODE with boundary condition..

> restart;
> with*plots;

then i got this error..

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

i dont know where i need to change.. could you help me..

Hello, I'm new here, I hope there is someone willing to give me a hand with my problem. I've been trying to solve a problem with a for loop and I have errors in my code, could you please tell me what's wrong with the  bold part?

restart;

with(plots);

g := proc (F, A, B, X0, g::evaln) local f, a, b, E, x0, erreur, N, n, x1, ps, pp, pic; f := F; a := A; b := B; x0 := X0; E :=

evalf(abs(f(x0)));

erreur := 0.1e-3;

N := 100;

n := 0;

for i to N do x1 := evalf(x0-f(x0)/(D(f))(x0));

E := abs(f(x1));

n := n+1; ps[n] := plot([[x0, 0.], [x0, f(x0)], [x1, 0.]], color = green, thickness = 2);

x0 := x1 end do; g := evalf(x1); pp := plot(f(x), x = a .. b, color = red); pic := seq(ps[i], i = 1 .. n); display(pp, pic) end proc;
Warning, `i` is implicitly declared local to procedure `g`
f := proc (x) options operator, arrow; sin(x)-cos(2*x)+.2 end proc;
x -> sin(x) - cos(2 x) + 0.2
a := 0; b := 6;

Thank you very much!

What is the weak solution integral equation for 

d² u(x)/dx² -(1+x² ) u(x)-1=0

This is better explained with simple example. Suppose I have these 2 equations

a:='a';b:='b';x:='x';
eqs:=[3*a+4*b = 3*x, 5*a+7*b = 7*x];

What I want is to solve for (a/x) and (b/x) from the above, and not for "a" and "b".

Currently, I do this manually as follows. First divide both equations by 'x' to obtain the ratio. I had to do this manually as writing eqs/x gave expression which was hard to work with for what I want to do. I could not simplify it for what I wanted. So I typed this

eqs2:=[3*(a/x)+4*(b/x)=3,5*(a/x)+7*(b/x)=7];   %manually typed !

Now used algsubs to replace (a/x) by new variable z1, and replace (b/x) by new variable z2, then solve for z1,z2 (btw, I had to do algsubs one at a time. (need to look at this later, but not impotant)

algsubs(a/x=z1,eqs2);
algsubs(b/x=z2,%);
solve(%,{z1,z2});

       {z1=-7,z2=6}

Is there a way or function(s) to automate this. The hard part is telling Maple to use the ratio as new variable starting from set of equations.

Hy i need your help experts,

I have à small program for some variables calculations and im intrested to plot loop result,

Here is the program:

I just want to assign the results of my loop to a matrix to be easily ploted.

Thak you for help.

Dear all;

Please how can I plot the error between the two function.

Hi:

can I use Laplace Transform Method for solve Nonlinear ODE in maple?

This has been buggin me all day, but if i define a variable with an expresion say "a:=3 x+5" and i then want to assign this expression to a function like so "f:=x->a". if i then call f(7), maple return with "3 x+5" so just returning the value of a without substituting x with 7.

Is there any way in which i can define a function this way?

Hello everybody,

Since I am geting older, I need to fix the font size form 12 to 14.  I have been tryong to find where to do that permanent change.  But without success.  Could someone help to to resolve this problem.  Thank you in advance for your help.

--------------------------------------
Mario Lemelin

Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

I have a function. It's calculating coefficients of a polynomial.

So when it's done, I have the variables a0, a1, a2,a3 ... aN

where N is an input value and could be any positive integer.

Now, I just want to show them on the screen. So if I have, say , a0=1.3 and a1=5.76, I want this to appear:

1.3x + 5.76x²

That's all. I just want a loop that does this:

for term from 1 to N do:

    print( a0 ); print (x^N); print ("+");

od:

so all the terms come out on the same line, as in my example above. I've tried so many ways, but each time it ignores the numerical values of a0 and just writes "a0" on the screen, or breaks it up into a set, or puts || marks everywhere, or does it on sequential lines, or something or something or something.

I have Maple 16 for Windows 64-bit. How to get a printable Maple 16 Quick Reference Card? Also, how to save the Maple 16 Quick Reference Card as a separate file? (When I try to copy the Card I get en error message "selection too large").