Carl Love

Carl Love

28035 Reputation

25 Badges

12 years, 319 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

@Stretto If is your series, do

eval(S, sin= 0);

To understand how and why this syntax works: Any number, such as the 0 in this case, can be used as if it were a constant function.

Please supply your procedure in plaintext form. And it's about time that you learned to format and indent them properly. So instead of using some automatic plaintext conversion, please retype your procedure as indented plaintext. Why should I need to do that?

@Kitonum Yes, I was just trying to emphasize that the default definition of series convergence seems to be different between Maple and Mathematica, and that Maple's is the standard, Calculus II, definition.

I have no problem with there being an alternate definition named formal, or even multiple alternate definitions, as long as it's clear which definition is being used. 

You can't make a variable independent just by declaring it to be so. The equation that you show can be solved for r in terms of phi and Z. Doesn't that violate r's supposed independence?

Why can't you adapt the Answer given for your previous extremely similar Question to this new function? What trouble are you having with that? Show an attempt.

And why can't you stop posting poorly formatted output, like I asked before?

@Magma Yeah, of course, I've already read it a few times. Still thinking about it. 

@acer The fact that formal works as you show suggests that Mathematica's default definition of series convergence is nonstandard. Perhaps we should call it A New Kind of Convergence.

@AzU Thank you for reposting the question. Here is my dsolve(..., numeric) solution. I've tried several other solution methods that confirm that the plot is correct. At the bottom of my worksheet, I show how to get a Taylor-polynomial solution. But getting a polynomial solution accurate over the interval 0..1000 will take a huge number of terms.
 

restart:

params:= [               #  Units
                         #------------
   C[th]    = 18,        #  J/K
   R[th]    = 16,        #  K/W
   A        = 1.5e-3,    #  m^2
   sigma    = 5.6704e-8, #  W/m^2/K^4
   T[a](t)  = 295,       #  K
   P[el](t) = 20         #  W
]:

dl1:= C[th]*diff(T[ovn](t),t) =
    P[el](t) - (T[ovn](t) - T[a](t))/R[th] - A*sigma*(T[ovn](t)^4 - T[a](t)^4);

C[th]*(diff(T[ovn](t), t)) = P[el](t)-(T[ovn](t)-T[a](t))/R[th]-A*sigma*(T[ovn](t)^4-T[a](t)^4)

IC:= T[ovn](0) = T[a](t):

NumSol:= dsolve(eval({dl1, IC}, params), numeric):

plots:-odeplot(NumSol, t= 0..1000);

Order:= 10:

SymbSol:= dsolve(eval({dl1, IC}, params), T[ovn](t), series)

T[ovn](t) = series(295+(499999999981/450000000000)*t-(1113036845457704599871/506250000000000000000000)*t^2+(6441121405411812445085050147777/3417187500000000000000000000000000000)*t^3-(2347655292260113721049166962857979230657/3075468750000000000000000000000000000000000000000)*t^4+(131208795361063611624804313442381111345483974469581/34599023437500000000000000000000000000000000000000000000000000)*t^5-(302325084826550881101749707338386065872088785176098082624471/29192926025390625000000000000000000000000000000000000000000000000000000000)*t^6+(929309305224972476960223479032630077755331210552390088591299354693259/459788584899902343750000000000000000000000000000000000000000000000000000000000000000000)*t^7+(140165775258414248849855327966762249447795082286687479526593600681255477993893591/4138097264099121093750000000000000000000000000000000000000000000000000000000000000000000000000000000)*t^8-(204522916792385017204666839604013922862934219651427786304817666760359425602767258073396631/3351858783920288085937500000000000000000000000000000000000000000000000000000000000000000000000000000000000000000)*t^9+O(t^10),t,10)

 


 

Download IVP.mw

 

@Christian Wolinski You are conflating two distinct definitions of constant (noun):

  1. A symbol which for the purposes of a specific problem or solution process is not to be considered a variable or an independent variable.
  2. A symbol representing a specific number or other fixed, well-defined mathematical entity.

When you speak of redeclaring constants in the context of a solver, you are using the first definiton. However, the predefined (yet modifiable) global constants refers to the second definition. There's no need to mix the two. Maple lets you create your own environment variables to your heart's content. So, if you want to use _EnvConstants for the solver purpose that you describe (which I agree is potentially useful), you may do so. 

@ondrejkubu The only part of that process that doesn't seem totally automatic to me is choosing which independent variable to declare in the identity clause. Let me know if you have any trouble with that. 

@Christian Wolinski I do think that there should be a built-in type var equivalent to And(name, Not(constant)). I get tired of typing that. 

Thank you for the additional information about your problem. Now I understand it, and I'm about to post a new Answer below.

@AzU  You've done a despicable thing by deleting your question after I helped you. It's like you spit in my face. I suspect that not many people here would be willing to help you if they suspected that you were going to delete the question. 

@ActiveUser Again, the example that you show is not a permutation. It doesn't have any representation in disjoint cycle form. When a permutation is applied to a list, the result is a list with the same elements as the original, each occuring the same number of times as in the original; only their order may be different.

Neither list that you show is itself a permutation, nor are those lists permutations of each other.

First 240 241 242 243 244 245 246 Last Page 242 of 708