TransWikia.com

Numerical solution of differential equations, accuracy is not enough

Mathematica Asked by ClyrionChen on June 13, 2021

Clear["Global`*"];
M1=17969/10000000;
M2=9/1000;
q=2/5;
Λ=143/50000;
c=q*(M2/M1)^2/2
v=Piecewise[{{Λ^4*((1+ph[t]^2/M1^2)^(q/2)-1), ph[t]>0},
             {Λ^4*(c*Tanh[ph[t]/M2]^2), ph[t]<=0}}];
vph=D[v,ph[t]];
ro=1/2*Λ^4*ph'[t]^2+v;
p=1/2*Λ^4*ph'[t]^2-v
s1=Flatten[NDSolve[
      {Λ^4*ph''[t]+3*Λ^4*A'[t]/A[t]*ph'[t]+vph==0,
      (Λ^4*A''[t])/A[t]==-(1/6)*(ro+3p),
      ph[0]==1321/250,
      ph'[0]==-(208213/200000000000),
      A[0]==1,
      A'[0]==57/2500000},
     {ph,A},{t,0,34.42349},
     MaxSteps->[Infinity],WorkingPrecision->70]]

if I add the number of workingprecision, It will prompted to me

"Error test failure at
t==13.51983911114427669481141395951743899849861225404442210074803325679357077321916`60.;
unable to continue."
Additional supplementary:
I ultimately need an image with such a physical parameter B[t]=-H'[t]/H[t]^2, and the H[t_]:=A'[t]/A[t]/.s1. If I don’t adjust the accuracy,The Images of the B[t] is intermittent and inaccurate.If I can get as accurate as possible, The graph of B[t] will be smoothed out.So I want to improve the accuracy as much as possible, but I can only improve WorkingPrecision->52 by using the above methods.

Can you help me with the answer? Thank you very much.

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP