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,

thanx for ur answer.it reallly work for float data.but some of the data are complex,like:0.0536472874I,and some are float like:0.043443322222

the float numbers are written in the file and the program is worked until the result is complex number.it stop the program and print the error again!

what should I do with this?!

thank you

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    

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 

Hello, I'm creating a rolling thin disk (rolling without slip) in the MapleSim, the model like the picture below (Fig. 1).

       ...

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 ...

I'm trying to solve the following

> de := proc (omega, l, M) options operator, arrow; diff(y(r), `$`(r, 2))+(2*r-2*M)*(diff(y(r), r))/(r*(r-2*M))+(omega^2*r^2/(r-2*M)^2-l*(l+1)/(r*(r-2*M)))*y(r) = 0 end proc;

How to use maple get the nth term of sequence A: 7 + 77 + 777+ ... to n terms.

please help me!!

Thanks!

I am trying to fit some data to a model.

The model is given by this equation.

A*(C/B)^(2*B)*(B-1)^(2*(B-1))/(205*10^9*(2*B+1)*(sigma^2-(650*10^6)^2)*sigma^(2*(B-1)))+(C-sigma)^2*A/(205*10^9*(2*B+1)*(sigma^2-(650*10^6)^2))

> with(Statistics);
> X := Vector([819.4, 795.6, 788.8, 782.0, 776.56, 763.64, 748.28, 724.42, 717.40, 711.28, 707.20, 680], datatype = float);
> Y := Vector([5.3*10^4, 7.8*10^4, 9*10^4, 9.5*10^4, 10^5, 1.2*10^5, 1.37*10^5,...

hello,

I have been asked to show that the parametric representation of this curve is given by 

x(t)=cost/1+sin²t and y(t)=sintcost/a+sin²t.

The curve is defined by (x²+y²)²=x²-y² 

The polar representation of this curve is given by r²=cos(2theta)

Help!

Thankyou