janhardo

700 Reputation

12 Badges

11 years, 68 days

MaplePrimes Activity


These are replies submitted by janhardo

@Alfred_F 
Thanks, yes you have corrected now  your mistake what i noticed earlier, good.
"If I understand your question correctly, you are investigating a nonlinear ordinary differential equation of the structure
a*y´´+b*y^3+c*y=0."
I see there exists a compendium for non-limear odes, maybe i can track the ode there too ?

If you own this book it seems to me, then you can look for solutions to odes in a more structured way and not have to reinvent the wheel every time

Found another resource for ode solving : 12000.org/index.htm

I'm not a big fan of deleting a question.

Often repetition can make a topic clearer for the questioner and also for the person trying to help .
The ideal case is for the questioner to continue with the question they are working on and ask new questions about it , until some things become clearer to the questioner .

This does require patience on the part of the person helping and a helpful attitude

This is how @salim-barzani' s question was put away and I wanted to respond to it further and had responded to it, well..

@salim-barzani 
Copy the code and paste then in a empty worksheet and it will work when you run it
Easy going.... 

 

 

 



Is it for a 1Dimensional  wave equation or 2 Dimensional wave equation this soliton wave?
Note: i am a beginner with this 

@salim-barzani 
Probably is this the most compact expression you can get, but useless because you don't know the values for A and B, belonging to the solutions 
 

restart; with(PDEtools): with(DEtools):
sol := [seq(dsolve(eval(diff(G(xi), xi) = G(xi)^2 + A*G(xi) + B, {seq(x = 0, x in s)})), s in combinat:-powerset({A, B}))];

with a set calculation  ( set A and set B )
 

sol := [seq([s, dsolve(eval(diff(G(xi), xi) = G(xi)^2 + A*G(xi) + B, {seq(x = 0, x in s)}))], s in combinat:-powerset({A, B}))];

The power set of a set A consists of all subsets of A.
The power set of is written as P(A) 
Proposition: If set A has n elements, then P(A) has 2^n elements
P(A,B) ....

@salim-barzani Can you use deepl for translating your language into Englisch ?
Vertalen met DeepL Translate - 's werelds meest accurate vertaler

@salim-barzani 
Explain what exactly you want to compute about an pde and whether that is a repetitive computation ? 

With odetest 

restart;
with(PDEtools):
with(DEtools):

# Define a procedure to solve the ODE and check the solution
SolveODE := proc(A, B)
    local S, Solution, Check;
    
    # Set up the differential equation
    S := diff(G(xi), xi) = G(xi)^2 + A*G(xi) + B;
    
    # Solve the ODE
    Solution := dsolve(S, G(xi));
    
    # Check if the solution satisfies the ODE
    Check := odetest(Solution, S);
    
    # Return the solution and the validation check
    return Solution, Check;
end proc:

# Example usage
# Solution and test for specific values of A and B
sol1, test1 := SolveODE(A, B);  # A = A, B = B
sol2, test2 := SolveODE(0, B); # A = 0, B = B
sol3, test3 := SolveODE(A, 0); # A = A, B = 0
sol4, test4 := SolveODE(0, 0); # A = 0, B = 0
 

"maple.ini in users"

 

G(xi) = -(1/2)*A-(1/2)*tanh((1/2)*(A^2-4*B)^(1/2)*(c__1+xi))*(A^2-4*B)^(1/2), 0

 

G(xi) = tan(B^(1/2)*(c__1+xi))*B^(1/2), 0

 

G(xi) = A/(-1+exp(-A*xi)*c__1*A), 0

 

G(xi) = 1/(-xi+c__1), 0

(1)
 

 

Download mpromes_salim_procedure_dsolve_A_B_odetest_26-10-2024.mw

restart;

with(plots):

# Define a custom `rmod` function for periodic switching
rmod := (t, n) -> frac(t * n):

# Create a function for the chessboard pattern that switches periodically
f := (t, p) -> if (floor(4 * t/Pi) + floor(4 * p/Pi)) mod 2 = 0 then 1 else NULL end if:

# Plot the black chessboard squares on the sphere
p1 := plot3d(f, 0 .. 2*Pi, 0 .. Pi, coords = spherical, scaling = constrained, 
             color = "black", grid = [200, 200], style = surface, adaptmesh = false):

# Create a white sphere to make the pattern clearly visible
p2 := sphereplot(0.999, theta = 0 .. 2*Pi, phi = 0 .. Pi, 
                 color = "white", style = surface):

# Combine both plots
plots:-display(p1, p2);

If you follow these axioms in Maple , then it is proved
1 is neutral element for multiplication
0 is neutal element for addition 

Tsol := Re(value(eval(rhs(sols), infinity = 12)))


 

On school we learned also to bring terms to one side ...:-)

ode := diff(u(x), x, x)/u(x) - diff(u(x), x)^2/u(x)^2 = x^2 - diff(u(x), x)^2/u(x)^2;

# Move all terms to one side
eq := lhs(ode) - rhs(ode);

# Simplify the result
simplified_eq := simplify(eq);

        

 

@C_R 

Indeed, the cpu load was the same  before and after hibernation.
You must not have created restore points in window 10 to revert to an earlier configuration (impractible, so )?

@C_R 
Amd ryzen 7 5700x- 8-coreprocessor +SDD +16 mb 
Cpu load about average 10 % ? or lower , with or without maple in windows 11 started  , memory 40 % , network 0 % 

Seems there is very little effect from maple on the windows 11 system 
When get things cpu hungry : graphics ? 
Procesverkenner - Sysinternals | Microsoft Learn

First 24 25 26 27 28 29 30 Last Page 26 of 75