Message Boards Message Boards

0
|
3761 Views
|
4 Replies
|
0 Total Likes
View groups...
Share
Share this post:

Error "List encountered within..." in DSolve[ ]?

Posted 2 years ago
c=3*10^8
T=0
h=1.0545718* 10^(-34)
u={1,0,0,T}/Norm[{1,0,0,T}]
v={0,1,0,T}/Norm[{0,1,0,T}]
a={0,1,0,T} ⋅{0,0,c,T}
b={0,0,1,T} ⋅{0,c,0,T}
d= {0,0,1,T} ⋅{c,0,0,T}
pde =I*h^2*(u⋅a-u⋅b)*D[psi[x, y,z,t],x] -I*h^2*(u⋅a)*D[psi[x, y,z,t],y]+
I*h^2*(u⋅b)* D[psi[x, y,z,t],z]-I*h^2*(v⋅d)*D[psi[x, y,z,t],z]+
I*h^2*(v⋅d)*(D[psi[x, y,z,t],y]) == 0

Then I try to solve it:

DSolve[pde,{psiT[x,y,z,t],psiX[x,y,z,t],psiY[x,y,z,t],psiZ[x,y,z,t],psiT[x,y,z,t],psiX[x,y,z,t],psiY[x,y,z,t],psiZ[x,y,z,t]},{t,x,y,z}]

But get the message:

"List encountered within <<1>>==0. there should be no list on either side"..

What is the cause of this, and how should I write it correctly?

Thanks

POSTED BY: Ser Man
4 Replies
Posted 2 years ago

Probably the result of preexisting values for some symbols. Try clearing all symbols in the current context and reevaluating the code.

ClearAll[Evaluate[Context[] <> "*"]]
POSTED BY: Rohit Namjoshi
Posted 2 years ago

Hi Rohit, I found some inconsistencies with cdot and asterisk symbols, after correcting these, I get a working code, but a new and quite different problem:

c=3*10^8
T=2
h=1.0545718* 10^(-34)
u={1,0,0,T}/Norm[{1,0,0,T}]
v={0,1,0,T}/Norm[{0,1,0,T}]
a={0,1,0,T}*{0,0,c,T}
b={0,0,1,T} *{0,c,0,T}
d= {0,0,1,T} *{c,0,0,T}
pde =I*h^2*(u*a-u*b)*D[psi[x, y,z,t],x] -
I*h^2*(u*a)*D[psi[x, y,z,t],y]+I*h^2*(u*b)* D[psi[x, y,z,t],z]-
I*h^2*(v*d)*D[psi[x, y,z,t],z]+I*h^2*(v*d)*(D[psi[x, y,z,t],y]) == 0

Then when solving it:

dsol=DSolve[pde,{psiT[x,y,z,t],psiX[x,y,z,t],psiY[x,y,z,t],psiZ[x,y,z,t],psiT[x,y,z,t],psiX[x,y,z,t],psiY[x,y,z,t],psiZ[x,y,z,t]},{t,x,y,z}];

Supplied equations are not differential or integral equations of the given functions.

This happens no matter what values the vectors attain.

This should not be the case, if the vectors don't cancel out.

POSTED BY: Ser Man
Posted 2 years ago

Evaluate the variable pde. You will see the source of the problem.

enter image description here

POSTED BY: Hans Milton
Posted 2 years ago

As written in the discussion, this problem is there what ever the coefficients are. In your case they give 0=0. If you change them, they don't give this equality, and still it doesn't solve. Try this:

c=3*10^8
T=2
h=1.0545718* 10^(-34)
u={1,0,0,T}/Norm[{1,0,0,T}]
v={0,1,0,T}/Norm[{0,1,0,T}]
a={0,1,0,T}*{0,0,c,T}
b={1,1,1,T} *{c,c,c,T}
d= {1,0,1,T} *{c,c,c,T}

    pde =I*h^2*(u*a-u*b)*D[psi[x, y,z,t],x] -
I*h^2*(u*a)*D[psi[x, y,z,t],y]+
I*h^2*(u*b)* D[psi[x, y,z,t],z]-I*h^2*(v*d)*D[psi[x, y,z,t],z]+
I*h^2*(v*d)*(D[psi[x, y,z,t],y]) == 0

    dsol=DSolve[pde,{psiT[x,y,z,t],psiX[x,y,z,t],psiY[x,y,z,t],psiZ[x,y,z,t],
psiT[x,y,z,t],psiX[x,y,z,t],psiY[x,y,z,t],psiZ[x,y,z,t]},{t,x,y,z}];

DSolve", "deqx", "\"Supplied equations are not differential or integral equations of the given functions.\

So this 0=0 is not the problem. Something else is.

POSTED BY: Ser Man
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract