## interpolate and plot larger data sets...

Maple needs commands to interpolate and plot large data sets. Consider the following:
```# generate data
data := convert(LinearAlgebra[RandomMatrix](10000,2,generator=-100000..100000),listlist):
data := sort(map([op], [op(op(table(map(`=`@op, data))))]), (a,b)->evalb(a[1] < b[1])):
data := Matrix(map(proc(a) [a[1],a[2]/200.] end,data), datatype=float[8]):

# now suppose I take this data and...
f := CurveFitting[BSplineCurve](data, x);
plot(f, 0..100);
```
Someone trying to do this with real data (and 10000 points is small) will experience the following problems:

## sparse linear algebra...

Maple needs proper sparse linear algebra. That means a data structure where row or column operations (good luck getting both) can be done in linear time and algorithms such as structured Gaussian elimination or iterative methods run in O(A) time, where A is the number of non-zero elements in the matrix.

## plot smoothing...

Maple 11's implcitplot command contains a cool option to reduce the line count of a plot. I would like a separate command for this, along with the corresponding generalization (level of detail) to 3d plots. I think it would be good to apply this sort of thing automatically to make plots faster and also to make larger more detailed plots possible. I also think that the 3d plot facility is far too limited in terms of the number of polygons it can handle. It should be replaced with something capable of "millions of triangles per second", because that's what any \$10 integrated graphics chip ca

## Descrete Math...

How do I use Maple to apply depth-first search to a graph to obtain a depth-first search of numbering of the verticies. Once again I need to input my graph in order to retrieve this information. How does Maple provide this info based off my graph?

## Integro-Differential Equations...

This is my system: rho * diff(diff(delta(x, y, t), t), t) = (y ^ 2 / 0.2e1 - d ^ 2 / 0.8e1) * Y * diff(diff(diff(diff(xi(x, y, t), x), x), x), x) + (((2 * G + lambda) * v - lambda) * diff(diff(xi(x, y, t), x), x)) + 0.16e2 * Q(t) ^ 2 * y * (d ^ 2 - 0.2e1 * y ^ 2) / d ^ 4 / b ^ 2 / l ^ 2 / e; diff(delta(x, y, t), y) = v * y * diff(diff(xi(x, y, t), x), x); i = 0.8e1 * Df * Q(t) / d ^ 2 + l * Lp * (int((2 * G + lambda * (v - 1)) * diff(diff(xi(x, y, t), x), x), x == (0..b)) + Y * d ^ 2 * int(diff(diff(diff(diff(xi(x, y, t), x), x), x), x), x == (0..b)) / 0.12e2) + diff(Q(t), t); My objective is to obtain the expression of delta,xi and Q to the space coordinates(x,y) and time,t.The other letters are constants. Can the maple do this? I try with pdsolve but obtain this error:

## Hoe to solve a boundary value problem...

I need to solve following ODE: ODE:= y(f)*(diff(y(f), f, f))+50.0*f -50.0*f^2= 0 Subject to boundary conditions: BC:= (D(y))(.5) = 0, y(0) = 0 I used BVP method : dsol2 := dsolve({ODE, BC}, numeric, method = bvp[midrich], y(f), range = 0 .. .5) But I got " Error, (in dsolve/numeric/bvp) division by zero" Can someone please tell me what I can do. Thanks, http://maplenet.maplesoft.com/maplenet/worksheet/mapleprimes/4384_BVP.mw

## Maple - easy talking to DLL...

Something I hadn't tried until today was calling into a DLL with Maple. I was surprised at how easy it is. Once inside a DLL, the sky is pretty much the limit as far as passing the data around between applications for realtime plotting and even passing to device drivers for other processing or even I/O in WinXP.

## querying the 'view' option of a graph...

Is it possible to get/set the view option of a plot component using the DocumentTools package? I'm trying to plot a function containing a parameter whose value is controlled by a slider component. When I adjust the view using rightclick->manipulator->scale and then move the slider, the graph reverts back to the hard-coded scale in the original plot command. I'd like to set a default scale, then allow a user to set the scale using the scale manipulator and have that new scale persist through the next slider manipulation.

## Last Year's (2006) Maple Conference...

Did anyone here attend the 2006 Maple Conference ? If so, is anyone willing to share some of the papers from the collected proceedings ? I would just like to see what new exciting developments were created from the global Maple community. v/r,

## markov model & plot a derivative function...

i create an availability function for my markov model in which use constraits like 2 maintenance tresholds and some maintenance duration, mean time between failures. i want to determine the mean time between inspections value which max. the availability and plot its graph... should i have to give exact numerical values or is there any other ways to solve it..

## Who is behind ScientificComputing blog?...

Does anyone know who the author of the Scientific Computing blog on blogspot?

## Relaxtion method to solve Ax=b...

implement a program for a matrix a a vector x and numbers eps>0, omega e(0,2) approximates the solution of the system Ax=b by relaxation with parameter omega and accuracy eps the program should stop when last iterations x[i-1],x[i] and x[i+1] satisfy norm(x[i+1]-x[i]> norm(x[]-x[i-1]

## Plotting 3D...

Hi, Sorry but i ve repeated the same question a few days ago but because it was quite down the line i thought to give it another shot I would like to make a 3D plot of some density functions(F(X)) where each lets say cdf, corresponds to a time value.So far Ive tried using the transform((x, y) command.The example below where I show how I have worked out for 4 cdfs corresponding to 10,20,30 and 40 ms, hopefully it will make it clearer of what i mean.(I have a thousand of different distributions up to 10 secs). with(stats): > C001:=plot(statevalf[cdf,lognormal[0.3792,0.9556]], 0..4, colour=green):