Good morning
I'm looking for the best way to create an array with four index, let's call it V_{ijkm}, knowing that it is completely independent from the order of the index (eg V_{1234}= V_{2134}= V_{4213} =etc) and that I must use nested for cycles to initialize it, instead of a indexing function
I can't guess how to insert the knowledge that the order of index in unimportant in the code
Thank for any hints
s.v.

Hi boys

Here I'm trying to solve this problem.

#A bar of 1meter of length is placed in a room at 10ºC. At a time t=0 its borders (x=0 and x=1) is placed in contact with Ice at 0ºC and then the system is kept isolated. Show the temperature profile of the bar at each time until t=0.2s#

I want to display the output from printf ( columns as seen below ) in a MathContainer (embedded component) . How ?

seq(printf("%20a %10a %10a %10a %10a %10a \n", i, i^2, i^3, i^4, i^5, rt[i]), i = 1 .. 10);

Dear all,

I'm trying to obtain a 3D graph with Maple starting with 3 experimental variables : IC50, REDOX, LIPO

Here are my list of points and the procedure I wrote.

Hi
I've remarked something strange in a program I wrote.
It's a mix of analytical and numerical computation, so that I use a value of Digits higher than the default, to improve and stabilize the results against floating numbers errors.
What happens is that for Digits from 10 to 15 included, the time the PC needs to run the program increases reasonably, so that for Digits=15 it needs 10 minutes or so.
But if I set Digits=16 the time rises exponentially, and becomes more than two days!!

January 21 2010
EUIX 44
I am having problems making histograms in maple. It seems that the default Maple command plots a Histogram (that is the command Histogram), with which the area of the bars is overall normalized to 1. How do I get rid of this feature? If I plot the data set (1,2,3,4,5,6,1,2,3,4,5) with a binwidth of 1. I should expect 6 bars each of count 2 except the last 1. However in maple I get something completely strange using the command Histogram(DATASET, binwidth=1). I am using Maple 12.

1) compute the first ten derivatives of f(x)=sin(x)cos(x) using a for...do...od loop. Describe the pattern for the derivatives. Use your description to predict the 2500th derivative form.
2) Expand the do-loop code in the first question to include a command line inside the loop that plots each derivative. The output shows that each derivative is plotted separately. To plot the derivatives simultaneously after the loop has terminated, save the plot inside the loop using a label such as pl[i]:=plot(df[i],x=-10..10_.

Hi, all,
I want to sort hundreds of variables which are read from a file.
Basically, this is what I want to use:
Basis(L,plex(V));
"L" is a variable of list which is read froma file. This part is done.
"V" is a variable of list which is read from a file. Well, this part is not done yet. Because I have no idea how to do this.
For example, in a file, there are some variables: a, b ,c,d,e,f,g
and I can read them into a variable of list: V=[a,b,c,d,e,f,g]; But the problem here is you cannot use it as plex(V) in Groebner basis package.
Any suggestions?
Thanks in advance.
Gepo

January 20 2010
djc 546
Hello,

I 'd like to use j as imaginary unit, and I want to display it in the title of a plot using *typeset*, but it doesn't work (in the expression F):

restart ;interface(imaginaryunit=j);sqrt(-1);

F:=2*sin(w)/w*exp(-j*w*2);F;

plot(abs(eval(F,w=2*Pi*f)),f=0..10,title=typeset('F(j*w)'=F));

Thanks for help in advance.

January 20 2010
kayj 20
Hi All,

I am new to maple and I am having problems writing to a file, the problem that i am having is that maple rounds

the small numbers into 0 when I write into a file but does not round if it outputs to the screen,

how can I get the screen output written in a file? Also, is it possible to write to a comma delimited file ".csv". the code is below

I apprreciate your help

##############################################

file := fopen("table.txt",'WRITE');

f:=(k,m,a)->expression;

Someone has implemented Levenberg-Marquardt algorithm. That is implemented in the Dr. DE Holmgren's procedure in

http://www.maplesoft.com/applications/view.aspx?SID=4159

www.maplesoft.com/applications/view.aspx

But it returned the error: __Error, (in ln) numeric exception: division by zero__ for the function a * x ^ b.

Gracias

HerClau

January 20 2010
emoh 24
I have an area delimited by y=sqrt(x), y=0 and x=4.
I want to rotate this area around the y axis to get the solid of revolution, but instead I get another solid.
I am using:
VolumeOfRevolution(sqrt(x),x=0..4,axis='vertical','distancefromaxis'=0,'output'='plot');
If you look at the resulting solid, you'll see that is the product of rotating the area limited by y^2(reverse function of sqrt(x)), x=0 and y=2 around the y axis, but not of the area delimited by y=sqrt(x), y=0 and x=4.
How can I plot it with VolumeOfRevolution command?

Hi,

i want to fit the parameters of a function for a set of points in 2D, problem is, the function is not explicitly given as y=f(x), but in parametric form.

Example:

restart:with(Statistics):

N := 200:

U := Sample(Normal(0, 1), N):

X := <seq(cos(2*Pi*i/N)+U[i]/5, i = 1..N)>:

Y := <seq(sin(2*Pi*i/N)+U[i]/15, i = 1..N)>:

ScatterPlot(X, Y);

Lets say i want to fit a,b of the parametric model function [a*cos(t),b*sin(t)], t=0..2*Pi

I have now used quite some time trying to find out how to display an isotope in Maple in 2D-math notation. In He-4-2 I want the nucleon (or mass) number 4 to stay as a superscript, but to the left of He, and the proton number 2 as a subscript also to the left of He. How can this be done in Maple in 2D-math notation? I hope someone can help.

Erik V.