MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

want to write in maple code

to generate all commutative diagram 

with adjacency matrix

 

however, i only know a -> b, b-> c , a->d , d-> c

google no information about all commutative diagram, 

 

another problem is

would like to enrich theory , however, do not know how to connect property such as equations with diagram

Hello!

I am trying to write the following function in maple:

f(P(x)) = (x-1)*P'(x)-x*P(1) 

Where P(x) = 1 + x + x^2 

 

My approach thus far has yielded limited results; i am able to get some, but not all, of the results expected.

This is my code,( the relevant bits ):

p := x-> 1+x+x^2

f:= P -> (x-1)*P'(x)-x*P(1); 

if i then type f(1); in maple i get the expected result of -x

if however i type f(x) i get an equation that is unable to produce the expected result when simplified - seemingly, it doesn't parse for maple. The expected result would be -1 

Using unapply make's no difference. Either i have made a mistake typing in my functions - in which case i would like an ëxplanation as to what i did wrong, or seemingly there is a different approach here? Perhaps solving this problem with maple is simply not doable? - anyway my line of reasoning leads me to finally ask how i write a function within a function in this nice program :) 

Any help in this matter is greatly appreciated - i only just recently started using this program. 

I have the following d.e.:

I need to change the s variable into a different one, where the new variable is defined by

(the old s shows up in the limit of the integral)

I tried dchange, but it chokes on this as I don't have an explicit representation of s in terms of Theta.

(I know the overall solution as other people smarter than me have solved this a long time ago, but I 'd like to have the derivation to understand it).

Mac Dude

Let A and B be matrices of dimensions 3,2 and 2,3 corresondingly. It is
known Multiply(A,B)=Matrix([[8, 2, -2], [2, 5, 4], [-2, 4, 5]])


 What is Multiply(B,A)?

limit(t*(int(exp(-t*tan(x)), x = 0 .. (1/2)*Pi)), t = infinity)?

Is it possible to find it in Maple? The command

MultiSeries:-limit(t*(int(exp(-t*tan(x)), x = 0 .. (1/2)*Pi)), t = infinity);
outputs
.

 

Dear All,

I have problems connected with a double integration in Maple.

The simplified task with its description is in the uploaded file.

Any help will be appreciated.

Thank You! 

 

mp_qn_double_integra.mw

Dear all,

I would like to solve the Fredholm Integral equation, using numerical method.
This is my code.

there is a problem with subs, does not working here.
# Then, we obtain from (9) the coeficient A[n] and B[n].

Then I woulk like to recompute (2), and then compute (1).
# Puting x=m*h, in (1), how can we generate a linear Matrix from (1).

 

Fred.mw

Thanks

f:=C->sum(x[cat(new,c)],c=1..C);

f(4) # does not give me what i want , which is

x[cat(new,1)]+x[cat(new,2)]+x[cat(new,3)]+x[cat(new,4)];

 

seq(x[cat(new,c)],c=1..4);
convert([%],`+`);

This works. But it's just a simple example to illustrate my problem.

 

Is there a way to fix this 'cat' problem in the procedure 'f'?

 

Thanks,

 

casper

infolevel[dsolve]...

April 10 2014 mcypert0725 5

Hi everybody,

I am using maple to solve differential equations.  I have been using maple for some time and came across the infolevel[dsolve] command or option.  I tried finding the answer to my question online but cant seem to come across a straight answer.

My question is what does the integer do that you have to equal the infolevel[dsolve] command?  

For example I am using a book that states

"An even more important diagnostic tool is the infolevel[dsolve] command, which will give information on what methods are used in attempting to solve the ODE when dsolve is applied, even if unsuccessful.  An integer between 1 and 5 must be specified, with generally more detailed information being provided as the number is increased. On applying the dsolve command to ode2, the method of attack is summarized in the following output and, in this case, the general solution y(x) given with two arbitrary coefficients C1 and C2.

What does the integer between 1 and 5 mean or do? In the example above the book uses 5, why 5?.

Thanks,

Matt 

Dear all;

Please in this code Maple, I plot using pointplot3d, some points. I would like to add a name of each point.

The name of each point  P:=[x(i),y(j),t(k)] is u[i,j]^k.  Please how can I add the name at each point. here is the maple code. 

Thanks.

pointplot3d.mw

v1:=[seq(tau[lambda,i]^2,i=1..5)];

latex(v1);

 

gives:

[{\tau_{{\lambda,1}}}^{2},{\tau_{{\lambda,2}}}^{2},{\tau_{{\lambda,3}}
}^{2},{\tau_{{\lambda,4}}}^{2},{\tau_{{\lambda,5}}}^{2}]

which produces a non natural output in latex, see below.

I wonder it could be improved in the future release?

'latex' command is really useful when i have large algebric output. But having to change this tiny thing at many places is just very time consuming. I'd really appreciate if this can be fixed in updates or future versions of Maple.

 

Casper

 

Hi all,

I will use the following dummy example.

with function,f

f:=(xid,yid)->sum(x[i],i=1..xid)*sum(y[i],i=1..yid);

 

and a complicated term, myterm


myterm:=(f(3,4)+f(2,2))*f(1,1):
myterm:=expand(myterm);

 

'if' i have some previous knowledge, or know a bit of the term, i can find the structure by doing this


repar:=[f(1,1),f(2,2),f(3,4)];  # Or with more f(xid,yid) terms

tmp:=seq(repar[i]=ff[i],i=1..3);
simplify(myterm,{tmp});          # This is fine, gives me what i want

 

But, can we go further, and more 'obvious'

 

Given the fucntion f, same as before, and the same 'myterm'

can I have this
restart:
iwant:=(f(3,4)+f(2,2))*f(1,1);  # as a result, straightforward

so I dont have to go back to 'repar' and find that the terms exactly are.

 

Thanks,

 

Supressing Warnings...

April 09 2014 xcyborg 15


restart; printf("Factory Simulation: %s\n\n", Ghana*Chocolates); printf("Chocolate Factory Simulation Project .\n\n"); printf("%5a  %10a  %10a  %10a  %10a  %10a   \n", W, O(t, t+1), B(t), T(t), F(t), R(t)); printf(" ---------------------------------------------------------------------------\n"); _local(T, R, F, B, O, Q); initialise := proc (n) local i, t, a; T := proc (t) options operator, arrow; 250 end proc; R := proc (t) options operator, arrow; 150 end proc; F := proc (t) options operator, arrow; 250 end proc; B := proc (t) options operator, arrow; 50 end proc; for i to n do t := i; if B(t) < F(t) then Ds(t, t+1) := B(t) else Ds(t, t+1) := F(t) end if; M(t, t+1) := 50; O(t, t+1) := 50; P(t, t+1) := 50; if t = 2 then O(t, t+1) := 100 elif t = 3 then O(t, t+1) := 0 end if; P(t, t+1) := T(t)-F(t)+Ds(t, t+1); if R(t) < P(t, t+1) then P(t, t+1) := R(t) elif R(t) < 0 then P(t, t+1) := 0 end if; R(t+1) := R(t)+M(t, t+1)-P(t, t+1); F(t+1) := F(t)+P(t, t+1)-Ds(t, t+1); B(t+1) := B(t)+O(t, t+1)-Ds(t, t+1); R(t) := R(t+1); F(t) := F(t+1); B(t) := B(t+1); for a to t do Q(t, t+1) := O(t, t+1)+O(t-t+a, t-t+a+1) end do; T(t+1) := (t+1)*Q(t, t+1)/t; printf("%5a  %10a  %10a  %10a  %10a  %10a   \n", i, evalf[3](O(t, t+1)), evalf[3](B(t)), evalf[3](T(t)), evalf[3](F(t)), evalf[3](R(t))) end do; printf(" ---------------------------------------------------------------------------\n\n\n"); Restart; printf("%5a  %10a  %10a  %10a     \n", W, M(t, t+1), P(t, t+1), Ds(t, t+1)); printf(" ---------------------------------------------------------------------------\n") end proc; initialise(10)

    1         50.         50.        250.        250.        150.   
    2        100.        100.        200.        200.        200.   
    3          0.          0.        300.        300.         50.   
    4         50.         50.          0.          0.        400.   
    5         50.        100.        125.        125.        325.   
    6         50.         50.        120.        120.        280.   
    7         50.         50.        117.        117.        283.   
    8         50.         50.        114.        114.        286.   
    9         50.         50.        112.        112.        288.   
   10         50.         50.        111.        111.        289.   
 ---------------------------------------------------------------------------


    W          50      875/18          50     
 ---------------------------------------------------------------------------

 

NULL

NULL


Download suppress_warnins.mw

Hi All,

 

I need to suppress the warning messages when i try to pring my Table. Any help will be much apreciated.

 

Regards,

Xcyborg

I’m trying to figure out how to find a basis for a subspace, V, of Rdefined by V = {(x, y, z)l(2x-3y+6z=0)}

 

I’m using the student linear algebra module for maple 17

 

I’ve tried defining the subspace and asking for the basis of V but I always get an error code.

 

I’ve tried consulting the maple website and looking through their help menu, but can’t find anything that answers how to find a basis... At least a basis from the subspace defined in my problem.

I know how to find a basis for the subspace by hand but not with maple.

Any help will be greatly appreciated. 

(1/2)*(-x-y+sqrt(-3*x^2-6*x*y-3*y^2))/(x+y)

 

the complex value is always a constant whatever x and y 

i suspected that this is a wrong function

5 6 7 8 9 10 11 Last Page 7 of 996