I am using Odesolve to solve a system of ODE's with a boxcar forcing function that starts at 1 goes to 10 and and I would like to return to 1. If I set the return value at 1.5, everything is fine. I I reduce it to, say, 1.2 I get the (to me, cryptic) error "The return value of this function must match the problem size." (see attached worksheet). Can someone help me out.
Right-click on Odesolve and choose, say Adams/BDF (or anything other than Radau) and it works ok with llim = 1
(also, perhaps vhat(t):=1+9*(phi(t-1)-phi(t-5)) might be a neater way of defining vhat - though this doesn't help Radau to work).
Well, actually you have a bigger problem there: you plot for v(t) starts in 0, but in the solve block it is stated that v(0)=1. Same thing for U.
> you plot for v(t) starts in 0, but in the solve block it is stated that v(0)=1.
This only looks that way in the plot because the ordinate is just showing beginning with 1. Seems to be all OK here:
> Same thing for U.
Why do you think so? The plot for U is not shown, but it starts at (0/0) as demanded:
You're right, U is OK, I confuse U with q4. But I see some problems in v even can't see where, but first place to check is the set of equations and initials conditions.
Thanks. (I think you gave me similar advice a while ago on a similar issue. My excuse: I use MathCad sporadically and my memory is not what it used to me)