Skip to main content
edited tags
Link
bbgodfrey
  • 62.1k
  • 18
  • 92
  • 160
edited body
Source Link
Stratus
  • 3k
  • 13
  • 25

enter image description here

I'm trying to solve the above differential equation using the shooting method. When i set Method -> {"Shooting", "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1} i get a solution, but it is equal to zero. I need a solution that is different of zero. When i increase the value of z'[t]z'[0], say to z'[t]=1000z'[0]=1000 Mathematica stay running for a long time and i get no solution. How can i solve this?

Here is the code:

n = 3.;
xf = 10000.;
s = NDSolve[
   {
    -z''[t] == (1 + (n - 2) t)^(1/(2 - n) (2 (n - 1) + 0.001)) z[
        t]^2 (Sin[z[t]+2])
    , z[0] == 0, z'[xf] == 0}, z, {t, 0, xf},
   Method -> {"Shooting", 
     "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1}
     }
   ] // Chop
Plot[Evaluate[z[t] /. s], {t, 0., 10}, PlotRange -> All]

enter image description here

I'm trying to solve the above differential equation using the shooting method. When i set Method -> {"Shooting", "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1} i get a solution, but it is equal to zero. I need a solution that is different of zero. When i increase the value of z'[t], say to z'[t]=1000 Mathematica stay running for a long time and i get no solution. How can i solve this?

Here is the code:

n = 3.;
xf = 10000.;
s = NDSolve[
   {
    -z''[t] == (1 + (n - 2) t)^(1/(2 - n) (2 (n - 1) + 0.001)) z[
        t]^2 (Sin[z[t]+2])
    , z[0] == 0, z'[xf] == 0}, z, {t, 0, xf},
   Method -> {"Shooting", 
     "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1}
     }
   ] // Chop
Plot[Evaluate[z[t] /. s], {t, 0., 10}, PlotRange -> All]

enter image description here

I'm trying to solve the above differential equation using the shooting method. When i set Method -> {"Shooting", "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1} i get a solution, but it is equal to zero. I need a solution that is different of zero. When i increase the value of z'[0], say to z'[0]=1000 Mathematica stay running for a long time and i get no solution. How can i solve this?

Here is the code:

n = 3.;
xf = 10000.;
s = NDSolve[
   {
    -z''[t] == (1 + (n - 2) t)^(1/(2 - n) (2 (n - 1) + 0.001)) z[
        t]^2 (Sin[z[t]+2])
    , z[0] == 0, z'[xf] == 0}, z, {t, 0, xf},
   Method -> {"Shooting", 
     "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1}
     }
   ] // Chop
Plot[Evaluate[z[t] /. s], {t, 0., 10}, PlotRange -> All]
deleted 2 characters in body
Source Link
Stratus
  • 3k
  • 13
  • 25

enter image description here

I'm trying to solve the above differential equation using the shooting method. When i set Method -> {"Shooting", "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1} i get a solution, but it is equal to zero. I need a solution that is different of zero. When i increase the value of z'[t], say to z'[t]=1000 Mathematica stay running for a long time and i get no solution. How can i solve this?

Here is the code:

n = 3.;
xf = 10000.;
s = NDSolve[
   {
    -z''[t] == (1 + (n - 2) t)^(1/(2 - n) (2 (n - 1) + 0.001)) z[
        t]^2 (Sin[z[t]] + 2Sin[z[t]+2])
    , z[0] == 0, z'[xf] == 0}, z, {t, 0, xf},
   Method -> {"Shooting", 
     "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1}
     }
   ] // Chop
Plot[Evaluate[z[t] /. s], {t, 0., 10}, PlotRange -> All]

enter image description here

I'm trying to solve the above differential equation using the shooting method. When i set Method -> {"Shooting", "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1} i get a solution, but it is equal to zero. I need a solution that is different of zero. When i increase the value of z'[t], say to z'[t]=1000 Mathematica stay running for a long time and i get no solution. How can i solve this?

Here is the code:

n = 3.;
xf = 10000.;
s = NDSolve[
   {
    -z''[t] == (1 + (n - 2) t)^(1/(2 - n) (2 (n - 1) + 0.001)) z[
        t]^2 (Sin[z[t]] + 2)
    , z[0] == 0, z'[xf] == 0}, z, {t, 0, xf},
   Method -> {"Shooting", 
     "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1}
     }
   ] // Chop
Plot[Evaluate[z[t] /. s], {t, 0., 10}, PlotRange -> All]

enter image description here

I'm trying to solve the above differential equation using the shooting method. When i set Method -> {"Shooting", "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1} i get a solution, but it is equal to zero. I need a solution that is different of zero. When i increase the value of z'[t], say to z'[t]=1000 Mathematica stay running for a long time and i get no solution. How can i solve this?

Here is the code:

n = 3.;
xf = 10000.;
s = NDSolve[
   {
    -z''[t] == (1 + (n - 2) t)^(1/(2 - n) (2 (n - 1) + 0.001)) z[
        t]^2 (Sin[z[t]+2])
    , z[0] == 0, z'[xf] == 0}, z, {t, 0, xf},
   Method -> {"Shooting", 
     "StartingInitialConditions" -> {z[0] == 0, z'[0] == 1}
     }
   ] // Chop
Plot[Evaluate[z[t] /. s], {t, 0., 10}, PlotRange -> All]
Source Link
Stratus
  • 3k
  • 13
  • 25
Loading