I'm working on a 3D program that plots stars, star clusters and other objects of interest against a model of the Milky Way. The guts of this program create a few hundred points (< 2000) that represent stars, and move them a little bit in random directions to create a cloud effect. Then the animate procedure rotates them about the center of the galaxy. I would like to add more stars and do more things with them, but I am encountering a very sluggish editor, suggesting that Maple is imposing a memory limit. The Activity Monitor says that Maple is using 140 MB, so gigabytes of RAM remain available. Furthermore, if I save the plot as an animated gif, I can run half a dozen copies simultaneously in separate browser windows without encountering noticeable slowdown. But I can't elaborate the program without better performance from the editor.

The program is written in the Document interface. I tried saving it to the classic interface, but that didn't seem to help, and I found it difficult to edit.

Ideas sought.

The problem has been discussed in mapleprimes before, but I have not come across a neat solution. Perhaps there is one. What is the state of the art for exporting 3D graphics (plot3d, etc.) with Maple 15?

What I think I know is this: standard GUI has a completely new way of treating 2D graphics but has an outdated way of rendering 3D, while the classic GUI, though no longer developed, for some reason typically renders better 3D graphics.

I used to export graphics...

> restart;
> with(LinearAlgebra);
> with(linalg);
> with(ArrayTools);
> n := 2; p := vector([2, 3, 3]);

> x := vector(n, []);
> x0 := vector(n, []);
> st := vector(n, []);
> for i to n do x0[i] := 0; st[i] := -1 end do;
> f := (x[1]-p[1])^2+p[3]*(x[2]-p[2])^2;
                          ...

Dear all,

I have encoutered a funny problem with Optimization package. I minimize some long stuff (lets call it A) and get an answer, e.g.

[.116553990209052046e-19, [n = 4.24568325257496060, rhocut = 1.66125593854788400, m = 2.00000021862965971, rho = .956914761577446150]]

then i just do

evalf(subs(n = 4.24568325257496060, rhocut = 1.66125593854788400, m = 2.00000021862965971, rho = .956914761577446150, A));

and i get

I hope I've got the terminology right...

I've got a set of functions that I want to display in a maplet table. The tables shows up fine, except all the functions are written like:

x^2+1

What I'd like is them to be shown like:

I've tried playing with MathMLViewer, but that doesn't seem to work.

> restart;
> with(LinearAlgebra);
> with(linalg);
> with(ArrayTools);
> n := 2;
> x := vector(n, []);
> x0 := vector(n, []);
> st := vector(n, []);
> for i to n do x0[i] := 0; st[i] := -1 end do;
> f := proc (x) options operator, arrow; x[1]^2-4*x[1]+x[2]^2-2*x[2] end proc;
         2                2 ...

Has anyone been able to successfully install the Maple-Linbox package?

Basically when trying to make the maple configure linbox package one gets errors like:

 error: no match for call to ‘(LinBox::BlasBlackbox<LinBox::Modular<double> 

There is already a duplicate post on the linbox support forums about this w/out any answer:

https://groups.google.com/d/topic/linbox-use/5ILTopngYGc/discussion

Hi all,

I want to solve the following 2 odes simultaneously:

ode1:=(1/r^2)*diff(r^2*diff(u(r),r),r) =v(t); (#boundary value problem with 2 boundaries)

ode2:= diff(v(t),t)=u(r); (#initial value problem)

with conditions

IBC={u(0.2)=1,u(0)=0.5,v(0)=0.01}

I am getting error if i use dsolve/numeric. Is there any way/method to solve these numerically using maple packages like dsolve/fsolve...

Thanks,

I tried a partial differential equation example from a Robert Lopez example.  The answer shows a nice smooth graphical representation of heat transfer. 

However why does my method here, using numeric pdsolve, not produce the same results? 

Also pdsolve using numeric does not use the default timestep of 1/20.  Explicitly telling Maple to use timestep of 1/20 produces a different graph telling me that the default value is not being used in pdsolve.

Hello,

I have recently been toying around with nonlinear fitting for as part of a research job and I've been having problems getting nonlinear fit to:

1) work

2)create a good fit.

Here is my input:

> with(Statistics); with(plots); with(LinearAlgebra);
> list1 := Vector([0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70,
 75, 80, 85, 90, 95, 100, 105, 110, 115, 120, 125, 130,
 135, 140, 145, 150, 155, 160, 165, 170,...

In coding theory, there are parity-check codes whose parity-check matrices H are generated via column permutations. For instance, the LDPC codes constructed in Gallager's 1962 IRE Trans paper uses the following H matrix:

[ X1 ]

[ X2 ]

[ .... ]

[ Xn ]

where submatrices X2 .. Xn are just random column permutations of X1. However, to make the codes efficient in decoding, there is one restriction which requires that any two row vectors in...

Hi,

After having Maple solve a large system of equations, I'd like to convert the resultset into functions programmatically.

Something along these lines:

for p in resultset do functionname=lhs(p); functionname=unapply(p); od;

This of course doesn't work, this assigns the function to 'll'. I tried something like this as well:

for p in resultset do functionname=lhs(p); subs(FN=functionname, FN=unapply(p)); od;

Hi all,

 I am having a very hard time with numerical stability. I am solving system of ode's (7-coupled ode's) using dsolve(stiff) and then using spline function for interpolation and finally solving system of pde's (2-coupled pde's) using pdsolve for one time step and solving all again for the next step. the solution is not stable and it requires very fine/small time step. Is there any procedure/method to improve the stability? 

 What is the stability criteria of dsolve...

I have to evaluate several complex integrals involving square roots. Maple automatically assume the branch cut is all negative reals. Is there a way to change this? This would make my life alot easier! thanks!

I want to get the formula for frequency at resonance peak (   ω_n=ω*sqrt(1 - 2*ζ²)   ) in Maple 15, starting from the transfer function


                      G(s) := 1/(1+2*zeta*s/omega[n]+(1/(omega[n])^2)*s^2);


I wrote these lines but I haven't obtained the results I wanted:

                    'abs(G(j*omega))'=subs(s=i*omega,abs(G(s)));