Hello!

I have a systematic  and 100% reproducible error  by using Maple 16 in my calculations:

Error,(in evala/Factors) Bad inputs to diophant

Running the same job with Maple 15, 14,   I don't have such an error.

Would it be possible to correct something in Maple 16 and to make it useable?

If you need code of my program I will send it to you.

Best regards,

     ...

I'm trying to figure out the forces on this arm crane when alpha angle is changed from 0 degrees, 30 degrees, 45 degrees, 60 degrees and 90 degrees. What I'm confused about is how to input the angle in maple. Do I add 45 degrees and alpha angle together?

Apparently the load (L) is not unknown but I'm not sure how its not unknown when I am not given any number for it. How do you find out the load (L)?!

Here's the picture:
http://postimage.org/image/x4s13q3b7/

Hi,

I am new to Maple and I am trying to write a CFD code. I want my code to be as modular as possible, therefore, I want to use subroutines  in my code. (Because I will need to do some steps several times )

Here is one of the main subroutines that I want to call for it in various times with different parameters. I am familiar with matlab so let be show what I mean first in matlab notation.   Lets say my subroutine is called MAT 

I am very new to Maple. I learned how to graph polar coordinates. The next thing I want to learn how is to highlight the area (not to calculate, but to show) between two polar graphs. Can you tell me what is the code for it? This is my set up:

with(plots);

polarplot([2+cos(x), 2], x = 0 .. 2*Pi, color = [red, blue]);

Thanks.

I'm trying to add the output of my function FunctionToTree,which gets all elements from an expression, to a graph using the Graph function but I'm getting erros even when I add it as a list.

with(GraphTheory):
with(SpecialGraphs):
FunctionToTree:=proc(expr)
local elem;
local Tree:=NULL;
local exprType := convert(whattype(expr),string);
local newExpr:=convert(expr,string);
if (nops(expr) <> 1) or (nops(op(1,expr)) <> 1) then
for elem in op(expr) do :

Hi, I am having a little trouble finding the roots of Bessel's function J0 between -4 and 4 using the newton Raphson method, the main issue I am having is not knowing what to input at a certain point.

f:=x->x-evalf(cos(x));
                     f := x -> x - evalf(cos(x))

> low:=-10;

                              low := -10

> up:=10;

                               up := 10

> step:=1;

                              step := 1

Hi I am new to using Maple,

Mybe this is a very basic question.

I am trying to find solutions to the following nonlinear system of equations using fsolve.

This is the right hand side of a system of ODE.

nonlinearsystem.mwnonlinearsystem.mw

I do get two solutions for the given parameter values. My question are:

As discussed in http://www.mapleprimes.com/questions/144586-Op-Command-Isnt-Working-As-Expected , the indices change some time at run time so previously working command is not working for me.

I am trying to use 

eq3 := applyop( u->Parts( u, phi[i](y) ), { 1,10,13,20 }, eq2 );

which clearly use the 'op' to get operands from eq2, but since the indices change at every run or on some run, the code is useless. The problem is...

Hi all,

I am new to Maple and I have (possibly a simple) question. I'm solving 2 coupled differential equation, it all goes well, until I want to plot them.

I have the output in the following form:

sol = {w1(x) = [ -a*x^4+b*x^3-c*x]}, {w2(x) = [-d*x^4+e*x^3-f*x^2+gx-h]}

My goal is now to plot the solution in 1 figure with w1(x) between x=0..100000 and w2(x) between x=100000..150000

I have searched on the web for a way to get the formula's of w1(x...

Hello

I am trying to plot solution of ode0 together with the maximum and minimum values but I am having difficulty since the first plot is a solution and second is values. I should have a plot with two line one represent the solution of ode0 and second (the max an min). Any advise or suggestion?

This is the code:

> restart;
with(DEtools); with(plots); Nsols := 5; Ntstep := 10;
 k := 0; A := 0.37e-1; B := 0.2e-6;
ode0 := diff(U(t), t) = -(A+B*U(t))*U(t);

((x-6.05)/sqrt(.0004/9) < (6.035-6.05)/sqrt(.0004/9))

Maple comes up with P(150.000000000x<905.2500000), but I cant get Maple to tell me

P(Z<-2.25).. How can I get it to calculate and tell me to use Z tables ?

I am trying to get the kron of two vectors as in;

N=6;

p:=Vector([$0..N]):
pn := KroneckerProduct(p, transpose(p));

but I get this error ==>

Error, invalid input: LinearAlgebra:-KroneckerProduct expects its 1st argument, A, to be of type Matrix, but received Vector(7, {(1) = 0, (2) = 1, (3) = 2, (4) = 3, (5) = 4, (6) = 5, (7) = 6})  

Why I can not get kron of two vectors? 

I don't understand why 'expand(x+y)^18' the smart popup work, but with 'expand(x+y)^19'  (or more then 19) i have no smart popups!

Any idea?

Thank you all

Hi 

I am a Matlab user and just started using maple because I need more precision for my calculations. I failed to make Matlab interface work, therefore, I am writing my codes from the begining. And I am having a problem with a nested for loop.

The thing that I am trying to make is (in Matlab)

for j=1:N;
    for i=mod(j-1,2):2:j;
D(i+1,j+1)=2*j;
    end
end

In maple I tried something like

Hi Maple Community!!

Using dsolve, Mathieu 15 gives me the following solutions:

Given that I have specific values for b, e and g. How can I plot x(t...