salim-barzani

reply...

Answered: salim-barzani 1560

January 01 2025

0 0

@mmcdara i want do the same like paper did factoring (m+F)=(m+G'/G)

reply...

Answered: salim-barzani 1560

December 31 2024

0 0

@mmcdara

Previously, you solved a problem for me involving an equation like (m+1/G′)=ϕ. Your explanation was excellent, and I enjoyed learning the new method you used to solve it. Inspired by that, I attempted to apply the same technique to factor (m+G′G)=phi, but I couldn't figure it out this time.

I don’t fully understand the technique you demonstrated earlier, but I’d really appreciate it if you could explain it again. If possible, could you also provide a code implementation for this problem? Your previous explanation was incredibly clear and helpful. Thank you!

Thanks, you did a great job...

Answered: salim-barzani 1560

December 30 2024

0 0

@nm that is 2 day i am stuck becuase of this my second odetest make a problem i knew problem is parameter becuase in a lot place have a gamma anyway thank you so much is great it is work .

thank you so much...

Answered: salim-barzani 1560

December 30 2024

0 0

@Oliver K but union not work but try the identical is work.

don't use chatGpt...

Answered: salim-barzani 1560

December 30 2024

0 0

@janhardo i want figure out by maple.

add some more detail to question post...

Answered: salim-barzani 1560

December 29 2024

0 0

@dharr if you watch the file part in indet we have name but i remove the name instead i write down the parameter but it say is not exists how this is possible

How about upload code...

Answered: salim-barzani 1560

December 29 2024

0 0

@dharr it is easy to find all case really is a good code @ mmcdara write the code i asked letter for change to arbitrary parameter but he didn't response ,becuase i think my PDE is not going to zero becuase all parameter there must be a problem i should found i am stuck becuase my pde which i transformed to ode is not going to zero have a problem with parameter but i don't know what is it no one answered untill now i think thus parameter make a problem

A secret...

Answered: salim-barzani 1560

December 29 2024

0 0

@janhardo have a long story

if have any other question ask...

Answered: salim-barzani 1560

December 29 2024

0 0

@nm i try to explain in this file

>

(1)

>

Warning, A new binding for the name `gamma` has been created. The global instance of this name is still accessible using the :- prefix, :-`gamma`. See ?protect for details.

>

(2)

>

pde1 := I*(diff(Omega(x, t)^m, t))+alpha*(diff(Omega(x, t)^m, `$`(x, 2)))+I*beta*(diff(abs(Omega(x, t))^(2*n)*Omega(x, t)^m, x))+m*sigma*Omega(x, t)^m*(diff(W(t), t)) = I*gamma*abs(Omega(x, t))^(2*n)*(diff(Omega(x, t)^m, x))+delta*abs(Omega(x, t))^(4*n)*Omega(x, t)^m

(3)

>

ode := 4*V(xi)^4*n^2*sigma+(-4*beta*k*m*n^2+4*gamma*k*m*n^2)*V(xi)^3+(4*alpha*k^2*m^2*n^2-4*delta^2*m*n^2+4*m*n^2*w)*V(xi)^2-2*V(xi)*(diff(diff(V(xi), xi), xi))*alpha*m*n+(-alpha*m^2+2*alpha*m*n)*(diff(V(xi), xi))^2 = 0

4*V(xi)^4*n^2*sigma+(-4*beta*k*m*n^2+4*gamma*k*m*n^2)*V(xi)^3+(4*alpha*k^2*m^2*n^2-4*delta^2*m*n^2+4*m*n^2*w)*V(xi)^2-2*V(xi)*(diff(diff(V(xi), xi), xi))*alpha*m*n+(-alpha*m^2+2*alpha*m*n)*(diff(V(xi), xi))^2 = 0

(4)

>

(5)

>

(6)

>

(7)

>

(8)

>

G(xi) = -exp(-(2*m*mu+lambda)*xi)*c__1/(2*m*mu+lambda)-mu*xi/(2*m*mu+lambda)+c__2

(9)

>

diff(G(xi), xi) = -(-2*m*mu-lambda)*exp(-(2*m*mu+lambda)*xi)*c__1/(2*m*mu+lambda)-mu/(2*m*mu+lambda)

(10)

>

K := V(xi) = a[-1]/(m+1/(diff(G(xi), xi)))+a[0]+a[1]*(m+1/(diff(G(xi), xi)))

V(xi) = a[-1]/(m+1/(diff(G(xi), xi)))+a[0]+a[1]*(m+1/(diff(G(xi), xi)))

(11)

>	$case1 := {alpha = alpha, beta = gamma, delta = delta, gamma = gamma, k = k, lambda = 0, m = 2n, mu = mu, n = n, sigma = 32alphamu^2n^4/a[-1]^2, w = -2alphak^2n-4alphamu^2n+delta^2, a[-1] = a[-1], a[0] = 0, a[1] = 0}$

${alpha = alpha, beta = gamma, delta = delta, gamma = gamma, k = k, lambda = 0, m = 2*n, mu = mu, n = n, sigma = 32*alpha*mu^2*n^4/a[-1]^2, w = -2*alpha*k^2*n-4*alpha*mu^2*n+delta^2, a[-1] = a[-1], a[0] = 0, a[1] = 0}$

(12)

>

V(xi) = a[-1]/(2*n+1/(diff(G(xi), xi)))

(13)

>

128*V(xi)^4*n^6*alpha*mu^2/a[-1]^2+(16*alpha*k^2*n^4-8*delta^2*n^3+8*n^3*(-2*alpha*k^2*n-4*alpha*mu^2*n+delta^2))*V(xi)^2-4*V(xi)*(diff(diff(V(xi), xi), xi))*alpha*n^2 = 0

(14)

(15)

>

V(xi) = a[-1]/(2*n+1/(-(-2*m*mu-lambda)*exp(-(2*m*mu+lambda)*xi)*c__1/(2*m*mu+lambda)-mu/(2*m*mu+lambda)))

(16)

>

V(xi) = a[-1]/(2*n+1/(exp(-4*mu*n*xi)*c__1-(1/4)/n))

(17)

>

(18)

>

U(xi) = (a[-1]/(2*n+1/(exp(-4*mu*n*xi)*c__1-(1/4)/n)))^((1/2)/n)

(19)

>

Omega(x, t) = (a[-1]/(2*n+1/(exp(-4*mu*n*xi)*c__1-(1/4)/n)))^((1/2)/n)*exp(I*(-k*(tau*c[0]+xi)+w*tau+delta*W(tau)-delta^2*tau))

(20)

>

(21)

>

(22)

>

4*U(xi)^2*delta^2-(32*I)*U(xi)*(diff(U(xi), xi))*alpha*k-16*U(xi)^2*alpha*k^2-4*U(xi)^2*(-4*alpha*k^2-8*alpha*mu^2+delta^2)-(4*I)*U(xi)*(diff(U(xi), xi))*c[0]+(4*I)*U(xi)^5*(diff(U(xi), xi))*gamma+12*(diff(U(xi), xi))^2*alpha-512*alpha*mu^2*U(xi)^10/a[-1]^2+4*U(xi)*(diff(diff(U(xi), xi), xi))*alpha = 0

(23)

>

U(xi) = (a[-1]/(4+1/(exp(-8*mu*xi)*c__1-1/8)))^(1/4)

(24)

>

(25)

>

(26)

>

Download explain.mw

a question?...

Answered: salim-barzani 1560

December 28 2024

0 0

@mmcdara in the code i don't want all parameter maybe 4 or 5 is enough How i do that?

F-P.mw

thank you so much...

Answered: salim-barzani 1560

December 28 2024

0 0

@acer last night my work stoped becuase of this physic i will not used anymore untill i something strange appear.

nothing,...

Answered: salim-barzani 1560

December 27 2024

0 0

@acer just some time a new sign appear with physic they remove , why problem is physic?

i like your code style...

Answered: salim-barzani 1560

December 27 2024

0 0

@mmcdara is different no one can undrestand except you best of best thank you so much.

Thanks you, you did a great job...

Answered: salim-barzani 1560

December 27 2024

0 0

@mmcdara you are legend of the coding

How you did that?...

Answered: salim-barzani 1560

December 27 2024

0 0

@mmcdara i look to my old code and find another thing but is not give my answer even

How you do that is unbliavable , just this is unclear for me why the author just write a[-1],a[1],c[4]

>

(1)

>

(2)

>

(3)

>

F := sum(a[i]*(m+1/(diff(G(xi), xi)))^i, i = -1 .. 1)

>

(4)

>

NULL

>	# rewritting rule RR := isolate(m+1/(diff(G(xi), xi))=Phi, diff(G(xi), xi));

(5)

>	# Apply RR and collect wrt Phi subs(RR, L): normal(%): PhiN := collect(numer(lhs(%)), Phi): PhiD := denom(lhs(%%));

(6)

>	with(LargeExpressions): LLE := collect(PhiN, Phi, Veil[phi] ): LLE / PhiD = 0;

(3*Phi^8*phi[1]+6*Phi^7*phi[2]+Phi^6*phi[3]+Phi^5*phi[4]-Phi^4*phi[5]-Phi^3*phi[6]+Phi^2*phi[7]-6*Phi*phi[8]+3*phi[9])/Phi^4 = 0

(7)

>	# phi[i] coefficients CoefficientNullity:=print~( [ seq( phi[i] = simplify(Unveil[phi](phi[i]), size), i=1..LastUsed[phi] ) ] ):

phi[3] = -8*(mu*c[3]*(m^2*mu^2+m*mu*lambda-(7/8)*lambda^2)*c[1]^3+(3/4)*c[3]*((m^2*a[1]+a[-1])*mu^2+m*mu*lambda*a[1]-(1/2)*lambda^2*a[1])*c[1]^2+(1/2)*mu*c[1]*c[4]-(3/8)*mu*c[2]^2)*a[1]

phi[4] = -8*lambda*(c[3]*(m^2*mu^2+m*mu*lambda-(1/8)*lambda^2)*c[1]^3+(3/4)*(m*lambda*a[1]+mu*(m^2*a[1]+2*a[-1]))*c[3]*c[1]^2+(1/2)*c[1]*c[4]-(3/8)*c[2]^2)*a[1]

phi[5] = -2*(m^2*mu^2+m*mu*lambda-(1/2)*lambda^2)*((m^2*a[1]+a[-1])*mu+m*lambda*a[1])*c[3]*c[1]^3-3*((m^4*a[1]^2+4*m^2*a[-1]*a[1]+a[-1]^2)*mu^2+2*m*lambda*a[1]*(m^2*a[1]+2*a[-1])*mu+lambda^2*a[1]*(m^2*a[1]-2*a[-1]))*c[3]*c[1]^2-4*c[4]*((m^2*a[1]+a[-1])*mu+m*lambda*a[1])*c[1]+3*c[2]^2*((m^2*a[1]+a[-1])*mu+m*lambda*a[1])

phi[6] = -8*a[-1]*(c[3]*(m^2*mu^2+m*mu*lambda-(1/8)*lambda^2)*c[1]^3+(3/2)*(m*lambda*a[1]+mu*(m^2*a[1]+(1/2)*a[-1]))*c[3]*c[1]^2+(1/2)*c[1]*c[4]-(3/8)*c[2]^2)*lambda

phi[7] = -8*a[-1]*(mu^2*c[1]^2*c[3]*(mu*c[1]+(3/4)*a[1])*m^4+2*mu*lambda*c[1]^2*c[3]*(mu*c[1]+(3/4)*a[1])*m^3+((1/8)*mu*c[3]*c[1]^3*lambda^2+(3/4)*c[3]*(lambda^2*a[1]+mu^2*a[-1])*c[1]^2+(1/2)*mu*c[1]*c[4]-(3/8)*mu*c[2]^2)*m^2+(3/4)*lambda*(-(7/6)*c[1]^3*c[3]*lambda^2+mu*a[-1]*c[1]^2*c[3]+(2/3)*c[1]*c[4]-(1/2)*c[2]^2)*m-(3/8)*lambda^2*a[-1]*c[1]^2*c[3])

(8)

>	COEFFS := solve({phi[1],phi[2],phi[3],phi[4],phi[5],phi[6],phi[7],phi[8],phi[9]}, {a[-1],a[1],mu,c[4]})

(9)

>

(10)

>

sols := solve(CoefficientNullity, [a[-1], a[0], a[1], c[4]]); sols := `assuming`([eval(sols)], [b > 0]); whattype(sols); print(cat(`$`('_', 120))); `~`[print](sols)