> 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)));

I have a problem.

I don´t know how can i do for working with a dynamic model equation was done in maple in Maplesim?

Help me please

At one point in earlier versions of Maple I was able to plot expressions such as y=4*a*x and the y axis would be incremented in units of a. Is there a way to do this in Maple 14?