65 Reputation

One Badge

1 years, 209 days

MaplePrimes Activity

These are replies submitted by Zeineb


I agree with you, there is something not clear which option 1 or 2 is correct and why when we change the boundary condition one get the solution 


The system of equation is correct. It comes from Finite difference method and we add some boundary condition.

I we change sin(Pi*i*h) by zero we get the solution.

Please see modified code


With zero as boundary condition i obtain a solution 

what do you think


Thnak you for your idea, as fsolve can not give the solution, maybe one can use iterative method to get an approximate solution


@Preben Alsholm 

Its a good answer, I verified directly i find the same solution.

Using maple can we simplify diff(g(x),x) to get an answer using piecewise function maybe we can use convert diff(g(x),x) using piecewise function

Many thanks

@John May 

I agree with you, I can make test only with small vale of n =6 or 7 

Maybe there are others methods to solve this system...



@Markiyan Hirnyk 

Many thank for your remarks but i tried the code , I get the following error :

Error, `DirectSearch` does not evaluate to a module



Yes that I mean. But maximization and not minimization

I trided your code using maximization but there is display 

@Rouben Rostamian  


If we work using discrete dynamical system

we have S_t(x)=x^(0.5^t);

using this solution depend on t how can we say if a set if local attractor or not


I tried your code, i get the following error

Error, `DETools` does not evaluate to a module


DETools:-DEplot*(diff(x(t), t)-x(t), x(t), t = -1 .. 1, {x(0) = 1});



Many thanks


Thank you

I agree with you.

But, I which case that the ODE may not have a solution.

Many thank



@Mariusz Iwaniuk 

Thank you for your idea in this subject.

transonic gaz equation

I tried to solve the equation in a square with homogeneous boundary conditions

Can we put non-homogeneous boundary condition for the same PDEs, like


bc[2] := u(0, y) = a, u(1, y) = b, u(x, 0) = c, u(x, 1) = d;

where a,b,c and d are given real number

or maybe we can put a general function as boundary condition like

bc[2] := u(0, y) = f1(y), u(1, y) = f2(y), u(x, 0) = g1(x), u(x, 1) = g2(x)

where the four functions f1 g2 g1 ang g2 are given

Many thanks




@Mariusz Iwaniuk 

Thank you for the previous ideas presented to solve this problem.


Its not less than one but you say  its less than one

In the vertical axis there are many point that the finction is not less than one

but in the horizontal axis the moduls of the function is less than one


 how can i make a loop on row or colunm to construct Anew matrix from A matrix

1 2 Page 1 of 2