## 30 Reputation

3 years, 357 days

## How to plot q vs A?...

Maple

hi all

i have tow equations and i want plot q vs.A , that A change beetwin 0..phi. vn and a are definite and positive.

plz plz help me

(diff(f,y))^2-(f[y]^4)/4+q^2/(f[y]^2)+(1+2*q^2)*(f[y]^2)/2=piecewise(y<-a,1/4+2*q^2,-a<y<a,vn*(f[y]^2)-vn+1/4+2*q^2,y>a,1/4+2*q^2);

phiB[infinity]=q*int(f,y=-infinity..infinity):
A:=phiB[infinity]:

## Problem with plot...

Maple

hi all

i have this code and want to plot(NphiB1,y). but it is not plotted.

help me plz

restart:
nB0:=1:
hb:=6.63*10^(-34):
aB:=1/(3*(3*nB0*Pi^2)^(1/3)):
gB:=4*Pi*aB/mB:
xiphase:=1/sqrt(4*mB*gB*nB0):
qcB:=8.26*10^(-3)/xiphase:
qB:=y*qcB:
v1:=-gB*nB0:
v0:=4*gB*nB0:

L:=1.5*xiphase:
k2:=sqrt(qB^2+2*mB*v1):
k3:=sqrt(qB^2-2*mB*v0):
beta:=k3*L:
W:=2*cos(beta)*(qB/k2-k2/qB):
Q:=sin(beta)*(qB/k3-k3/qB+qB*k3/k2^2-k2^2/k3/qB):
P:=sin(beta)*(qB/k3-k3/qB-qB*k3/k2^2+k2^2/k3/qB):
theta1 := -arctan((Q*(-P*Q+sqrt(Q^2*W^2+W^4-W^2*P^2))/(Q^2+W^2)+P)/W):
b1:=(theta1+k2*L)/(2*k2):
Stheta1 := -(Q*(-P*Q+sqrt(Q^2*W^2+W^4-W^2*P^2))/(Q^2+W^2)+P)/(W*sqrt(1+(Q*(-P*Q+sqrt(Q^2*W^2+W^4-W^2*P^2))/(Q^2+W^2)+P)^2/W^2)):
Ctheta1 := 1/sqrt(1+(Q*(-P*Q+sqrt(Q^2*W^2+W^4-W^2*P^2))/(Q^2+W^2)+P)^2/W^2):
phiB1=evalf(-arctan((1/4)*(cos(2*qB*b1)*(sin(beta)*(-qB/k3-k3/qB+qB*k3/k2^2+k2^2/qB/k3)+cos(theta1)*sin(beta)*(-qB/k3-k3/qB-qB*k3/k2^2-k2^2/qB/k3)+2*sin(theta1)*cos(beta)*(-qB/k2-k2/qB))+sin(2*qB*b1)*(4*cos(theta1)*cos(beta)-2*sin(theta1)*sin(beta)*(k2/k3+k3/k2)))/(-sin(2*qB*b1)*(sin(beta)*(-qB/k3-k3/qB+qB*k3/k2^2+k2^2/qB/k3)+cos(theta1)*sin(beta)*(-qB/k3-k3/qB-qB*k3/k2^2-k2^2/qB/k3)+2*sin(theta1)*cos(beta)*(-qB/k2-k2/qB))+cos(2*qB*b1)*(4*cos(theta1)*cos(beta)-2*sin(theta1)*sin(beta)*(k2/k3+k3/k2))))):
NphiB1:=evalf(phiB1/Pi);
#SNphiB1:=simplify(NphiB1);
plot(NphiB1,y);

## How can i calculate s* ? ...

Maple

hi all

i have a complex numeric like s.

how can i calculate s

## Symmetry with respect to y-axis ...

Maple

hi all

i have a plot and i want  make it symmetry into axes y.

restart;

I5:=Int((x^2)/((x^4-2*x0*x^2+x0^2+1)^(3/2)),x=0..infinity);
I6:=(1/2)*Int(1/(x^4-2*x0*x^2+x0^2+1)^(1/2),x=0..infinity);
f:=((-4/Pi)*((x0*I6-I5)/(x0*I5+I6)^(1/3)));
A:=1/(x0*I5+I6)^(2/3): #A:=Delta/eF

plot(A,x0=-3..5);

## display plots together...

Maple

hi all

i have tow plots and i want to display them together, but i cant, plz help me

restart;

I5:=Int((x^2)/((x^4-2*x0*x^2+x0^2+1)^(3/2)),x=0..infinity);
I6:=(1/2)*Int(1/(x^4-2*x0*x^2+x0^2+1)^(1/2),x=0..infinity);
f:=((-4/Pi)*((x0*I6-I5)/(x0*I5+I6)^(1/3)));
B:=x0/(x0*I5+I6)^(2/3): #B:=mu/eF
C:=B/((16/Pi^2)*((x0*I6-I5)^2/(x0*I5+I6)^(2/3))): #C:=mu/epsilon0/2

p1:=plot(B,x0=0..4):
p2:=plot(C,x0=-3..0):
display({p1,p2},x0=-3..4);

 1 2 3 4 Page 1 of 4
﻿