My goal is to obtain a formula F:= x -> fa(x)*c1_1(x)+ fb(x)*c1_2(x)+fc(x)*c2_1(x)+... where every ci_j(x) is a function of x. fa, fb, fc are known functions of x. I do not need to have it printed, I just need it to return a numerical value for every value of x I throw in.
I have obtained the solutions of ci_j's in the form
Sol[1]:=[c1_1=f11(x), c1_2=f12(x), c1_3=f13(x), ...], N1 terms
Sol[2]:=[c2_1=f21(x,c1_j’), c2_2=f22(x,c1_j’), c2_3=f23(x,c1_j’), ...], N2 terms
Sol[3]:=[c3_1=f31(x,c1_j’,c2_j’), c3_2=f32(x,c1_j’,c2_j’), c3_3=f33(x,c1_j’,c2_j’), ...], N3 terms
Sol[4]:=[c4_1=f41(x,c1_j’,c2_j’,c3_j’), c4_2=f42(x,c1_j’,c2_j’,c3_j’), c4_3=f43(x,c1_j’,c2_j’,c3_j’), ...], N4 terms

For testing an external routine against Maple I have something like
Gamma_mpl:=define_external(
'Gamma_mpl',
'C',
'Z'::ARRAY(1..n,datatype=float[8]),
'result'::ARRAY(1..n,datatype=float[8]),
LIB=theDLL);
where the function in the DLL (compiled with Dev-C++) is like
extern "C" __declspec( dllexport ) __stdcall
void Gamma_mpl(double *Z, double *R){ ... some code ...}
It takes the input in Z and updates values in R.
The following works as expected:
z:= -1.1 + 2.1*I*0:
# provide memory
Z:=Array(1..2, order= Fortran_order, datatype= float[8]):

How do I instruct a worksheet to automatically store a sequence of plots as .gif images?
In my worksheet, I generate a sequence of say 50 plots and then look at them within the worksheet.
I really need to store these someplace else on my computer as .gif files and do not want to export each plot individually.

I need to get the slope of this plotted line that I did, and I have no clue how. Here is the points I used:
points := [[22, 100.19], [40, 200.19], [50, 300.19], [62, 400.19], [75, 500.19], [95, 600.19]]
I'm not sure what do to, I tried slope[[22,100.19],[95,600.19]] so try and get the slope of the line between those two points, but nothing came up, it just went to the next line and said "slope".

I want to write a maple program/script involving matrix computations. The program will consist of approximently 50 statements and will require user supplied input data. I have Maple 9 and would like to call the script from Matlab.
Matlab would provide the input data to the Maple script and then I would like to have Maple pass the output data to Matlab. Is there a way to do this using Maple 9 or do I need the new Maple Matlab toolbox?
Thank you,
Richard

Hello together,
I've just started to do some calculations with maple 9.5 and still in trouble. I would like to modify the step function below to a n-periodic function. In maple version 10 there's a function called "rept" which is not available in my version unfortunately. Does anyone have a ease of use solution?
f:=(MTBF,MTTR,t,n) ->
piecewise(t<>

Hi,
I have a small problem.
I want to use the Schurform computation in the LinearAlgebra-package. The main objective is to compute the stabilizing solution to a Algebraic Matrix Riccati equation. (perhaps someone knows some code for this?)
I call the SchurForm function as
(T,Z) := SchurForm(A,output=['T','Z'])
where A has the correct type (ie complex(sfloat)). This works well in a Maple sheet, but in a procedure (where T and Z are defined local) the kernel says that the output chosen is invalid. If I don't specify any output, it also works good (but then I get no Z). The Maple version is 10 and I am using Linux.

On a whim I decided to write a routine to numerically calculate the volume of a sphere and compare the result with the well known formula (4/3)*Pi*r^3. What struck me was how slow the MAPLE code ran compared to the same code in VB. I realize that a program like MAPLE which is geared toward symbolic mathematics would not be expected to be as fast at numerical calculation as other vehicles, but I was struck by the difference. Below is the code in VB and MAPLE. Both routines give the same answer. The MAPLE roputine runs about 10 times slower than the VB routine. I was wondering - perhaps I'm not being fair to MAPLE. Is there some way I could have written the code to run significantly faster?

I have data of the type xyz-coordinates. How do a fit a line to it?
Thanks

I was trying to use the LinearAlgebra package to discover the eigenvalues and eigenvectors of the 2x2 matrix representing a rotation in the x-y plane.
The matrix [[cos(theta)|-sin(theta)],[sin(theta)|cos(theta)]]
However the function, Eigenvectors does not work with this matrix, and the function Eigenvalues returns the eigenvalues of the trasformation but not in the usual complex number notation.
Can anyone help with explaining and or resolving this ?
v/r,
Dan

Hi I have tried installing Maple 10 on both by Desktop and laptop both of which are running XP SP2. It fails in both cases giving an error about an unknown user interface.
"Unable to preprare and load the installer in GUI mode"
Any ideas guys?
Cheers,
Moz!

I have the following maple file. Is it possible to instead of having the incoming rays at a vertical at intervals along the x axis, whether to have them emerging from a single point on the y axis, and also whether this point can moved up and down the y axis. Thanks

View 3259_mirrorplot.mw on MapleNet or

Download 3259_mirrorplot.mwView file details
I am trying to use the Mean command of the Statistics package to calculate the arithmetic mean of a transformed random variable symbolically.
I am using Maple 10.
In my case the Mean command gives an incorrect answer.
First a simple example, which gives me a correct result symbolically:
restart;with(Statistics):
X:=RandomVariable(Uniform(-Pi, Pi));
Mean(X+a);#it gives a, it is ok.
#and here is my more complicated problem:
restart;with(Statistics):X:=RandomVariable(Uniform(-Pi, Pi));
assume(Gt::RealRange(Open(0),Open(1)),Gr::RealRange(Open(0),Open(1)));
Mean(abs(1-Gt*Gr*exp(-I*X)));

Hi,
I have three equations Y=f(t) with each equation corresponding to a probability as shown below.
Y1=f(t) with probability of 0.05%
Y2=f(t) with probability of 5% and
Y3=f(t) with probability of 50%
If we assume that these probabilities belong to a normal distribution how can I work out the equations for 99.95% and 95%probability.In other words I would like to mirror Y1 and Y2.
The idea is to plot all the equations of 0.05%,5%,50%,95% and 99.95% probability in the same plot and create a 3D surface where the axis will be Y,t and P(probability) or even better generate the equation of P=f(Y,t) that will describe this surface.

I have a list of arrays which I would like to merge into a single array.
In the case of two arrays (say A and B) I basically want maple to dump them both into a new array in which wherever A[i,j]<>B[i.j] the new array takes whichever is nonzero and if A[i,j]=B[i,j] does nothing.
is there any short way of doing this without writing a function which runs through the arrays and compares each element induvidually and outputs the results into a new array?