Hello my dear friends:

i will obtain value f(4) from below equation:(f is funciton of t)

f[0](t)=a*cos(omega*t)
f[n+1](t)=f[0](t)+(1/omega)*int((omega^2*f[n](tau)-f[n](tau)^3)*sin(omega*(t-tau)),tau=0..t)

THANKS...

Hi:

how obtain value a4 from below equation in maple:

a[0]=1

a[n+2]=a[0]+a[n+1]

I calculate eigenvalues and  eigenvectors of a floating-point square matrix M with the command Eigenvectors(M, output = 'list').  How can I estimate errors of my results?

I have two 6x1 Matrices which are the results of a calculation process in Maple. One with a set of equations and the other one with a set of variables: 

A := [0, f(x6), f(x6), 0, 0, f(x6)];

b := [x1, x2, x3, x4, x5, x6];

I'd like to solve the following system:

for i from 1 to 6 do

eq[i] := A[i] = b[i]:

od;

which is

eq[1] := 0 = x1;

eq[2] := f(x6) = x2;

eq[3] := f(x6) = x3;

...

If I type in the eqations manually, and execute "s := solve({eq[1],..,eq[6]},{x1,..,x6})" everything solves fine.

If I use the "for i from..." - structure, and execute "s := solve({eq[1],..,eq[6]},{x1,..,x6})" I get an empty space as solution.

I've tried to convert both matrices into lists, but it doesn't work.

Could it be that Maple doesnt know that x6 has to be the x6 in the function f(x6) ?

Can anyone tell me how to solve this please?

Hi,

I need your help to fix the error in this code.  many thinks
restart:
pde:=diff(u(t,x),t$2)=diff(u(t,x),x$2)-sin(u(t,x)):
f:=x->x^2:
IBC := {u(0,x)=f(x),u(t,-50)=0,D[2](u)(t, 50)=0,D[1](u)(0, x)=-diff(f(x),x)}:

pds := pdsolve(pde, IBC, numeric):

p1 := pds:-plot(t = 0);
p2 := pds:-plot(t = 1/10);
p3 := pds:-plot(t = 1/2);
p4 := pds:-plot(t = 1);
p5 := pds:-plot(t = 2);
plots[display]({p1, p2, p3, p4, p5}, title = `Sine Gordaon at t=0,0.1,0.5,1,2`);

Hi Maple-Prime-ers!

I have a system of equations, containing 18 variables and 13 equations, making this a 5 degree of freedom (DOF) system.  I would like to analytically solve each of the equations in terms of each of these DOFs.  Normally I would use solve(system, dof_variables) to accomplish this, but it doesn't return anything.  Not even [].

I can solve this system by hand.  I've included a hand-solution involving isolate() and subs() in the attached worksheet.  I'm looking to incorporate this in an optimization algorithm with varying system, so I would like an automated way of doing this.

Does anybody have any suggestions to get solve to work as intended?

3driversys_FD_BRAKE_ICE_GEN.mw

Here is the system I am talking about:

The free variables are:  {FD_T, FD_W, ICE_T, EM2_T, BRAKE_T}

I'm looking for a solution in this form:

For context, I'm designing a work sheet based around quantum tunneling. Currently I'm looking at the boundary conditions.

What I want to be able to do is to set the expression psi[1] equal to psi[2], but only for the value x = 0. Is this possible? I've tried using if statements, and I've considered converting these expressions into functions for this purpose, but I'm not having much luck. 

Thanks

Blanky

Dear all,

Please, I need help to remplace in the code attached,  addadd(beta[i]*u(x+[i]*h),i=0..N)+Error; 
by add(beta[i]*u(x+alpha[i]*h),i=0..N)+Error;  I see, that we must add a function, alpha, but how.

I think I must modify this rows: but how:

stencil := add(a[ii[i]]*subsop(r=op(r,phi)+ii[i]*stepsize,phi),i=1..N);

Question2.mw

Many thinks

Hi,

I need to solve the linear system  AX=F  without using any  Linearalgebra package. 

A is a tridiagonal matrix ( band matrix);

X is vector of size n.

F is the right hand side of my equation AX=F.

F: given, A given, Find X???

restart;

A := n-> Matrix(n,scan=band[1,1],[[seq(b[i],i=2..n)], [seq(a[i],i=1..n)],
                [seq(c[i],i=1..n-1)]] ):
X:=n->[seq(x[i],i=1..n)]:
F:=n->[seq(d[i],i=1..n)]:

Thanks

R dJ(t)/dt+J(t)/C=f(t)

where f(t) is a driving electromotive force. Use the fourier transform to analyze this equation as follows.

Find the transfer function G(alpha)  then find g(t) .

 Thanks ....

Hi,

I try to solve this equation using pdsolve but there is no results.
restart:
with(PDEtools):
with(plots):

Eq:=diff(u(t,x), t$2) =diff(u(t,x),x$2)+sin(u(t,x));

pdsolve(Eq);

Thank you for your help.

I have system of equation contain of two equations

eq1:=a1*x^2+b1*y^2+c1*x*y+d1*x+e1*y+f1:
eq2:=a2*x^2+b2*y^2+c2*x*y+d2*x+e2*y+f2: 

I want to find the number k so that the equation

a1*x^2+b1*y^2+c1*x*y+d1*x+e1*y+f1 + k*(a2*x^2+b2*y^2+c2*x*y+d2*x+e2*y+f2) = 0

can be factor, where k satisfy

a:=a1+k*a2:
b:=b1+k*b2:
c:=c1+k*c2:
d:=d1+k*d2:
e:=e1+k*e2:
f:=f1+k*f2:

I tried

A:=a*x^2+b*y^2+c*x*y+d*x+e*y+f:
collect(A,x);
B:=collect(discrim(A, x), y);
C:= discrim(B,y);

and got the expression

-16*(4*a*b*f-a*e^2-b*d^2-c^2*f+c*d*e)*a=0.

For example 

> restart:

a1:=14:

b1:=-21:

c1:=0:

d1:=-6:

e1:=45:

f1:=-14:

a2:=35:

b2:=28:

c2:=0:

d2:=41:

e2:=-122:

f2:=56:

a:=a1+k*a2:

b:=b1+k*b2:

c:=c1+k*c2:

d:=d1+k*d2:

e:=e1+k*e2:

f:=f1+k*f2:

P:=c*d*e+4*a*b*f-a*e^2-b*d^2

-f*c^2:

eq1:=a1*x^2+b1*y^2+c1*x*y+d1*x+e1*y+f1:

eq2:=a2*x^2+b2*y^2+c2*x*y+d2*x+e2*y+f2:

with(RealDomain):

Q:=solve(P=0,k);

factor(eq1+Q*(eq2));

solve([

eq1=0,eq2],[x,y]);

Where, Q is k which I want to find. 

My question is, if I have the system of equations

eq1:=a1*x^3+b1*y^3+c1*x^2*y+ d1*x*y^2 + e1*x+f1*y+g1:
eq2:=a2*x^3+b2*y^3+c2*x^2*y+ d2*x*y^2 + e2*x+f2*y+g2: 

How can I get a similar to the 

-16*(4*a*b*f-a*e^2-b*d^2-c^2*f+c*d*e)*a=0?

In that attached file is a multip step problem that involves graphing a right circular cylinder using transtion matrices and orthonormal basis. I have completed the hole question minus the very last part which is asking for new parametric equations for the cylinder if its center point is located at (-2, 10, 3) instead of the origin.

Any ideas on how to do this will be greatly appreciated.

question_4.mw

the sequence is non-decreasing up to some point after which it is non-increasing. Note that i can be 1 or n. A constant sequence is considered to be unimodal.

Examples of unimodal lists:

[1, 1, 1, 1, 1],
[1, 2, 2, 3, 4, 5, 5, 5],
[5, 5, 4, 4, 3, 3, 1],
[1, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 1],
[1, 2, 2, 3, 3, 3, 4, 4, 2, 2, 1, 1, 1]

Examples of lists that are not unimodal:

[1, 0, 1, 0],
[1, 1, 1, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 1, 1],
[1, 1, 2, 2, 3, 4, 5, 2, 2, 6, 4, 2, 2, 1, 0]

i don't have a clue

I am trying to make a substitution:

> restart;
> f1 := conjugate(z)^2*z^2;
> f2 := applyrule(z*conjugate(z) = x, f1);

however the result is

f2 := conjugate(z)^2*z^2

Do you have any ideas how to solve the problem? Thank you.