C_R

1960 Reputation

19 Badges

5 years, 317 days

MaplePrimes Activity


These are replies submitted by C_R

An example of "the" equations would be helpful to better understand what you are trying to achieve

there are no determinants in your post

@Ronan 

The start page is another hard to reproduce trouble maker. Have you tried to disable it.

@acer Thanks for the 'Simple way' to check

What you are describing sounds very familiar to me (and I think I have even posted something here, that I cannot find straight away). It is most likely not a bug but related to Windows (7 and 10) and your personal setup.

A few questions:

What happens if you wait (for me Maple sometimes started after a few minutes)?

Any recent changes on your computer or network configuration?

Does Maple start always from the start menu (i.e. without double click on a file)

It is super annoying and not reproducible. For me it disappeared as it appeared several times over the past 5 years. I remember that this effect was more pronounced after a fresh system start-up. Also changing/ connecting/disconnecting networks seemed toplay a role.

PS.:My biggest suspicion at the moment is Java in combination with Windows.

@ That's what I get

An example would be helpfull.

Maple uses the method exact to solve the ode. Slighty better results can be obtained this way

factor(ode);
%*denom(lhs(%));
dsolve(%)

 

 

This way the first 5 solutions can be identified as the roots of the first term encricled in red. The sixth solution (in yellow) is now what we would expect since a different method (separable) was applied

 

Your original sixth solution is implicit which can be made explicit this way

simplify(isolate(diff(sol[6],x),diff(y(x),x)));dsolve(%)

I can't tell why Maple returns an implicit solution

@nm 

Interesting design decision that seems to save memory space, as no power statements are required in arithmetic expressions.

dismantle(a^b);

POWER(3)
   NAME(4): a
   NAME(4): b

 

@jganding 

This makes it a bit more generic

expr :=1* Unit('h');
evalf(convert(expr, 'unit_free'))*indets(expr, specfunc(Units:-Unit))[];

 

Do you have an example?

@SwissMapleBeginner 

You are welcome. Assuming that dharrs solution is what you want, I have added an extact and generic solution, that you can apply to any event in a piecewise fashion. Event by event you can not only redefine the ICs but also the (temperature dependend?) parameters of the ODEs. (Assembling from piecewise solutions a piecewise function for plotting is the better way for plotting as compared to plots,display.)

restart;

ode_sys:=diff(SA2(t),t)=SA1(t)*k12-SA2(t)*k21,
         diff(SA1(t),t)=-k1x*SA1(t);
ics:=SA2(t_ev)=0,SA1(t_ev)=Dose;

diff(SA2(t), t) = SA1(t)*k12-SA2(t)*k21, diff(SA1(t), t) = -k1x*SA1(t)

 

SA2(t_ev) = 0, SA1(t_ev) = Dose

(1)

sols:=dsolve({ode_sys,ics})

{SA1(t) = Dose*exp(-k1x*t)/exp(-k1x*t_ev), SA2(t) = -(-k1x*Dose*k12/(exp(-k21*t_ev)*(k1x-k21))+k21*Dose*k12/(exp(-k21*t_ev)*(k1x-k21)))*exp(-k21*t)/(k1x-k21)-Dose*k12*exp(-k1x*t)/((k1x-k21)*exp(-k1x*t_ev))}

(2)

params:=k12=0.0452821,k21=0.0641682,k1x=0.00426118;

k12 = 0.452821e-1, k21 = 0.641682e-1, k1x = 0.426118e-2

(3)

p:=subs(t_ev=1,Dose=484,params,sols)

{SA1(t) = 484*exp(-0.426118e-2*t)/exp(-0.426118e-2), SA2(t) = -365.8425407*exp(-0.641682e-1*t)/exp(-0.641682e-1)+365.8425406*exp(-0.426118e-2*t)/exp(-0.426118e-2)}

(4)

plot(rhs~(p),t=1..50);

 

 


 

Download Split_from_Event.mw

An excellent initiative. It's a shotgun approach, but why not try it instead of developing an algorithm that tries to classify before selecting a method.
That will be complicated. As an alternative to classification and to increase speed, perhaps machine learning could be applied to a large dataset to create a method selector.

Test case where combine,trig is key

expr2 := (a^8*(r^10 - 10*m*r^9 + 5*(a^2 + 8*m^2)*r^8 - 80*m^3*r^7 + 10*(a^4 + 12*a^2*m^2 + 8*m^4)*r^6 + 4*(-15*a^4*m - 40*a^2*m^3 - 8*m^5)*r^5 + 10*(a^6 + 12*a^4*m^2 + 8*a^2*m^4)*r^4 + 40*(-a^6*m - 2*a^4*m^3)*r^3 + 5*(a^8 + 8*a^6*m^2)*r^2 - 10*a^8*m*r + a^10)*cos(theta)^8 - 28*(-(30*sin(theta)^2*a^4*m*r^5)/7 + r^10 - 10*m*r^9 + 5*(a^2 + 8*m^2)*r^8 + 40*(-a^2*m - 2*m^3)*r^7 + 10*(a^4 + 12*a^2*m^2 + 8*m^4)*r^6 + 32*(-5*a^2*m^3 - m^5)*r^5 + 10*(a^6 + 12*a^4*m^2 + 8*a^2*m^4)*r^4 + 40*(-a^6*m - 2*a^4*m^3)*r^3 + 5*(a^8 + 8*a^6*m^2)*r^2 - 10*a^8*m*r + a^10)*a^6*r^2*cos(theta)^6/3 + 14*(-(260*sin(theta)^2*a^6*m*r^3)/7 + r^10 - 10*m*r^9 + 5*(a^2 + 8*m^2)*r^8 + 40*(-a^2*m - 2*m^3)*r^7 + 10*(a^4 + 12*a^2*m^2 + 8*m^4)*r^6 + 4*(-15*a^4*m - 40*a^2*m^3 - 8*m^5)*r^5 + 10*(a^6 + 12*a^4*m^2 + 8*a^2*m^4)*r^4 - 80*a^4*m^3*r^3 + 5*(a^8 + 8*a^6*m^2)*r^2 - 10*a^8*m*r + a^10)*a^4*r^4*cos(theta)^4 - 4*a^2*(-10*a^8*m*r*sin(theta)^2 + r^10 - 10*m*r^9 + 5*(a^2 + 8*m^2)*r^8 + 40*(-a^2*m - 2*m^3)*r^7 + 10*(a^4 + 12*a^2*m^2 + 8*m^4)*r^6 + 4*(-15*a^4*m - 40*a^2*m^3 - 8*m^5)*r^5 + 10*(a^6 + 12*a^4*m^2 + 8*a^2*m^4)*r^4 + 40*(-a^6*m - 2*a^4*m^3)*r^3 + 5*(a^8 + 8*a^6*m^2)*r^2 + a^10)*r^6*cos(theta)^2 + r^8*(a^2 - 2*m*r + r^2)^5/9)*m^3*r/((r^2 + cos(theta)^2*a^2)^9*(a^2 - 2*m*r + r^2)^5)

 

 

I just tried to upload and insert contents to your earlier question (I used reply).

This worked.

I am now deleting the upload

@one man 

with Maple 2024 I get:

Error, invalid input: too many and/or wrong type of arguments passed to allvalues; first unused argument is Pi*(2*_Z1+1)

With this change

restart;
f := sin(x)/x:
allvalues({solve(f, x, allsolutions)}[1 .. 2])[];
                          

                    Pi (2 _Z1 + 1), 2 Pi _Z1

the output is still the same.

For _Z1=0 I still get Pi and 0. 0 is what I want to have removed form the solutions.

1 2 3 4 5 6 7 Last Page 1 of 36