hi all

is there any (possibly simple) way to distribute a piece of maple software not allowing external user to access the source code ?

Translating to C, compiling and then including as an external function is very awkward procedure.

thanks in advance

Hello!

I'm trying to solve numerically an ODE system with piecewise. And this piecewise is very important for this task.

This system describes behavior of a pulley with friction. There are some constants: m, c, g, mu and J. Values of this constants are not important.

> sys := m*a(t) = piecewise(a(t) < a0, F0*time, a(t) >= a0, 0), v(t) = diff(x(t), t), a(t) = diff(v(t), t);
> m := 5; F0 := 10; a0 := 5;
> initialconditions := x(0) = 0, v(0) = 0;

Dear folks,

 In mapple programming while assigning a fraction to a variable...it directly simplifies it....

eg: x:=4/6 and it directly gives 2/3

can we make mapple not to simlify this fraction....

It is very important for me to find out....

please help me

thanks in advance

regards

uday

I have looked at the maple help for the ErrorPlot function, but i can't make it work for my plot. 

I have the following procedure:

K := 1.; Xrng := [seq(dt*i, i = 0 .. ntime)];

value := array(1 .. ntime+1);

ave := array(1 .. ntime+1);

for j to ntime+1 do:

discountfactor := exp(-r*j*dt);

Stot := 0.;

Stot1 := 0.;

for i to npath do:

Spath := max(Spot[i][j]-K, 0);

Stot := Stot+Spath;

What does it mean when I get the message, "Warning, model is not of full rank"?

Thank you.
Janice

Hi all,

Is there some way i can find the width of a plot?

I have the equation below:

Prob:=(w,T)->(B^2)*(sin(w*(T/2))^2)/((w/2)^2);

plot(Prob(w,6),wf=-3..3);

Which produces a gaussian looking plot, the plot reaches zero at some...

mem2.mw

Although i have defined K5 as a constant before solving the equations , maple returns an error message saying that

"Error, (in fsolve) k5 is in the equation, and is not solved for"

And i am alos trying t find the roots of the equations in a certain range as specified above.

Dear all,

Find here a part of my precedent post "2D finite element method" which represent a worksheet for 2D triangular Mesh.A part of the mesh can rotate with regard to the other by modifying theta.

restart: with(geometry): with(plots):

ms := 8: n := 6: mr := 5: theta := (1/180)*(23*Pi):

wrs := (1/1000)*<30.8, 33, 36, 39, 40, 43, 46, 49, 53.25>:

wrr := (1/1000)*<20,23,26,28,30,30.8>:

wt := (1/180)*Pi*<0,10,20,30,40,50,60>:

Dear All!

Could you explain how the coulditbe() function works?

> assume(q, integer, w, integer); about(q, w); coulditbe(q/w < 1); coulditbe(q < w);
Originally q, renamed q~:
  is assumed to be: integer

Originally w, renamed w~:
  is assumed to be: integer

                              FAIL
                              true

BR,

Zoltán Faigl

In a theme similar to the movie Matrix a little procedure to create falling random number animation, and in this example 0's and 1's hence the name binary rain.

MatrixRotaterows := proc (arows, bcolumns, randstart, randend, iterations) 
local a, i, b, c, d, e, f, g:
a := Matrix(arows, bcolumns, 0):
for i to iterations do
  b := Matrix(1, bcolumns, rand(randstart .. randend)):
  a := ArrayTools:-CircularShift(a, 1, 0):
  a[1,..] := b:

Good Morning.

I have problem when i want to graph my equations. I received this error: 

Error, (in plots/odeplot) curve is not fully specified in terms of the ODE solution, found additional unknowns {t}

I would like to know how i can solve it

Thanks in advance 
 

Dear All,

I ran into a problem, which I can not resolve. 

> Eqs := [x*0 = y*0]; Vars := [x, y];
                            [0 = 0]
                             [x, y]
> solution := solve(Eqs, Vars);
                        [[x = x, y = y]]
> Eqs := [x*b = y*b, b = 0]; Vars := [x, y];
                       [x b = y b, b = 0]
                             [x, y]
> solution := solve(Eqs, Vars);
                               []

Hi, my professor wants our graphs to be labeled (which is reasonable). How can I label the graph I made this way:

>inequal({x[1] >= 0, x[2] >= 0, x[1]-3*x[2] >= -3, 2*x[1]+3*x[2] <= 6}, x[1] = 0 .. 3.5, x[2] = 0 .. 2.5)

Thanks.

Hi all,

I noticed that there are not many applications which deal with the finite element method with Maple. I attached in the file below a code for the magnetostatic probems and which allows the calculation of the magnetic field in a permanent magnet synchronous machine.

mafem.zip

I hope this work will help many users of maple in the numerical analyzis with the finite element method

I posted...