Internal
problem
ID
[25460]
Book
:
Ordinary
Differential
Equations.
By
William
Adkins
and
Mark
G
Davidson.
Springer.
NY.
2010.
ISBN
978-1-4614-3617-1
Section
:
Chapter
7.
Power
series
methods.
Exercises
at
page
537
Problem
number
:
10
Date
solved
:
Friday, October 03, 2025 at 12:01:44 AM
CAS
classification
:
[[_2nd_order, _with_linear_symmetries]]
Using series method with expansion around
Order:=6; ode:=t^2*diff(diff(y(t),t),t)+3*t*(1+3*t)*diff(y(t),t)+(-t^2+1)*y(t) = 0; dsolve(ode,y(t),type='series',t=0);
ode=t^2*D[y[t],{t,2}]+3*t*(1+3*t)*D[y[t],t]+(1-t^2)*y[t]==0; ic={}; AsymptoticDSolveValue[{ode,ic},y[t],{t,0,5}]
from sympy import * t = symbols("t") y = Function("y") ode = Eq(t**2*Derivative(y(t), (t, 2)) + 3*t*(3*t + 1)*Derivative(y(t), t) + (1 - t**2)*y(t),0) ics = {} dsolve(ode,func=y(t),ics=ics,hint="2nd_power_series_regular",x0=0,n=6)