Skip to main content
added 58 characters in body
Source Link
Stratus
  • 3k
  • 13
  • 25

I need help to find a solution to the differential equation bellow. The solution obtained with this code is zeroa very small negative number in all the domain, and i think it should bewas mathematically proved that exists a polinomialpositive solution to this problem. Any hints?

solu = NDSolve[{-z''[t] == (Sin[2 + z[t]] z[t]^2)/(1 + t)^4.5,z[0] == 0., z'[10000] == 0.}, z, {t, 0., 10000}, Method -> {"Shooting","StartingInitialConditions" -> {z[0] == 0., z'[0] == 0.15}}]
LogLinearPlot[Chop[Evaluate[z[t] /. solu]], {t, 0.1, 10000}]

I need help to find a solution to the differential equation bellow. The solution obtained with this code is zero in all the domain, and i think it should be a polinomial solution. Any hints?

solu = NDSolve[{-z''[t] == (Sin[2 + z[t]] z[t]^2)/(1 + t)^4.5,z[0] == 0., z'[10000] == 0.}, z, {t, 0., 10000}, Method -> {"Shooting","StartingInitialConditions" -> {z[0] == 0., z'[0] == 0.15}}]
LogLinearPlot[Chop[Evaluate[z[t] /. solu]], {t, 0.1, 10000}]

I need help to find a solution to the differential equation bellow. The solution obtained with this code is a very small negative number in all the domain, and it was mathematically proved that exists a positive solution to this problem. Any hints?

solu = NDSolve[{-z''[t] == (Sin[2 + z[t]] z[t]^2)/(1 + t)^4.5,z[0] == 0., z'[10000] == 0.}, z, {t, 0., 10000}, Method -> {"Shooting","StartingInitialConditions" -> {z[0] == 0., z'[0] == 0.15}}]
LogLinearPlot[Chop[Evaluate[z[t] /. solu]], {t, 0.1, 10000}]
Source Link
Stratus
  • 3k
  • 13
  • 25

Need Help to solve a ordinary differential equation

I need help to find a solution to the differential equation bellow. The solution obtained with this code is zero in all the domain, and i think it should be a polinomial solution. Any hints?

solu = NDSolve[{-z''[t] == (Sin[2 + z[t]] z[t]^2)/(1 + t)^4.5,z[0] == 0., z'[10000] == 0.}, z, {t, 0., 10000}, Method -> {"Shooting","StartingInitialConditions" -> {z[0] == 0., z'[0] == 0.15}}]
LogLinearPlot[Chop[Evaluate[z[t] /. solu]], {t, 0.1, 10000}]