65 Reputation

4 Badges

3 years, 320 days

MaplePrimes Activity

These are replies submitted by mary120

If I have two different F(phi) functions (for example, see F1(phi) and F2(phi) in attached file):
1- How can solve these two ODE(s) and  plot phi vs x in one phi-x coordinate?

2- How can I plot phi-x in a symmetric interval of x (for example, from x=-100..100)
3- How can I export data of plot of these ODEs in the form of two distinct ascii table? (I couldn’t export data of the last figure in your last uploaded worksheet in the form of a ascii table)

I'll try to do your comment

Thanks for your answer. It was very useful.
How can I export the data of the phi-x plots in your attached file in  the form of an ASCII file or a two-column .txt or .dat file? 

Considering the problem as an ODE is a good suggestion. In fact, assuming (dphi/dx)**2+2F(phi)=0, I want to plot phi versus x.
(From your first plot, it is clear that F(phi) and dF(phi)/dx =0 at both phi=0 and nearly phi=0.09)
For an assumed F(phi), I expect the phi-x curve should be as follows(see the second Fig. in the attached file), but despite having F(phi), I could not draw the curve of phi  vs. x.

@Rouben Rostamian  
Thanks for your comment. There was a mistake in my uploaded worksheet and I correted it (Please, see the following attached file).

Tahnks a lot!

@Rouben Rostamian  
Thank you very much

@Preben Alsholm 
Thanks for your comprehensive comments.
Can I expect to come up with an answer similar to soliton waves?
How can I plot the curves for different values of delta[p] in the same frame?(for example, delta[p]=0.7,0.9)


@Preben Alsholm


You are right about giving up the IC. 

If I change the IC(s) to n(x)=0 when x=+/- infinity, what happen for curves corresponding to d[p] =0.2, 0.25? 

@Preben Alsholm 

How I can delete the shift of curves toward the RHS in such a way that the maximum of all curves locate on the y axis(n(x))?

@Preben Alsholm 
Thank you.
I try to plot RES1 for two different deta[p], namely 0.2, 0.25, by -display command, but I cauldn't. How can I do this?


I tried to plot the real and imaginary parts of a new parametric function, but I coulnd't. May you please help me.

I want to extract Re(x) and Im(x) vs y1 from the assumed eqaution and then plot Re(x) and Im(x) vs y1 for z=0.001, 0.002, 0.003 and r=0.1, 0.2, 0.3.

Thank you very much !

@Carl Love 

#   How to solve the following set of nonlinear equations and plot variation of u[i01],u[i02] and phi[d0] versus delta[d] for alpha=0.01,delta[d]=1e-4? 

Eq1:=1.359375000*10^36*delta[d]*phi[d0]-delta[e]+delta[i]+1 = 0:

Eq2:=7.217742610*10^14*u[i10]*(1-2*phi[d0]/u[i10]^2)+2.282450621*10^15*delta[i]*u[i20]*(1-2*phi[d0]/u[i20]^2)-2.206601886*10^17*sqrt(2)*delta[e]*(4*alpha*phi[d0]^2/(3*alpha+1)-16*alpha*phi[d0]/(3*alpha+1)+(24*alpha+1)*exp(phi[d0])/(3*alpha+1)) = 0:

Eq3:=5.882341298*10^8*u[i10]^2+1.203367844*10^9*delta[e]/u[i10]-8.823511948*10^7*u[i20]^2-9.633614643*10^7*delta[e]/(delta[i]*u[i20]) = 0:

Eq4:=2.406735687*10^10*delta[e]/u[i10]+9.633614644*10^8*delta[e]/u[i20]-(delta[e]*(4*alpha/(3*alpha+1)-1)+1.510416667*10^37*delta[d]*phi[d0]^2)*(5.882341298*10^8*u[i10]^2+1.203367844*10^9*delta[e]/u[i10])-1.359375000*10^36*delta[d]*phi[d0]*((1.157401021*10^(-25)*Pi)*u[i10]*sqrt(u[i10]^2+8/Pi)*(1-2*phi[d0]/(u[i10]^2+0.4244131815e-1))^1.0+1.395921345*10^(-54)*u[i10]*sqrt(u[i10]^2+8/Pi)*phi[d0]^2*ln((phi[d0]^2/(1000000000000*(u[i10]^2+0.4244131815e-1)^2)+2.005078125*10^14/(1+0.1666666667e-1*delta[e]+.5000000000*delta[i])^1.0)/(phi[d0]^2/(1000000000000*(u[i10]^2+0.4244131815e-1)^2)+(1/1000000000000)*(1-2*phi[d0]/(u[i10]^2+0.4244131815e-1))^1.0))/(u[i10]^2+0.4244131815e-1)^2+(1.157401021*10^(-25)*Pi)*delta[i]*u[i20]*sqrt(u[i20]^2+8/Pi)*(1-2*phi[d0]/(u[i20]^2+0.8488263632e-1))^1.0+1.395921345*10^(-55)*delta[i]*u[i20]*sqrt(u[i20]^2+8/Pi)*phi[d0]^2*ln((phi[d0]^2/(1000000000000*(u[i20]^2+0.8488263632e-1)^2)+2.005078125*10^14/(1+0.1666666667e-1*delta[e]+.5000000000*delta[i])^1.0)/(phi[d0]^2/(1000000000000*(u[i20]^2+0.8488263632e-1)^2)+(1/1000000000000)*(1-2*phi[d0]/(u[i20]^2+0.8488263632e-1))^1.0))/(u[i20]^2+0.8488263632e-1)^2-4.484955822*10^(-6))+6.240139645*10^(-20)*delta[d]*(7.217742610*10^14*u[i10]*(1-2*phi[d0]/u[i10]^2)+2.282450621*10^15*delta[i]*u[i20]*(1-2*phi[d0]/u[i20]^2)-2.206601886*10^17*sqrt(2)*delta[e]*(4*alpha*phi[d0]^2/(3*alpha+1)-16*alpha*phi[d0]/(3*alpha+1)+(24*alpha+1)*exp(phi[d0])/(3*alpha+1))) = 0:

Thanks for your reply. 
Please, see the attached file.

1 2 3 Page 1 of 3