Hi all,

I am having some problems when solving high order polynomials. When oders are 'high' and systems involve multiple variables, Maple is slow to do so.

See in worksheet for detailed codes.

slow.zip

kappa - gives the vector that i am interested in

sC - a reparameterisation which has 'summation' in it

Aim is to verify that it works (it should!)

I currently ask...

Hi,

How do I get ride of these Rootof?

I tried simplify,evala,value,Simplify and ect. Didnt really find anything useful.

Download rootof.mw

 

I dont care which root they actually take, all I want is one of the roots. So I can then use subs for substitution.

Casper

Hi all,

What else could I try if I encounter this error?

Warning:  Incomplete separation.

--------------------------------

 

I read the pdsolve help page, where it says I should have a solution with this warning, but I got nothing but the error.

 

ptr_mplprims.mwHello,I am trying solve  a travelling force problem with pdsolve, but I cannot obtain the second derivative of the solution. I obtained the first derivative, with help from another post (http://www.mapleprimes.com/questions/80997-Derivative-From-Numeric-PDE-Soln#comment81512). When I tried doing the same thing for the second derivative it didn't work. Any sugestions would be helpful. Thank you!!!

Say I have the equation:

f:=x^5+a*x^4+b*x^3+c*x^2+d*x+e

and I want to introduce an identity namely

x^3:=1

Since

subs(x^3=1,f)

or

eval(f,x^3=1)

doesnt work I wonder what is the most elegant way to implement this?

Why does maple not simplify this expression any further:

-4*cos((1/5)*Pi)*m^3-4*cos((1/5)*Pi)*n^3-4*cos((1/5)*Pi)*p^3-4*cos((1/5)*Pi)*q^3+4*cos((2/5)*Pi)*m^3+4*cos((2/5)*Pi)*n^3+4*cos((2/5)*Pi)*p^3+4*cos((2/5)*Pi)*q^3-10*cos((1/5)*Pi)*p*q^2+10*cos((2/5)*Pi)*m^2*n+12*cos((2/5)*Pi)*m^2*p+12*cos((2/5)*Pi)*m^2*q+12*cos((2/5)*Pi)*m*n^2+10*cos((2/5)*Pi)*m*p^2+12*cos((2/5)*Pi)*m*q^2+12*cos((2/5)*Pi)*n^2*p+10*cos((2/5)*Pi)*n^2*q+12*cos((2/5)*Pi)*n*p^2+12*cos((2/5)*Pi)*n*q^2+12*cos((2/5...

Hi,

Say for a 126 by 18 symbolic matrix, substituting numerical values into the matrix DD1 takes about 10 seconds (using time() );

for a 62 by 16 matrix, take just about 1/5 second.

By when the matrix is about 254 by 20, it take much longer! More than a minute!

I never actually waited how long, but I am wondering if there is something wrong? Am I using the right command?

Say the matrix is DD1 with numpar:=[x=0.3,y=0.2,z=0.98,....]

How do i solve this PDE with initial and boundary condition

1. δθ/δζ = δ^2θ/δy^2

At y=0, δθ/δy=0

At y=1, δθ/δy=-µ

At ζ=0, θ=1

Sorry to bother you again. I want my computations(results) to carry/have more decimal places so that I can choose how many decimal places to show myself. This is because, my result is always displayed in 9 d.p but I want more decimal places to be displayed. Just like the way we change the way maple display matrix in-line with command interface(rtablesize= N), is there a similar way i can do the decimal place thing?

Secondly, how can i suppress the value A amd B from...

How do I plot the imaginary or real part of a complex valued function

e.g.

E:=y^2+x^2=1

I tried:

E1:=eval(E,{y=y1+I*y2,x=x1+I*x2})

implicitplot3d(E1,x1=-3..3,x2=-3..3,y1=-3..3)

but which does not work

Hi

I want to simulate taking a card from to decks at the same time, and by doing it sufficiently many times, show that the probability is 1-(1/e).

I've implemented this in Python, but i'm new to maple. Here is what I've done so far.

 

randomDeck:=proc(n::integer)

  local deck:= [seq(1..52)];
local randomDeck:=[];

Hi,

I need to model the MapleSim ground as a hard ground. For example when the robot swing leg hits the ground, it should not have negative vertical position. Basically, the robot will move on a surface.

Thanks...

Hi

 

In earligere versions of maple (14 and 15) I yoused this code (program) to plot the area between functions, but it does not work in maple 17 (or 16).

integralplot2:=proc(f1,f2,a,b,{[color, colour] := blue,transparency := 0.5,numpoints:=500})     

local mingr,maxgr,variable,p,c,n,pts1,pts2,polys,i,xkoor,ytop,ybottom,yvals,plotf1,plotf2;          

if evalf(a)<evalf(b) then       

Hello I have the following system of equations and I want to solve for rho1 and z1

E := (1/2)*(1/rho1^2+1/rho1^2)-1/sqrt(z1^2+rho1^2)-1/sqrt((R-z1)^2+rho1^2)-1/sqrt(z1^2+rho1^2)-1/sqrt((R-z1)^2+rho1^2)+1/R+1/sqrt((2*rho1)^2+(R-2*z1)^2);

Er1 := diff(E, rho1);

Ez1 := diff(E, z1);

solve([Er1, Ez1], {z1, rho1});

The system is a function of R.

The solve command spits out Kernel lost after some time. So my question is how can I solve it ?