Hi, 

Struggling to get my head round how to plot the results of the follwoing, as it seems I have 4 variables...

I have a data file containing the results of the following problem:

1) Consider a cube. Each side (a,b,c) has a length of 1 unit.

2) Place an ion (0,0,0) where [a,b and c intersect] and find the associated energy at that point.

3) Move the ion to a point 0.05 units along side a to (0.05,0,0) and re-find the associated energy.

Reply With Quote

The goal is to get Kalman Form

https://docs.google.com/file/d/0Bxs_ao6uuBDUT25BeVhtWjZwRUU/edit?usp=sharing

https://docs.google.com/file/d/0Bxs_ao6uuBDUVXU5dXh6Rl9zUUU/edit?usp=sharing

 

I am doing the example in page 707 of Partial Differential control theory Volume 2 by J.F.

 

when using Robertz Daniel's Ore and Involutive packages

 

with(Involutive):

with(OreModules):

Alg := DefineOreAlgebra(diff=[D,t...

Find a(n)of sequence:5, 7, 7, 2, 1, 5, 6, 6, 4, ...............

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...

Hi,I want to compute some expressions as binomial(d-2,r)*binomial(v-d-1,r-1) and something like that. In the process, there are many such expressions but with differernt parameters. Some of them can pass successfully, but some of them will generate a mistake. That is: Error, (in gcd/LinZip) variable in input expression must be in either the evaluate or retained variables.Could you please what happened to it? Thanks.

Hello

 I wish to solve nonlinear PDE equation numerically  in u & v coordinate . where u=t-x, v=t+x . we note x=(v-u)/2.. 

pde := diff(u(t,x),t,t) - diff(u(t,x),x,x) + x^2*u(t,x);  

 

 We need only initial condition let f(u)=sech(u) & g(v)=sech(v). So u((i,0)=f(ih), u(0,j)=g(jh) 

we summe step size of time j =  step size of spatial i .

the stencil we use , 

 

I want to use Maple (dsolve command) to solve this problem, how do I input the command? Part (a) is enought.

I have tried the command, but Maple says "Warning, for {u(t)} arbitrary the system is inconsistent".

Thank you very much for your help.

Hello All 

I looking to write Maple code procedure to solve this equation 

restart;
with(PDEtools):
with(plots):
with(LinearAlgebra):
pde := diff(psi(t,x),t,t) - diff(psi(t,x),x,x) + V(x)*psi(t,x);

where V(x)->x^2;

I have to use this stencil 

u[i+1,j]=-u[i-1,j1]+(1-1/8*h^2*v[i-1])*u[i,j-1]+(1-1/8*h^2*v[i+1])*u[i,j+1]; where i is spatial step and j is time step

Thanks 

 

int(y^2/2 - lambda*(u*diff(y,t)-diff(u,t)-a), t);

 

with(Physics):

diff(y(t)^2/2 - lambda(t)*(u(t)*diff(y(t),t)-diff(u(t),t)-a), u(t)); # u -> something

 

diff(lambda(t),t) disappear

 

should be

diff(lambda(t),t)+Physics[`*`](lambda(t), diff(y(t), t))

 

diff(y(t)^2/2 - lambda(t)*(u(t)*diff(y(t),t)-diff(u(t),t)-a), y(t)); # y -> something

 

u*diff(lambda(t...

Hi, I want to know how to change the rectangular form to polar form, and how to manipulate the operation of them, for example the picture below. The angle of polar form contain radian and degree, how to represent them in Maple?Thank you.