Below is the code found on this site for generating the mandelbrot set.

with(ImageTools):

MyMandel := proc(pt)

local i,z,c;

c := pt;

if abs(c) > 2 then

return 1;

else

z := c;

for i from 2 to 100 do

z := z^2 + c;

if abs(z) > 2 then

return i;

end if;  end do:

return 0;

end if; end proc;

GenerateComplex := proc(pts,X1,X2,Y1,Y2)

option evalhf;

local X::Vector, Y::Vector, f1::Vector, Z::Vector;

Forgive my ignorance here but let say we have some data that I want to

model with some partial differential equations. What would be most efficient way

to come up with the specific equations inorder to be able to calibrate such system ?

I guess I could just calculate the partial elasticities and deduct some relationships from that

but there has to be a more beautiful approach ?!

Okay, I'm trying to solve a person's expected utility given a budget constraint using the Lagrangian method on Maple. I keep getting the error 'Warning, solutions may have been lost'

Here is what I have so far. Can anyone tell me what I'm doing wrong?

Z:=(EU+λ*(Px*X+Py*Y+f))

solve({diff(Z,X)=0,diff(Z,Y)=0,diff(Z,f)=0,diff(Z,λ)=0}.{X,Y,f,λ})

Download 7845_notes.pdf
View file details

Hi,

I am trying to go through an example in my notes.

p:=2038074743;
q:=1190494759;
n:=p*q;

x:=141213562373095048;

y:=Power(x,2) mod n;

then I do the following,

Power(y,(p+1)/4) mod p;

how to find an expansion of: (we may assume m::posint)

(sum(a^k*b[k], k = 1 .. infinity))^m

for example: the trinomial theor...

(a+b+c)^m = sum((m,k1)*(m-k1,k2)*(m-k1-k2,k3)*c1^(k1)*c2^k2*c3^k3,k1+k2+k3=m)

here (m,k1) is the binomial coefficeitn

how to expand the following expression;

((1/6)*e*(6*(D(f))(gamma)+3*((D@@2)(f))(gamma)*e+((D@@3)(f))(gamma)*e^2))^m

Especially in old books (but not only there) the authors prefer the notion
of 'argument' ( = polar coordinates and using the angle).

Sometimes it comes to something like 

                | arg(-z) | < Pi and | arg(1 - z) | < Pi

which I prefer to have in terms of interval notation 'z in ... ' or its
complement  'not (z in .. )'.

I do not even know how to get | arg(z) | < Pi  <==> z is not a negative Real. 

How can I do it in Maple?

1. how to replace in the expression:

I'm trying to do an integral like this:

F(y) = int(x^-3 delta(1-y/x^1/2),x=1..2)

any advice?

Thanks

I spent the last 2 hours searching in the Help commands for the command that returns the length or size of a given ordered array or list.

e.g. I have Thing:=[3,7,6,5]:

I want the function that returns the value 4 for the size of Thing.

I am really surprised that such an operation is so difficult and hard to compute that it is never in any of the Help menus.

Hello,

I am working on a homework assignment for a Maple class and am hung up on creating adjacent lists within a list.  I know that I need to use the seq command with a i=1..(end) condition, but that code only creates the list, not the lists within the lists..

Here is the actual problem statement:

Consider the list

monthList := [ january, february, march, april, may, june, july, august, september, october, november december ];