45 Reputation

15 years, 79 days

you're right...

Yes, you are right. Thanks for the good remarks about the problem.

unsuccessfully...

I did not succeed . Maybe the Heaviside function can help.
Anyway, thank you very much for your help.

Another problem...

Well, after fixing some errors, I have one more problem. Now in relation to the amount of boundary conditions. I believe this is related to an earlier discussion (https://mapleprimes.com/questions/225219-Why-Does-Pdsolve-In-Numeric-Mode-Return).

Test_1.mw

Perfect...

It is true. You're absolutely right. I was careless in my placements.

Thanks a lot for the help.

Problem...

Your response exceeded my expectations. Thank you.
However, in my version 18 the shadebetween function does not seem to exist. I checked the online help and saw that this function has been available since version 15. Can this be a bug in my program?

Solved the problem...

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

almost that...

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

Perfect!...

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...

an error message in resNew: = dsolve ({s...

@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?

Yes, you're right ... But ......

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

An example......

Well, in this example (ODE.mw) 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 (https://www.mapleprimes.com/questions/205037-Compiletrue-In-Dsolvenumeric) and I confess that I was a bit confused ...

Anyway, I thank you for your clarification.

@Carl Love Perfect! Thank you!...

@Carl Love Perfect! Thank you!

@tomleslie  The answer is simple: ...

The answer is simple: I am a beginner.