restart
infolevel[dsolve] := 4;
IiIl
sys_1 := {diff(x(t), t$2)=sin(t)-x(t), x(0)=0, D(x)(0)=0}; sol_1 := dsolve(sys_1)
PCUvLUklZGlmZkclKnByb3RlY3RlZEc2JC1GJTYkLUkieEc2IjYjSSJ0R0YsRi5GLiwmLUkkc2luRzYkRiZJKF9zeXNsaWJHRixGLSIiIkYqISIiLy1GKzYjIiIhRjkvLS1JIkRHRjI2I0YrRjhGOQ==
Methods for second order ODEs: --- Trying classification methods --- trying a quadrature trying high order exact linear fully integrable trying differential order: 2; linear nonhomogeneous with symmetry [0,1] trying a double symmetry of the form [xi=0, eta=F(x)] -> Try solving first the homogeneous part of the ODE checking if the LODE has constant coefficients <- constant coefficients successful -> Determining now a particular solution to the non-homogeneous ODE building a particular solution using variation of parameters <- solving first the homogeneous part of the ODE successful
Ly1JInhHNiI2I0kidEdGJSwmLUkkc2luRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlRiYjIiIiIiIjKiYtSSRjb3NHRitGJkYvRidGLyMhIiJGMA==
sys_2 := {(M__p+M__a)*diff(x(t), t$2)=M__p*sin(t)-x(t), x(0)=0, D(x)(0)=0}; sol_2 := dsolve(sys_2)
PCUvKiYsJkklTV9fcEc2IiIiIkklTV9fYUdGJ0YoRigtSSVkaWZmRyUqcHJvdGVjdGVkRzYkLUYrNiQtSSJ4R0YnNiNJInRHRidGM0YzRigsJiomRiZGKC1JJHNpbkc2JEYsSShfc3lzbGliR0YnRjJGKEYoRjAhIiIvLUYxNiMiIiFGPi8tLUkiREdGODYjRjFGPUY+
Methods for second order ODEs: --- Trying classification methods --- trying a quadrature trying high order exact linear fully integrable trying differential order: 2; linear nonhomogeneous with symmetry [0,1] trying a double symmetry of the form [xi=0, eta=F(x)] -> Try solving first the homogeneous part of the ODE checking if the LODE has constant coefficients <- constant coefficients successful -> Determining now a particular solution to the non-homogeneous ODE building a particular solution using variation of parameters <- solving first the homogeneous part of the ODE successful
Ly1JInhHNiI2I0kidEdGJSwmKiotSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IyomLCZJJU1fX3BHRiUiIiJJJU1fX2FHRiVGMyMhIiIiIiNGJ0YzRjMsKEYyRjNGNEYzRjZGM0Y2RjJGM0YxI0YzRjdGMyooRjJGMy1GK0YmRjNGOEY2RjY=
eval(sol_2, [M__p=1, M__a=0])
Error, numeric exception: division by zero