differentiaalvgl oplossen met machtreeksmethode
differentiaalvgl oplossen met machtreeksmethode
Hallo,
Ten eerste lukt het me niet te copy pasten van de equation editor naar het document dus ik moet het even zonder doen:
Ik probeer de volgende dv op te lossen met de machtreeksmethode:
y' = 1 + y^2
Nu is het zo dat je dan een machtreeks van y en y' krijgt van de vorm y = som van 0 naar inf van Am*(x)^m en y' = som van 1 naar inf van Am*m*(x)^m-1. Ik krijg het alleen niet voor elkaar nu om de som in de kwadratische y term in de dv zodanig in te vullen dat ik er wat nuttigs uit zie komen (Moet ik de haakjes gewoon uitwerken?). Met scheiding van variabelen is de dv prima op te lossen met als antwoord tan(x+c) maar ik kom er met de MRM niet uit. Ik herken ook niet echt een taylor reeks van de tangens in mijn eindantwoord.
Alvast bedankt voor hulp.
Ten eerste lukt het me niet te copy pasten van de equation editor naar het document dus ik moet het even zonder doen:
Ik probeer de volgende dv op te lossen met de machtreeksmethode:
y' = 1 + y^2
Nu is het zo dat je dan een machtreeks van y en y' krijgt van de vorm y = som van 0 naar inf van Am*(x)^m en y' = som van 1 naar inf van Am*m*(x)^m-1. Ik krijg het alleen niet voor elkaar nu om de som in de kwadratische y term in de dv zodanig in te vullen dat ik er wat nuttigs uit zie komen (Moet ik de haakjes gewoon uitwerken?). Met scheiding van variabelen is de dv prima op te lossen met als antwoord tan(x+c) maar ik kom er met de MRM niet uit. Ik herken ook niet echt een taylor reeks van de tangens in mijn eindantwoord.
Alvast bedankt voor hulp.
Re: differentiaalvgl oplossen met machtreeksmethode
Dit probleem is niet geschikt om met machtreeksen op te lossen.
Re: differentiaalvgl oplossen met machtreeksmethode
Hmm dat is gek want ik werk voor mijn natuurkundestudie uit het boek van kreyszig advanced engineering mathematics en daarin wordt dit expliciet gevraagd. Heb je toch misschien nog een idee?
Re: differentiaalvgl oplossen met machtreeksmethode
Deel door 1+y^2, en integreer vervolgens. (Komt de 1+y^2 je niet bekend voor?)
Given that, by scientifical reasons, the state of an object is completely determined by the physical influence of its environment, the probability to roll six with a dice is either one or zero.
Re: differentiaalvgl oplossen met machtreeksmethode
Zie mijn eerste post barto. Nogmaals met scheiding van variabelen is deze dv prima met substitutie op te lossen maar het betreft nu de vraag een oplossing te vinden met de machtreeksmethode. Iemand tips?
Re: differentiaalvgl oplossen met machtreeksmethode
Dat gaat niet.
Neem de oplossing .
De functie bestaat niet eens in 0.
Neem de oplossing .
De functie bestaat niet eens in 0.
Re: differentiaalvgl oplossen met machtreeksmethode
Hier een stukje uit het boek en ze staat er inderdaad tussen. Oefening 12.
Het is wel een feit dat de oefening zich niet echt leent tot een oplossing met machtreeksen, maar toch...
dus
...
Er zal wel uitgezocht moeten worden wat het convergentiegebied is...
Het is wel een feit dat de oefening zich niet echt leent tot een oplossing met machtreeksen, maar toch...
dus
...
Er zal wel uitgezocht moeten worden wat het convergentiegebied is...
Laatst gewijzigd door wnvl op 15 sep 2013, 23:47, 1 keer totaal gewijzigd.
Re: differentiaalvgl oplossen met machtreeksmethode
Matematica levert deze oplossing op
initconds = {y[0] == 0};
odeOperator = D[#, x] -1 - #^2 &;
yy = Series[y[x], {x, 0, 15}];
soln = SolveAlways[Join[{odeOperator[yy] == 0}, initconds], x];
truncatedSol = Normal[yy /. soln]
initconds = {y[0] == 0};
odeOperator = D[#, x] -1 - #^2 &;
yy = Series[y[x], {x, 0, 15}];
soln = SolveAlways[Join[{odeOperator[yy] == 0}, initconds], x];
truncatedSol = Normal[yy /. soln]
Laatst gewijzigd door wnvl op 16 sep 2013, 00:29, 1 keer totaal gewijzigd.
Re: differentiaalvgl oplossen met machtreeksmethode
In het blauw de oplossing met de machtreeksen. In het rood een andere numerieke methode.
ode = y'[x] - y[x]^2 - 1 == 0;
approxSol = NDSolve[Join[{ode}, initconds], y[x], {x, 0, 1}];
Plot[{y[x] /. approxSol[[1]]}, {x, 0, 1}]
Zou de convergentie ophouden bij ???
ode = y'[x] - y[x]^2 - 1 == 0;
approxSol = NDSolve[Join[{ode}, initconds], y[x], {x, 0, 1}];
Plot[{y[x] /. approxSol[[1]]}, {x, 0, 1}]
Zou de convergentie ophouden bij ???
Re: differentiaalvgl oplossen met machtreeksmethode
In het boek staat een héél ander probleem. Hier is gegeven dat y(0)=0.
In dat geval is het mogelijk; in het algemene geval volstrekt onmogelijk.
In dat geval is het mogelijk; in het algemene geval volstrekt onmogelijk.
Re: differentiaalvgl oplossen met machtreeksmethode
Het is natuurlijk andersom, blauw is een tangens is de echte oplossing.wnvl schreef:In het blauw de oplossing met de machtreeksen. In het rood een andere numerieke methode.
ode = y'[x] - y[x]^2 - 1 == 0;
approxSol = NDSolve[Join[{ode}, initconds], y[x], {x, 0, 1}];
Plot[{y[x] /. approxSol[[1]]}, {x, 0, 1}]
Zou de convergentie ophouden bij ???