G=6.674*^-11 m^3/kg/sek^2;
c=3*^8m/sek;
M0=2*^30kg;
M1=36 M0;
M2=29M0;
rs=2 G (M1+M2)/c^2;
r0=5.689rs;
μ=M1 M2/(M1+M2);
L0=μ r0^2 ω0;
ω0=Sqrt[G (M1+M2)/r0^3];
vi=ω0 r0;
τ=5c^5 r0^4/256/G^3/(M1 M2)/(M1+M2);
z=0.093;
d=1*^25m;
   
v[r_,t_]:=Sqrt[G (M1+M2)/r[т]^3]r[т]Sqrt[1-rs/r[t]] ;
f[r_,t_]:=v[r,t]/r[t]/π;
P[t_]:=F Evaluate[(ε'[t]Sqrt[1-rs/r[t]])/.solNum][[1]];
λ[t_,r_]:=Evaluate[-(4 P[t] G)/(f[r,t] d c^4)/.solNum][[1]];
 
kg=1;m=1;sek=1;

(*Rechnung mit evaluiertem r,farbige Kurve*)

solNum=Quiet[NDSolve[{r'[t]==-64/5 G^3/c^5(M1 M2) (M1+M2)/r[t]^3,r[0]==r0,L'[t]==-32/5 G^(7/2)/c^5 M1^2 M2^2 Sqrt[M1+M2]/r[t]^(7/2),L[0]==L0,ε'[t]==-32/5 G^4/c^5 M1^2 M2^2 (M1+M2)/r[t]^5,ε[0]==0},{r,L,ε},{t,0,τ}]];

(* Fitting Kurve *)

т=Quiet[t/.FindRoot[Evaluate[r[t]/.solNum][[1]]/rs-2, {t,τ/2}]];
T=If[t<т,t,2т-t]; F=If[t<т,1,1+т-t]^500;

(*Proberechnung mit gefittetem und interpoliertem r,gestrichelte Kurve*)

ř=Interpolation[{{0.25,180/38rs},{0.275,174/38 rs},{0.3,168/38 rs},{0.325,160/38 rs},{0.35,150/38 rs},{0.375,138/38 rs},{0.4,112/38 rs},{0.425,55/38 rs}},InterpolationOrder->1];

solFit=Quiet[NDSolve[{L'[t]==-32/5 G^(7/2)/c^5 M1^2 M2^2 Sqrt[M1+M2]/ř[t]^(7/2),L[0]==L0,ε'[t]==-32/5 G^4/c^5 M1^2 M2^2 (M1+M2)/ř[t]^5,ε[0]==0},{L,ε},{t,0,τ}]];

(*Plot*)

Quiet[Grid[{
{
Plot[{(L0-Evaluate[L[t]/.solNum][[1]])/L0},
{t,0.25,0.424},Frame->True,PlotRange->All,PlotStyle->{Magenta},
ImageSize->260,ImagePadding->{{30,8},{14,8}},PlotStyle->"Classic", 
GridLines->{{0.3,0.35,0.4},{}}],

Plot[{-(Evaluate[ε[t]/.solNum][[1]])/M0/c^2},
{t,0.25,0.424},Frame->True,PlotRange->All,PlotStyle->{Cyan},
ImageSize->260,ImagePadding->{{30,8},{14,8}},PlotStyle->"Classic",
GridLines->{{0.3,0.35,0.4},{}}]
},{
Plot[{Evaluate[r[t]/.solNum][[1]]/rs},
{t,0.25,0.424},Frame->True,PlotRange->All,PlotStyle->{Orange},
ImageSize->260,ImagePadding->{{30,8},{14,8}},PlotStyle->"Classic", 
GridLines->{{0.3,0.35,0.4},{}}],

Plot[{Evaluate[f[r,t]/.solNum][[1]]},
{t,0.25,0.424},Frame->True,PlotRange->All,PlotStyle->{Blue},
ImageSize->260,ImagePadding->{{30,8},{14,8}},PlotStyle->"Classic",
GridLines->{{0.3,0.35,0.4},{}}]
},{
Plot[{-Cos[Pi t Evaluate[f[r,t]/.solNum][[1]]](λ[T,r]1*^21)^(1/3)},
{t,0.25,0.424},Frame->True,PlotRange->All,PlotStyle->{Red},
ImageSize->260,ImagePadding->{{30,8},{14,8}},PlotStyle->"Classic",
GridLines->{{0.3,0.35,0.4},{}}],

Plot[{-Cos[Pi t Evaluate[f[r,t]/.solNum][[1]]]λ[T,r]1*^21},
{t,0.25,0.424},Frame->True,PlotRange->All,PlotStyle->{Darker[Red]},
ImageSize->260,ImagePadding->{{30,8},{14,8}},PlotStyle->"Classic",
GridLines->{{0.3,0.35,0.4},{}}]
},{
Sound[Play[-Cos[Pi t Evaluate[f[r,t]/.solNum][[1]]]λ[t,r]1*^21,{t,0,0.424},SampleRate->44100]]
}}, Alignment->Center]] (*yukterez.net*)

({{, }, {, }, {, }})