conisder

>eq1:=n=(k[f1]^3+k[f2]^3)/6/Pi^2;


> eq2:=el*V=h_bar^2/2/m*(k[f2]^2-k[f1]^2);


> eq3:=solve({eq1,eq2},{k[f1],k[f2]});

> eq4:=allvalues(eq3):

Solving this i get a very very massive result. Let me take a branch:

eq5:=subs(eq4[2],k[f1]);

I get

1/h_bar*(-(1/2*n*Pi^2*h_bar^3/el/V/m+1/2*...

What is the next two terms for this pattern...

1, -2, 2, -4, 0, ...

Hi/

the code below wroten by Maple 13. I use h=0.29669 and N:=15 and digit=80 fsolve works perfect! but when I changed N:=16 and another number greater than 15 fsolve didnt return result and return some equation.

whould you help me how to resolve this problem?

sorry i upload this code. to run easilySinc_MBVP_V2-temp.mw

thanks

> restart;
> N := 15;

I am playing around with the PolynomialFit command in maple.
If I want to get a return of 1% then I can simply plug in the return at time t
for v in the expression below. That works for all different returns.

Hello,

(sorry if this appeared elsewhere, but it seems my first try to ask this failed)

I tried to use Maple last night to do some heavy computational work, but it turns out (I let him print time()) that he was only active for about an hour.  He did not do anything when my Macbook Pro was on stand-by.

This morning, when he woke up, he was unwilling to stop!  I pushed the "Interrupt the current operation" button...

seriestoratpoly(series(simplify(LaguerreL(5, x)),x,6),['egf']);

seriestoratpoly(series(simplify(LaguerreL(5, a, x)),x,6),['egf']);

inverse z transform of sqrt(1-a^2)/(1-a/z)

restart;
alpha := 0;
s := x;
tau := alpha + 1 - x;
source := s*Diff(y(x),x$2)+tau*Diff(y(x),x)+n*y(x)=0;
sol := dsolve(source,y(x),method=laplace);
zsol := ztrans(sol,x,z);
convert(zsol, ratpoly);

sol := y(x) = _C1*invlaplace((_s1-1)^n*_s1^(-1-n), _s1, x)
zsol := ztrans(y(x), x, z) = _C1*invlaplace((_s1-1)^n*_s1^(-1-n)*z/(z-1), _s1, x)
Error, (in convert/ratpoly) expecting a Laurent series or a Chebyshev series

what are _C1 and _s1?

When I execute several times my maple code the program completely freezes, even though in the beginning of the code teh command "restart" is inserted. I work with matrices, perform integrations and derivations for their components and calculate eigenvalue problems. The freeze happens at no particulat block of execution.

Any thoughts?

I need program of multiple shooting method with maple or algorithm?

Thanks alot.

When I use ImportMatrix to import a MATLAB mat file, I get a list of names and values. Should I write a "for" loop to assign them in Maple or there is an easier way to do it?

I am working over F_p, where p is prime. When p = 2, the number of 2x2x2 arrays (aka hypermatrices or tensors) over this field is 2^(2*2*2) = 256. Of those 256 arrays, I only want the symmetric ones: that is, if x_{ijk} is the ijk-th element of the array X, then x_{ijk} = x_{ikj} = x_{jik} = x_{jki} = x_{kij} = x_{kji}. Is there a quick loop that does this? 

Question 1. After googled, x*diff(F,t)+t*diff(F,x) for two variables, how about three variables or general ?

Question 2.

if deal with on variables case using the example in

http://demonstrations.wolfram.com/AFunctionInvariantUnderAGroupOfTransformations/

F := (x^2-x+1)^3/((x^2)*(x-1)^2)-y;
G:=y*diff(F,x)=0;
sol:=solve(G,y);

F(H)-F = 0
sol:=solve(F,x);
H := simplify(subs(y=F, sol[1]));

all above code can...

I would appreciate a hint why Maple 16 Classic under Win XP 32 bit crashes due to the following sequence of commands

> dir1:="C:\\Maple16\\P16\\test1\\";
> path:=cat("dir ",dir1);
> ssystem(path);

with the lasst ssystem command doing the crashing while the same sequence of commands under Maple 13 Classic does not: Instead, system returns to Maple 13 worksheet the content of the directory "dir1" as expected.

Is there another way to issue that ssystem command?

hello,

after defining N_1 N_2 and N_3 with long numbers, and defining M by N_1*N_2*N_3 i tried to step backwards.

By entering something like: M/N_1 (what is exactly like N_2*N_3) it does not show me full numbers, but instead 1/something*M.

After refreshing M and afterwards refreshing M/N_1 it totally worked out and showed me the product of N_2*N_3

Now my question is: how can i prevent this odd display format without refreshing everytime? (i tried...

similar to

http://www.singular.uni-kl.de/Manual/html/sing_951.htm
http://www.singular.uni-kl.de/Manual/html/sing_753.htm

i would like to find all modules in the complex
for example

0 -> A -> B -> C -> 0

would like to find B and C if know A

which book describe this algorithm for polynomial system