bradders

32 Reputation

2 Badges

15 years, 26 days

MaplePrimes Activity


These are questions asked by bradders

how would i get maple to ingnore wht is a mathematical singularity and not a physical one, and continue picking up the solution the other side of the point??

see code below;

> eqn1 := diff(e(z)^2, z) = ((3/2)*e(z)^2-3*(v[0]*exp(-u*f(z)))^2/(2*(e(z)^2-O*(1+z)^(3*g)))+(1/2)*((g-1)*3)*O*(1+z)^(3*g))/(2*(1+z));
> eqn2 := diff(f(z), z) = sqrt(abs(1-(v[0]*exp(-u*f(z)))^2/(e(z)^2-O*(1+z)^(3*g))^2))/(e(z)*(1+z));
>
> eqn3 := subs(O = .27, g = 1, v[0] = .1, u = 0.5e-2, eqn1);
> eqn4 := subs(O = .27, g = 1, v[0] = .1, u = 0.5e-2, eqn2);
>
> sol1 := dsolve([eqn3, eqn4, f(0) = 1, e(0) = 1], numeric, maxfun = 1000000000);
> plots[odeplot](sol1, [[z, e(z), color = red]], 0 .. 2);

below is the system of ODE's im trying to solve, however i get this error which i cannot get round, can anyone help?

 

Page 1 of 1