Advanced Search

Question: plot 2 functions versus each other

Question:plot 2 functions versus each other

golnaz 40 Maple
odeplot detools differential-equations deplot plotting
August 06 2017
0

Hello every body i need help with my program. how can i plot r&s versus each other?



 

> 

restart
> 

with(plots):
> 

with(PDEtools, casesplit, declare)

[casesplit, declare]

 (1)
> 

with(DEtools, gensys)

[gensys]

 (2)
> 

declare((Omega, H)(z));

Omega(z)*`will now be displayed as`*Omega

  

H(z)*`will now be displayed as`*H

 (3)
> 

`ε` := -1

-1

 (4)
> 

h := 70

70

 (5)
> 

r := 0.3e-3

0.3e-3

 (6)
> 

n := 42

42

 (7)
> 

``

``

> 

ode1 := -(diff(Omega(z), z))+Omega(z)*((3*(1-2*`ε`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n)))*(1+h*`ε`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`ε`*h*sqrt(r)/H(z))))-2*sqrt(Omega(z))/n)/(1+z) = 0:
> 

ode2 := -(diff(H(z), z))-(3*H(z)*(1/2))*(1-2*`ε`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n))*(1+h*`ε`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`ε`*h*sqrt(r)/H(z))))/(1+z) = 0:
> 

q := -1+(3/2*(1-2*`ε`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n)))*(1+h*`ε`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`ε`*h*sqrt(r)/H(z)))):
> 

``
> 

sys := {ode1, ode2}


ics := {H(0) = h, Omega(0) = .77}:
sol := dsolve(`union`(sys, ics), numeric, output = listprocedure, stiff = true):``

 

  
> 

 
 

 
> 

r2 := q+2*q^2+(1+z)*(diff(q, z)):

``

``

> 

s := (r2-1)/(3*(q-1/2)):
 

 

 

  



 

Download rs.mw
 

 

 



 

Download rs.mw

 
Please Wait...