I am running Maple 6 on a Mac G3 and I am very pleased with its

performance. However, I do have a few questions about dsolve and someone

out there will probably have answers.

Here is the command I enter

(ic is a set of initial conditions.)Quote:> ans:=dsolve({eq1,eq2,eq3}union ic,[f(y),A(y),w(y)],type=numeric);

To plot my solutions I enter, for example,

1) That is very awkward; is there a better way?Quote:>plot('rhs(ans(y)[5])',y=0..8.7,title=`w vs y`);

2) If my next command is

the output isQuote:> ans(2);

[y = 2, f(y) = 4.27870735371794453, diff(f(y),y) = 1.66784699355773735,

A(y) = .350729197486994881, w(y) = .161303888719663890]

However, the system is singular and so, for large y (with infolevel set

to 2 ) I get

ans: Last values returned:Quote:>ans(20);

ans: t = 5.06476724690444068

ans: h = .126926496554092524e-13

ans: y =

ans: 11.2284126795603996

ans: 1616788.64420042886

ans: .371735536516273383e-6

ans: .99999891610531232

Error, (in ans) rkf45 is unable to achieve requested accuracy

I want to get hold of the last value of t returned; how do I do that? (Of

course, I can read it off but I want to assign that number to a variable.)

3) Does dsolve have the remember option? I want to enter, for example,

Does Maple have to recalculate the values?Quote:> plot('rhs(ans(y)[4])*rhs(ans(y)[5])',y=0.001..8.5,title=`A w vs y`);

4) I usually use aliases to simplify typing; for example, I would like to

say alias(f=f(y)). I cannot make dsolve accept aliases when I have

initial conditions. The condition f(0)=1 is misinterpreted as f(y)(0)=1

which makes no sense. Asny suggestions?

Thanks.