@Markiyan Hirnyk 

You'll notice that I asked the OP several followup questions before I Answered. That was so that I could find out, despite the OP's limited English, what their true issue was. That true issue has only a little to do with what you quoted. The OP is obviously happy with the Answer that I gave, so what's your point? Why are you criticizing?

Speaking of criticizing, I'd like to point out that your so-called "precedent" had nothing to do with what was asked.

@Markiyan Hirnyk All that was asked is how to produce the delay continuously, which is all that I Answered. The OP presumably knows the plot command for what they want to plot. That command wasn't provided in the Question.

@Markiyan Hirnyk 

Maple's random numbers are reproducible. For example, if you execute the code below, you should get exactly the same three random numbers as I have below.

restart:
'rand()' $ 3;

(As always, make sure to put restart in its own execution group.)

The same is true (or should be true) of all of Maple's "random" commands. Any random command that deviates from this is a bug.

@Markiyan Hirnyk 

The following very simple example generates a new plot every five seconds. Use the stop button when you want it to stop.

do
     print(plot(randpoly(x), x= -2..2));
     Threads:-Sleep(5)
end do;

Please don't post duplicate Questions. I deleted your duplication of this Question from earlier today. If you have information to add to a Question, then either edit the original Question or post a Reply to it. Either action will put the Question back at the top of the Active Conversations stack, where it'll attract more attention.

In your duplicate, you put forth the idea that the trouble was caused by the appearence of erf in the integrals. This isn't quite true. The problem is that the integral can't be computed (to the requested precision) even numerically.

So what's the point of your solution technique? to obtain an approximate solution? There's not much point in an approximate solution if it's too complicated to use to get numeric values. On the other hand, standard numerical technique for IVPs (e.g., dsolve's rkf45) works fine on this equation.

@abbeykabir 

Matrix(n,m) makes an n x m matrix of zeros; LinearAlgebra:-ZeroMatrix is superfluous.

If matrices A and B have the same number of columns, then <A,B> joins them vertically; if they have the same number of rows, then <A|B> joins them horizontally.

@tkeith022

Yes, I want to help you. So, if I understand you correctly, you want the program to update a graph ONCE every minute. Is that correct? And the precise times that the program runs are not important; it's the interval of time between runs that matters. Is that correct? 

Can you do what you want from a Maple worksheet, or from the command-line interface, which will be left running, but in an essentially idle state while it waits for the timer? Do you need the timer to be relative (e.g., "one hour from now") or absolute (e.g., "at 11:53:00 pm")? Do you need the timer to reset itself after the code has executed (e.g., "execute this code every day at 11:53:00 pm")? 

@mapleq2013 Oops, I see that you're using Maple 12. There's two examples that need slight changing. I've already changed them in the Answer, but I wanted to warn you in case you'd already tried executing them. The issue is identical in both examples. I used

[x,y]=~ [X,Y]

In Maple 12, that needs to be changed to

[x= X, y= Y]

@mapleq2013 

In worksheet "not_changing", comment out the line that gives g a value. Then follow along and see which matrices depend on g and which don't. You'll see that, in the end, B and C don't depend on g, only Bm depends on g.

Update: I see that my answer is essentially the same as Acer's simultaneous answer. However, I'll leave this answer here because my point is that you can easily trace the variable dependencies yourself by commenting out one line.

@ Digits:= 15 was just a lucky guess and my first guess; it's my favorite value for Digits. It had nothing to do with the number of equations. I wish that there were a better way to choose the Digits for fsolve. This is a very common complaint about fsolve.

@Boby 

In the code that Markiyan referred to, notice the end of the second line of code:

A:= table([0= subs(u= u[0], G)])

This sets the value of A[0].

@Allan Wittkopf 

Great answer, but your interface(prettyprint= 1) style output turns into unreadable garbage when posted to MaplePrimes. I guess that you're using the Command-Line Interface? Using prettyprint=2, Classic worksheets and Standard worksheets with either 1D or 2D input transcribe well to MaplePrimes, and they'll also provide the reader with a downloadable link to the worksheet, so that they can easily continue the work.

(I also strongly discourage the use of 2D input in MaplePrimes Answers because the reader cannot be sure of what you typed, but that's another issue.)

Your line

RES:= res[(D@@4)(f)(0)];

makes no sense to me because (D@@4)(f)(0) isn't in the set extra_bcs.

@EverteQ As long as these results for L are acceptable to you, it's fine:

term:= 2*x^3*y + 5*x^2*y^2 - 6*x*y^3:
coeffs(term, [x,y,z], 'L');

L;

I thought that you wanted the symbolic names in term, which are {x, y}.