Skip to main content
14-Alexandrite
February 8, 2026
Solved

Odesolve Block won't solve

  • February 8, 2026
  • 1 reply
  • 259 views

Hi, I have been trying to get this non-linear gear dynamics system of ODE's to solve but for the life of me cannot seem to correct the errors.

The attached sheet has several attempts at the end to write the equations and initial conditions in different formats but all give "too few initial conditions" or a"floating point error". Like most programming errors its probably glaringly simple!

Any assistance would be appreciated.

Regards,

Ross

 

Best answer by Werner_E

The error tracing feature is misleading and incorrect in this case.

 

Your amplitude vector A starts with a value zero!

Using this the function Tzz will throw a division by zero error and so does Δkf used in the solve block.

 

So I guess you have three options:

1) Omit the first entry in vector A and use i := 1 .. N3 instead of i := 0 .. N3

2) define vector A so that it does not start with value AA := 0  but rather AA := 1 or AA := 10-6 or whatever seems appropriate

3) Rewrite function  Tzz so that it will return a meaningful result if A=0.

 

Being lazy by nature I used option #1

grafik.png

and got this result

grafik.png

1 reply

25-Diamond I
February 8, 2026
It seems you haven't attached a file!
 
Due to a bug in the forum software, which has persisted already for an unbearably long time and still hasn't been fixed, directly attaching Mathcad and Prime files fails (depending on the browser and the attachment method).
 
The workaround is to put the file into a ZIP archive and then attach that here.
Yes, that's quite annoying, and not without a certain irony that the Prime forum doesn't accept Prime files, but accepts other file types like ZIP or PDF.
 
remslie14-AlexandriteAuthor
14-Alexandrite
February 8, 2026

Werner,

Thanks for that lightening fast response. I have attached a Zip file. Lets see how that goes.

Regards,

Ross

25-Diamond I
February 8, 2026

I must confess that your sheet is far too demanding for me so I could not go into the details.

What i noticed us that your odesolve-block function xv fails for any value in vector A and the error traces back to this expression

Werner_E_0-1770536917260.png

I am not sure if the problem might have to do with the fact the the magnitude of Deltak(0) is in the range of as high as 10^10 ?