Internal
problem
ID
[25467]
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
:
20
Date
solved
:
Friday, October 03, 2025 at 12:01:49 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)+2*t*diff(y(t),t)-a*t^2*y(t) = 0; dsolve(ode,y(t),type='series',t=0);
ode=t^2*D[y[t],{t,2}]+2*t*D[y[t],t]-a*t^2*y[t]==0; ic={}; AsymptoticDSolveValue[{ode,ic},y[t],{t,0,5}]
from sympy import * t = symbols("t") a = symbols("a") y = Function("y") ode = Eq(-a*t**2*y(t) + t**2*Derivative(y(t), (t, 2)) + 2*t*Derivative(y(t), t),0) ics = {} dsolve(ode,func=y(t),ics=ics,hint="2nd_power_series_regular",x0=0,n=6)