Hello,

I have defined a function(say fun) to do certain task, it works well when I call it from a script file(say script1). However, when I copy the complete script1 and paste onto another script file(say script2), the function, fun, doesn't execute. It's very peculiar as my project requires lot of modification, so I copy the old script and...

Hello every one,

I am gettting this error "unable to store HFloat", while solving a system of ODES numerically.

restart:

eq1:=n*(-diff(f(eta),eta$2))^(n-1)*diff(f(eta),eta$3)+(2*n/(n+1))*f(eta)*diff(f(eta),eta$2)

-diff(f(eta),eta$1)^2-M*diff(f(eta),eta$1)+M*epsilon+epsilon^2+lambda*f1(eta)=0;

eq2:=diff(f1(eta),eta$2)+P*(2*n/(n+1))*(f(eta)*diff(f1(eta),eta$1))=0;

bc:=f(0)=0,D(f)(0)=1,D(f1)(0)=1,D(f)(N)=epsilon,f1(N)=0;

If S= { <0,1>,<1,2>} and T={<1,1>,<2,3>} are ordered bases for R² and v=<1,5> and w=<5,4>, then can someone help me:

-find the coordinate vector of v and w with respect to the basis T by creating a procedure?

-I can make a procedure to find the transition matrix P_S<-T from T to S-basis, but can you help find the coordinate vector of v and w with respect to the basis S. Using P_S<-T?

Hi guys, I would like some help writing a procedure that checks whether the group of input vectors is an orthogonal basis. Any help would be great!

How does one change a previous Atomic setting back to non-Atomic?  This seems so simple, but I can't find anyway to Undo this change.  

Thanks,

Dear Mapleprimes,

I am trying to solve an Optimization problem analytically. I set up a lagrangian as follows:

L:= Objective function + mu_1 (constraint 1) + mu_2 (constraint 2)

I maximise with respect to two variables, call them x1 and x2

I want constraint one to bind but not constraint 2. Hence, I set mu_2=0. The I do:

diff(L,x1); diff(L,x2), diff(L,mu1), diff(L,mu2)

Hi!

I want to claculate a triple integral of a piecewise function. The piecewise function has 17 domains and I guess that this is a bit too much.

Is there any possibility to optimize it?

the piecewisefunction looks like this:

a := unapply(piecewise(
    x1 < 0 or x2 < 0 or x3 < 0 or h < x1 or h < x2 or h < x3, 0,
    0.5*h <= x1 and x1 <= h and h-x1 <= x2 and x2 <= x1 and h-x1 <= x3 and x3 <= x1, A1[13],

Hello,

I am trying to force Maple to read big matrices from files, produced by other programs (GAP in my case).

The output looks like this after it's been saved in a *.txt file:

[ [ 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
      0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
      0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0,

I have install Maple13 and GRTensorII on my PC (Ubuntu12.04) ~ but GRTensor dosent work ! When i start Maple13 ,there are the words "Error, "/home/herb/.mapleinit" is an empty file". 

Is anybody can help me out of it? I really have done the following in README:

3. Go to the INSTALLDIR. In this directory you will find a file called
 mapleinit.sample. This is a Maple configuration script which has

I am Beginner in maple , i cannot see why while this function cannot be defined at x=1, but the graph is plotting fine and continuesely in x=-2..2 ?!

and with this one we get another graph and yet different one !

sin((1-x)^0.5)/(1-x^2)^0.5

Can anybody explain to me why , your help will be appreciated !

 I need help to plot  a sequence of the form given in the attachment below     

seq.docx    

As usual I am going to use this blog post to complain!
How can anyone in their right state of mind create a Finance Package
with out procedures for return calulation ie first difference, percentage change, log returns etc

Hi,
I'm trying to import matlab function from m-file.

File looks like this:

    function out = xxx(omega)
    out = heaviside(omega);

It's only heaviside function.

And now, when I'm using "FromMFile" command in Maple it doesn't recognize Heaviside function. It imports that function as "heaviside" - witch small first character (in Maple we needs big "H").

Any help? 

thx 

I am trying to write a script that will read in vector files using ImportMatrix, perform a simple calculation then output a matrix file using ExportMatrix.

The difficulty I am having is that I would like to do this on a number of files,

eg.

for i from 1 to 600 do

Mi:=ImportMatrix(rowi):

do calc....

ExportMatrix(matrixi,Mi):

end do:

The variable i is in italics...

I want to model this game (I want an ensemble average and a time average)

50% chance you get head=win 50%
50% chance you get tails= lose 40%
Starting capital =100

First of all I need a very large binary matrix(1..100000,1..100000) however not
even that I can do fast enough in maple!

I tried:

rtable(1 .. 1000, 1 .. 2000, frandom(0 .. 1, .5), subtype = Array, datatype = float[8], fill = [0, 1])

it is somewhat fast but it does not do what I want ie binary ...