> local `+`;
Error, unable to parse

meet difficulty running script in maple 12

restart;  local `+`;  `+`:=proc(a,b) :-`+`(a^`~2`,b^`~2`) end proc;

how to override plus + operator so that apply on matrix operation such as linearalgebra package

for example if + =  a^2+b^2

1 + 2 result in 5

and also apply on matrix multiplcation for each element

Hi all,
I am writing a procedure that I would like to have several outputs. The code is as follows:

Initialize := proc (p, theta, PiV, PiU, n, m)

J := CompJ(PiU, m);

for i to n do

if member(i, V) then

Ai[i] := CompAi(PiU, PiV, i, n, J)

end if

end do

end proc

The function calls other functions previously defined that are not important here. Basically I would like to return the set J and the list Ai, but I am not sure on how to do it. 

Once the procedure returns these elements, I would like to assign them to some variables. I am not sure on how this can be done when a function returns more than one thing.

Thanks in advance for the help
Bests
Manuele 

I have a system of linear differential equations and am trying to solve them using Fourier transforms.

I can reduce the system to a final result (for the variable of interest) in Fourier space as this (note the frequency variable is 'w'):

restart:

with(inttrans):

vout_fourier_num := fourier(phi[3](t), t, w) = 6.63569999999998*10^(-15)*w^2*fourier(V(t), t, w)/(-5.69875218358308*10^(-40)*w^4+(9.19473390627057*10^(-29)*I)*w^3+2.15219369729956*10^(-18)*w^2-(4.14691648617110*10^(-8)*I)*w-700.8);

#the drive can be defined as:
drive:=5.70000000000000*10^(-6)*exp(-3.18877551020408*10^18*(t-2.0*10^(-9))^2)*cos(4.8*10^9*Pi*(t-2.0*10^(-9)))*10^(19/20);

#substitute the drive in - this is not necessary, but it should work!...
vout_fourier_num2:=subs(V(t)=drive, vout_fourier_num);

#now take the inverse... note this gives 0!...
invfourier(vout_fourier_num2, w, t);

the final results calculated is zero. It is wrong... it seems like an accuracy issue, but increasing the digits does not help. I should note that I can calculate the solution directly via dsolve, and get completely reasonable answer.

any ideas how to get the Fourier method to work?

thanks!

hello

this is my program and fsolve for low intensity solve the equations but for high intensity cannot solve why?

this is my code:

ep0 := 1/(4*3.14);

el := 8.54*10^(-2);

hbar := 1;

vf := 1/300;

kb := 1;

tem := 2.586*10^(-2);

ci := 1;

p := 1.458*10^16;

beta := 2;

ai := 7.1*10^(-4);

bi := ai/sqrt(3);

enph := .196;

d := enph/(kb*tem);

n0 := 1/(exp(enph/(kb*tem))-1);

gama := hbar*vf;

intensity:=9000000

w := 7.28;

impurity := 7.2*10^3;

g := hbar*beta/(bi^2*sqrt(2*p*enph));

aa := g^2*(n0+1)/(2*Pi*hbar*gama^2);

bb := g^2*n0/(2*Pi*hbar*gama^2);

cc := 2/(Pi*gama^2);

l := (1*hbar)*w/(2*kb*tem);

u := el^2*intensity/(32*w*hbar^2);

[fsolve({op([((enph*ln(1+exp(c+enph/(kb*tem)))/(kb*tem)-polylog(2, -exp(c))+polylog(2, -exp(c+enph/(kb*tem))))*enph*(kb*tem)^2-(enph^2*ln(1+exp(c+enph/(kb*tem)))/(kb^2*tem^2)+2*enph*polylog(2, -exp(c+enph/(kb*tem)))/(kb*tem)+2*polylog(3, -exp(c))-2*polylog(3, -exp(c+enph/(kb*tem))))*(kb*tem)^3+(-exp(b)*enph*ln(1+exp(c+enph/(kb*tem)))+exp(c+d)*enph*ln(1+exp(b-d+enph/(kb*tem)))+exp(b)*kb*tem*polylog(2, -exp(c))-exp(c+d)*kb*tem*polylog(2, -exp(b-d))-exp(b)*kb*tem*polylog(2, -exp(c+enph/(kb*tem)))+exp(c+d)*kb*tem*polylog(2, -exp(b-d+enph/(kb*tem))))*enph*(kb*tem)^2/((exp(b)-exp(c+d))*kb*tem)+(exp(b)*enph^2*ln(1+exp(c+enph/(kb*tem)))-exp(c+d)*enph^2*ln(1+exp(b-d+enph/(kb*tem)))+2*exp(b)*enph*kb*tem*polylog(2, -exp(c+enph/(kb*tem)))-2*exp(c+d)*enph*kb*tem*polylog(2, -exp(b-d+enph/(kb*tem)))+2*exp(b)*kb^2*tem^2*polylog(3, -exp(c))-2*exp(c+d)*kb^2*tem^2*polylog(3, -exp(b-d))-2*exp(b)*kb^2*tem^2*polylog(3, -exp(c+enph/(kb*tem)))+2*exp(c+d)*kb^2*tem^2*polylog(3, -exp(b-d+enph/(kb*tem))))*(kb*tem)^3/((exp(b)-exp(c+d))*kb^2*tem^2))*bb+u*(1/(1+exp(-l-c))-1/((1+exp(-l-c))*(1+exp(l-b))))-(((1*enph)*(enph-2*kb*tem*ln(1+exp(-b+enph/(kb*tem))))/(2*kb^2*tem^2)+2*kb^2*tem^2*(-polylog(2, -exp(-b+enph/(kb*tem)))+polylog(2, -cosh(b)+sinh(b))))*enph*(kb*tem)^2-(enph^2*(enph-3*kb*tem*ln(1+exp(-b+enph/(kb*tem))))-6*kb^2*tem^2*(enph*polylog(2, -exp(-b+enph/(kb*tem)))+kb*tem*(-polylog(3, -exp(-b+enph/(kb*tem)))+polylog(3, -cosh(b)+sinh(b)))))*(kb*tem)^3/(3*kb^3*tem^3)-(-exp(b)*enph^2+exp(c+d)*enph^2-2*exp(c+d)*enph*kb*tem*ln(1+exp(-b+enph/(kb*tem)))+2*exp(b)*enph*kb*tem*ln(1+exp(-c-d+enph/(kb*tem)))+2*exp(c+d)*kb^2*tem^2*polylog(2, -exp(-b))-2*exp(b)*kb^2*tem^2*polylog(2, -exp(-c-d))-2*exp(c+d)*kb^2*tem^2*polylog(2, -exp(-b+enph/(kb*tem)))+2*exp(b)*kb^2*tem^2*polylog(2, -exp(-c-d+enph/(kb*tem))))*enph*(kb*tem)^2/((2*(-exp(b)+exp(c+d)))*kb^2*tem^2)-(exp(b)*enph^3-exp(c+d)*enph^3+3*exp(c+d)*enph^2*kb*tem*ln(1+exp(-b+enph/(kb*tem)))-3*exp(b)*enph^2*kb*tem*ln(1+exp(-c-d+enph/(kb*tem)))+6*exp(c+d)*enph*kb^2*tem^2*polylog(2, -exp(-b+enph/(kb*tem)))-6*exp(b)*enph*kb^2*tem^2*polylog(2, -exp(-c-d+enph/(kb*tem)))+6*exp(c+d)*kb^3*tem^3*polylog(3, -exp(-b))-6*exp(b)*kb^3*tem^3*polylog(3, -exp(-c-d))-6*exp(c+d)*kb^3*tem^3*polylog(3, -exp(-b+enph/(kb*tem)))+6*exp(b)*kb^3*tem^3*polylog(3, -exp(-c-d+enph/(kb*tem))))*(kb*tem)^3/((3*(-exp(b)+exp(c+d)))*kb^3*tem^3))*aa-u*(1/(1+exp(l-b))-1/((1+exp(-l-c))*(1+exp(l-b)))) = 0, -cc*polylog(2, -exp(b))+cc*polylog(2, -exp(-c))-impurity = 0])}, {op([b, c])})];

thank you.

hi,I want to solve this equation with the following boundary condition numerically by maple:

hello evreybody i have these Error :

restart:with(plots):
M:=765 : m:=587 :I:=76.3*10^3 :Jp:=7.3*10^3 :e:=10.92: F:=0.42: omega:=0.56 :ka:=0.1:kb:=0.1:kc:=0.1: lambda1:=0.1 :lambda2:=0.1:lambda3:=0.1:
Error, illegal use of an object as a name

please help 
thank you !

Is there away to quickly execute a specific section within Maple? I know that I can execute a selection, but that sometimes requires me to select several lines of Maple code and can be tedious and tiresome.

Hello, 

    I'm finding Maple very useful. I'm a student but not very far along, only in MATH-1075 (intermediate algebra).

I have learned how to plot equations using the context menu, and the smartplot() command.

I have two questions I have not been able to discover an answer for.

When I plot an equation, how can I ask Maple to tell me coordinates of X & Y intercepts?

When I plot 2 equations, how can I ask Maple to tell me coodinates of intersections?

thanks in advance for any help. I tried to look up myself how to do this but was unsuccessful :(

Brent

Hi every body:

i will plot this function in maple in domain x=0..1,but when use the plot commond i see this warning:

Warning, unable to evaluate the function to numeric values in the region; complex values were detected

how can I plot it?

eq := 4.728139762*10^(-14)*sqrt((-5.905404581*10^83*x^4-1.382542004*10^88*x^2+2.271592177*10^92+1.631888578*10^79*x^6+2.795756904*10^23*sqrt(-1.317535223*10^121*x^8+6.370084310*10^125*x^6-7.852926774*10^129*x^4+2.223707894*10^132*x^2-1.592721566*10^134))^(1/3)*(-(1.114348319*10^53*I)*x^4-6.433693022*10^52*x^4+(2.688368867*10^57*I)*x^2+1.552130489*10^57*x^2+5.072945110*10^26*x^2*(-5.905404581*10^83*x^4-1.382542004*10^88*x^2+2.271592177*10^92+1.631888578*10^79*x^6+2.795756904*10^23*sqrt(-1.317535223*10^121*x^8+6.370084310*10^125*x^6-7.852926774*10^129*x^4+2.223707894*10^132*x^2-1.592721566*10^134))^(1/3)-3.723246850*10^61-6.448852713*10^61*I+(1.732050808*I)*(-5.905404581*10^83*x^4-1.382542004*10^88*x^2+2.271592177*10^92+1.631888578*10^79*x^6+2.795756904*10^23*sqrt(-1.317535223*10^121*x^8+6.370084310*10^125*x^6-7.852926774*10^129*x^4+2.223707894*10^132*x^2-1.592721566*10^134))^(2/3)+1.231339558*10^31*(-5.905404581*10^83*x^4-1.382542004*10^88*x^2+2.271592177*10^92+1.631888578*10^79*x^6+2.795756904*10^23*sqrt(-1.317535223*10^121*x^8+6.370084310*10^125*x^6-7.852926774*10^129*x^4+2.223707894*10^132*x^2-1.592721566*10^134))^(1/3)-1.*(-5.905404581*10^83*x^4-1.382542004*10^88*x^2+2.271592177*10^92+1.631888578*10^79*x^6+2.795756904*10^23*sqrt(-1.317535223*10^121*x^8+6.370084310*10^125*x^6-7.852926774*10^129*x^4+2.223707894*10^132*x^2-1.592721566*10^134))^(2/3)))/(-5.905404581*10^83*x^4-1.382542004*10^88*x^2+2.271592177*10^92+1.631888578*10^79*x^6+2.795756904*10^23*sqrt(-1.317535223*10^121*x^8+6.370084310*10^125*x^6-7.852926774*10^129*x^4+2.223707894*10^132*x^2-1.592721566*10^134))^(1/3)

with regards...

 Hello everybody, I need help please   

sign is different runnning the same script in maple 12 and maple 15 in different machine, just using matrixmultiply and matrix(xxx, shape=hermitian)

assume my window 7 infected by virus, can maple prevent virus change

its accuracy?

or

is there any change between maple 12 and maple 15? 

i feel that i will have to make large effort in order to find reason such as reinstall maple 12.

window 7 maple 15

[[a = 3.720799777 10 , b = -3.720817167 10 ]], 0.3469023622
[(-0.737729376724384 + 0. I) ((-0.737729376724384100 + 0. I) a

+ (-0.738311510115612690 + 0. I) b) + (-0.590656319609631

- 0. I) ((-0.590656319609630831 + 0. I) a

+ (-0.589459378339369122 + 0. I) b) + (0.326925800179577

- 0. I) ((0.326925800179577230 + 0. I) a

+ (0.327770888454982090 + 0. I) b) = -0.736196608749071 + 0. I,

(-0.590656319609631 + 0. I) ((-0.591008582207233624 + 0. I) a

+ (-0.589459378339369010 + 0. I) b) + (0.331003458223746

+ 0. I) ((0.331003458223746439 + 0. I) a

+ (0.327770888454981812 + 0. I) b) + (-0.735625969091165

- 0. I) ((-0.735625969091165288 + 0. I) a

+ (-0.738311510115613134 + 0. I) b) = 0.327869065042947 + 0. I

]

window 8 maple 12

[[a = 0.2249293777, b = 1.221244758]], 0.9888653482, "************"
[(0.737729376724384100 + 0. I) ((0.737729376724384100 + 0. I) a

+ (-0.738311510115612690 + 0. I) b) + (0.590656319609631053 - 0. I)

((0.590656319609631053 + 0. I) a + (-0.589459378339369122 + 0. I) b) +

(0.326925800179576676 - 0. I) ((0.326925800179576676 + 0. I) a

+ (0.327770888454982090 + 0. I) b) = -0.736196608749071002 + 0. I,

(0.590656319609631053 + 0. I) ((0.591008582207233624 + 0. I) a

+ (-0.589459378339369010 + 0. I) b) + (-0.331003458223746605 + 0. I) ((
-0.331003458223746605 + 0. I) a + (0.327770888454981812 + 0. I) b) + (
-0.735625969091165178 - 0. I) ((-0.735625969091165178 + 0. I) a

+ (-0.738311510115613134 + 0. I) b) = 0.327869065042946218 + 0. I]
7

please help me to find an analytical approach to the below equation:

> ode3 := diff(n(t), t)+(1/2)*(-(3.707186000*(0.815e-1*(diff(n(t), t, t))+diff(n(t), t)))/(0.815e-1*(diff(n(t), t))+n(t))^(3/2)-(.1428*(1+0.714e-1*n(t)))*(diff(n(t), t)))/sqrt(7.414372/sqrt(0.815e-1*(diff(n(t), t))+n(t))-(1+0.714e-1*n(t))^2)+n(t)+sqrt(7.414372/sqrt(0.815e-1*(diff(n(t), t))+n(t))-(1+0.714e-1*n(t))^2)-(2.518891688*(1+.3570*n(t)))*sqrt(0.815e-1*(diff(n(t), t))+n(t)) = 0;
                                                                           /
                                                                           |
/ d      \                                 1                               |
|--- n(t)| + ------------------------------------------------------------- |
\ dt     /                                                           (1/2) |
               /           7.414372                                2\      |
             2 |------------------------------- - (1 + 0.0714 n(t)) |      |
               |                          (1/2)                     |      \
               |/       / d      \       \                          |       
               ||0.0815 |--- n(t)| + n(t)|                          |       
               \\       \ dt     /       /                          /       
              /       / d  / d      \\   / d      \\
  3.707186000 |0.0815 |--- |--- n(t)|| + |--- n(t)||
              \       \ dt \ dt     //   \ dt     //
- --------------------------------------------------
                                     (3/2)          
           /       / d      \       \               
           |0.0815 |--- n(t)| + n(t)|               
           \       \ dt     /       /               

                                        \       
                                        |       
                              / d      \|       
   - 0.1428 (1 + 0.0714 n(t)) |--- n(t)|| + n(t)
                              \ dt     /|       
                                        |       
                                        |       
                                        /       

                                                           (1/2)
     /           7.414372                                2\     
   + |------------------------------- - (1 + 0.0714 n(t)) |     
     |                          (1/2)                     |     
     |/       / d      \       \                          |     
     ||0.0815 |--- n(t)| + n(t)|                          |     
     \\       \ dt     /       /                          /     

                                                             (1/2)    
                                   /       / d      \       \         
   - 2.518891688 (1 + 0.3570 n(t)) |0.0815 |--- n(t)| + n(t)|      = 0
                                   \       \ dt     /       /         
> ics := n(0) = 0, (D(n))(0) = 674.5142595;

thanks and regards

louiza

How to create animation f(x)=x^2, g(x)=sqrt(x), x=0, and x=1 if rotation about line y=-1 and x=-1?

Thank for your attention and helping.

Hello everybody,

I am trying to solve a partial differential equation (wave equation) with initial conditions that contain derivative and x-varying values

My equation
Uxx – Utt = 0 (hiperbólica) where xx and tt denote second partial derivatives with respect to x and t
0 <= x <= 10; t > 0;
u(0,t) = u(10,t) = 0;
Ut(x,0) = 0;
u(x,0)=
x for 0 <= x < 5
10-x for 5 <= x <= 10

I ve defined my equation as follow
> PDE := diff(u(x,t),x,x)=diff(u(x,t),t,t);

I m first trying without the x-varying initial condition.

> IBC := {u(0,t)=0,u(10,t)=0,u(x,0)=10-x,diff(u(x,0), t) = 0};
Returns the ugly

And that gives error when trying to solve
> pds := pdsolve(PDE,IBC,numeric);
Error, (in pdsolve/numeric/process_IBCs) initial/boundary condition does not contain the dependent variables of the PDE: 0 = 0

Another try with
> IBC := {u(0,t)=0,u(10,t)=0,u(x,0)=10-x,diff(u(x,t=0), t) = 0};

Looks like more promissing but returns (in a solving attempt)
> pds := pdsolve(PDE,IBC,numeric);

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions must depend upon exactly one of the independent variables: diff(u(x, t = 0), t) = 0

Regarding the x varying initial condition, I ve tried the `assuming` keyword in various attempts but without success.

Anybody could help on the matter?

Thanks in advance for anytip.
Have a good day
Nicolas