Hi,

While using pdsolve for a coupled system of pdes, the maple results doen't match with the desired ones.

Here is the system

restart:

Eq1:=diff(f(eta,tau),eta$3)+(f(eta,tau)+h(eta,tau))*diff(f(eta,tau),eta$2)-diff(f(eta,tau),eta$1)^2

+(1+epsilon*cos(Pi*tau))*theta(eta,tau)=Omega*diff(f(eta,tau),tau);

Eq2:=diff(h(eta,tau),eta$3)+(f(eta,tau)+h(eta,tau))*diff(h(eta,tau),eta$2)-diff(h(eta,tau),eta$1)^2

+c*(1+epsilon*cos(Pi*tau...

Hello everyone,

I have a system of PDEs

restart:with(PDEtools):with(plots,implicitplot):

Pr:=1:A:=1:lambda:=1:Omega:=1:epsilon:=1:C1:=1:N:=5:b:=1:fw:=1:a:=1:

Eq1:=diff(f(eta,tau),eta,eta,eta)-diff(theta(eta,tau),eta)*diff(f(eta,tau),eta,eta)+

f(eta,tau)*diff(f(eta,tau),eta,eta)-(diff(f(eta,tau),eta))^2-(1+A*(1-theta(eta,tau)))*diff(f(eta,tau),eta)

+lambda*(1+epsilon*cos(Pi*tau))*theta(eta,tau)-Omega*diff(diff(f(eta,tau),tau),eta);

Hello everyone,

I am dealing with an Eigen value problem, the equations are

restart:with(plots):

Eq1:=diff(f(y),y$2)-a^2*f(y)+a*(h(y)+R*q(y))=0;

Eq2:=diff(h(y),y$2)-a^2*h(y)+a*Z*y*f(y)=0;

Eq3:=diff(q(y),y$2)-a^2*q(y)+a*f(y)=0;

ic:=f(0)=0,f(1)=0,D(h)(0)=0,q(0)=0,h(1)=0,q(1)=0;

where f,h,q are Eigen functions, R, Z are dimensionless...

Hello everyone,

I saw a question on a different forum related to the stability of the numerical solution of a coupled ode system. The original question was about the convergence of the far field condition (infinity) for some specific values of the different parameters involved.  For the details please see the link below

http://mathematica...

Hello dear,

I want to plot a piecewise function.

The code is 

restart:with(plots):

lambda0:=PL-R3*x;

Psi0:=PL+R3*S0*(R1-1)-R1*R3*x;

S0:=(Q-1)/Q;

P:=piecewise(x<S0,lambda0,x>S0,Psi0);

R3=-5,k=4,R1=0.0006,Q=1.2

I want to plot P vs x (0..1) vs PL (0..1).

Thanks