Paulo Baumbach

40 Reputation

0 Badges

14 years, 71 days

MaplePrimes Activity

These are replies submitted by

@Preben Alsholm 
Perfect ... Everything is working ... Thank you very much.

@Preben Alsholm 

I think the problem is still not resolved.

The boundary conditions of the ode system depend on H, which is updated within the loop (H [i + 1] = H [i] + dH).

So, with each loop we have the same ode system with different ics.

Thank you

@Preben Alsholm 

Yes, now it works very well... Perfect result, with a very short processing time...

Well, just to clarify a point... As you well observed, the loop closes at i = 742... This is natural, since the purpose of the routine is to find a value of t that fulfills the condition "if fvp> = 0 then "... After observation on i, the value of dt (specifically Ts) can be adjusted, causing the loop to be executed close to the value 20000, which improves the discretization of the problem ...

Thanks a lot for the help...

@Preben Alsholm 
As always, you are helping me a lot. Not only with the maple but also with the math. And for that, I thank you.

Well, I'm looking at your code suggestion, line by line. However, when executing it, the following error occurs in the command "resNew: = dsolve ({sys_ode} ...": "Error, (in dsolve / numeric / process_input) input system must be an ODE system, got independent variables { x, H-47.3135421221118} ".

Any tips on this?


The variable nn is equal to 20, which causes this loop to execute at this speed ... But I usually need nn = 10000 or nn = 20000 ... In these situations the time processing is much larger ...

Also, I still deal with 8 more loops of this type, with the same problem ... So processing time really is my primary issue ...

About the Bcs you mentioned, they are correct ...

Thank you


Well, in this example ( I used the optimize = true option and this worked very well...

However, the compile option did not work... I saw a discussion on the topic ( and I confess that I was a bit confused ...

Anyway, I thank you for your clarification.

@Carl Love Perfect! Thank you!


The answer is simple: I am a beginner.
Thanks for your help.
I hope your next answers are more friendly.

@Preben Alsholm 

I have a final problem. It consists of a function defined as the integral of a function that is the result of a numerical solution of a PDE.

I need to calculate values for this function (obviously) and also plot it.

Sorry for the 2D entry. Thank you!

@tomleslie I copied the command lines from your file and pasted it into a new tab. Because it does not work?

@mahmood1800 no, because in my problem g (x, t) is a numerical function. It results from a numerical solution of a PDE.

@Preben Alsholm 

Yes, I understand what you have explained and I believe that is the key point.

The procedure I adopted to solve this problem is to discretize the temporal domain and incrementally solve the system at each step of time, taking the function H (t) = cte at each step of time. To do this, I transform the PDE system into an ODE system.

Well, thank you again for your precious help.


Well, I'll try to explain.

Removing the function H(t), the PDE system can be transformed into an ODE system (it is possible to remove the dependence of the time differential dt), which has an easy resolution. Finally, the function H (t) is obtained.

Yes, this IBCs is the problem.

Thank you for trying to help me.

@Preben Alsholm 

Sorry for the 2D input. I will change my way of working.

Very interesting the way you work and the way you solved the problem.

You have helped me a lot in the last few days and for that I am grateful. Thank you.


Yes, I understood what you explained.

Well, you asked a good question. The answer reflects my limitations in using Maple. Initially, the problem has a PDE system. The material domain varies over time, that is, it is an open physical system (the physical phenomenon is similar to that found in fluid mechanics). I could not resolve this PDE system with the maple, certainly because of my limitations in using the program. Well, to solve this difficulty I approached the problem in an incremental way, decoupling the PDE system from the temporal variable. The result is an ODE system that is calculated with each step of time. That's the reason for all the loops you mentioned.

I am very grateful for your help. Thank you.


1 2 Page 1 of 2