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.

Hi guys,

I would like to solve a system of equations and inequalities with symbolic parameters (D, E, M, N, O, P, Q, R, S, T, U, V, W, X, Y, Z) for several variables (a, b, c, d, e, f, i, j, l, m, n, p, q, r, s, t, u, v).

Unfortunately the command below does not work. I get the message "Warning, solutions may have been lost".

"Solve({ l-u=0, m-v=0, n-u=0, p-v=0, Y+q-u=0, Y+r-v=0, Z+s-u=0, Z+t-v=0, (a-P*M)*l=0, (b-Q*M...

Hello, everyone. I have been messed up by a system of 50 nonlinear parametric equations with 50 unknowns(see the code in the file named nonlinear ). Since I am new to Maple, I don't know whether it is possible to make such a system work out by software. I've tried the SOLVE and FSOLVE command but get no result. Are there any other ways or toolbox for such kind of systems? Is there any typing wrongs in my code? I am sincerely looking forward to  receive any ideas...