Skip to main content
1-Visitor
August 8, 2014
Question

Prime 3.0 can't solve BVP with units?

  • August 8, 2014
  • 6 replies
  • 2713 views

Prime 3.0 can't solve BVP with units?

Without units - no problem:

BVP.PNG

With units - error:

BVP1.PNG

Viktor

6 replies

Werner_E
25-Diamond I
August 8, 2014

C:=Odesolve..., not C(x):=Odesolve...!

Unfortunately this does not solve the problem.

1-Visitor
August 8, 2014

Thank you Werner, but the problem remains

BVP_1a.PNG

Werner_E
25-Diamond I
August 8, 2014

Yes, I already relaized that. It seems to be a bug in Prime and has to do with the units in combination with that kind of boundary conditions.

If I replace one of your conditions with a more "standard" one, Prime seems to be happy. Of course thats not a solution to your problem, though.

1.png

12-Amethyst
August 9, 2014

Viktor,

This problem is discussed at this link. It's a problem caused by sbval not attaching the units to the estimated initial condition. If you run sbval separately as Josue eventually did and then reattaching the unit to the new initial condition, it works.

Also, see here for other units problems with odesolve. It explains why odesolve sometimes works with units and sometime not. When it doesn't work, try using Rkadapt, line solver routine.

Harvey

1-Visitor
August 10, 2014

Thank you Harvey.

Viktor