SamuelTuvare

5 Reputation

2 Badges

12 years, 287 days

MaplePrimes Activity


These are replies submitted by SamuelTuvare

 

it works!

 

it works!

 

it works!

 

it works!

Hello,

 I have a small pde system 'sys' and Maple solves it correctly 'sol:=pdsolve(sys)'.  When I test 'pdetest(sol,sys)' I get '[0,0]' what is good. 

 But when I test with the same solution, only directly given by me, not found by Maple the expressions during the test are not simplyfied.

 

with(PDEtools,casesplit,declare);

with(DEtools,gensys);

declare(n2(y2,y3),n3(y2,y3),lambda,nu2,nu3,mu);

eq1:=((2*mu+lambda)*diff(n2(y2,y3),y2)+lambda*(diff(n3(y2,y3),y3)))=-lambda;

eq2:=((2*mu+lambda)*diff(n3(y2,y3),y3)+lambda*(diff(n2(y2,y3),y2)))=-lambda;

sys:=[eq1,eq2];

sol:=pdsolve(sys);

pdetest(sol,sys); #=[0,0]

# but if I take solMine, a solution entered by me, of the form  of sol and test I don't get explicitly zero! but a #non-simplifyied term!

solMine:={n2 = -(1/2)*lambda*y2/(mu+lambda), n3 = -(1/2)*lambda*y3/(mu+lambda)};

pdetest(solMine,sys);

# What shall I do to test my own candidte solutions?

 

Thanks,

Daniel

 

Page 1 of 1