I am working on this d2.pdf       problem,is ther any solution with maple?

 The coefficients are function of (T) and I would like Si's is also a function of T.

The well known(!) Heun's equation is,

df:=diff(y(z),z,z)+(gamma/z+delta/(z-1)+epsilon/(z))*diff(y(z),z)+(alpha*beta*z-q)/(z*(z-1)*(z-a))*y(z) = 0;

The solution is the general Heun's functions and it should be written as HeunG according to the help file of Maple 14 but my solution is found in terms of some function called HG which is not listed in the help system.

1. Are they the same?

I am using display command to get implicitplot. But I only get the graphs in one red colour. I want to get graphs in different colours. How to do that? 

How does one get the InequalityGraphics package? I don't see it as a standard package within Maple 14.

If I open help for pdsolve/numeric in Maple14 I can read

PDEsys - single or set or list of time-dependent partial differential equations in two independent variables

How can I solve numerical PDE system with more than two variable? If I try it, I get error message

Error, (in pdsolve/numeric/process_PDEs) can only numerically solve PDE with two independent variables, got {r, t, t0}

Hey everyone.

I have a wierd problem when using the "Units[Standard]" addon.

I'm trying to do some simple thermodynamics calculations, and need to raise certain values in an irrational power.

It looks like this:

restart:
with(Units[Standard]):

n:=1.38:
p1:=1.1[bar]:
v1:=3[dm^3]:
v2:=0.4[dm^3]:

p2:=(p1*v1^n)/(v2^n):

If I execute it, I get this error:
Error, (in Units:-Standard:-^) a unit can only be raised to a rational power

I wish to solve the following equations simultaneously:

(diff(g(t), t, t))*f(x)-(diff(f(x), x, x))*g(t) = 0

and

(diff(f(x), x))^2*g(t)-(diff(g(t), t))^2*f(x) = 0

I searched the help files and did the following:

I have tried using solve but get an error

solve({`&x`(diff(g(t), t, t), f(x))-`&x`(diff(f(x), x, x), g(t)) = 0, `&x`((diff(f(x), x))^2, g(t))-`&x`((diff(g(t), t))^2, f(x)) = 0}, [f(x), g(t)])

I really like this package:

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

But it seems to be impossible to intersect more than two areas:

restart;
with(InequalityGraphics);
inequalityplot(sin(x)+4*y>8 and x*y>1 and y<3-x,x=-2..4,y=-2..4,feasiblepoints=30);

It works if I leave out a condition:

inequalityplot(sin(x)+4*y>8 and x*y>1,x=-2..4,y=-2..4,feasiblepoints=30);

1. How do I 
2. Write a maple program to calculate the probability that k cards in 3-dimensonal set game (out of 27 cards) contain a set for k=3…9.

      3  randomly generate 12 cards in 4-dimensonal set game, either find a set or report “no set”.

question1.22.mwMy system differential equations are sys:=[deq1,deq2]

and the linearization point is given by linp:=[...]

Then I use linearization command to get a linear model of the system

linmod:=Linearize(sys,[input],[output],linp)

then maple returns a state-space equation of the linearized system.

I want to get the system transfer function, my method is ,first,extract...

i have the list of figure that produced by bellow maple.

for i from 0 to 10 do
  y := i/2;
  f := 2*x + y;
  plot( f, x, color=black); od;

i want the value of y print corresspond to  each line  of figure. 

for Example, the first value of y =0 should be  printed  correspond to the first  line and so on.

 Thus, in dispaly plot, each line correspond to its value of y.

could please help me how to do this?

being stalled myself here why don't these work?

for i from 1 to 200 do
  if ithprime(i)=prevprime(500) then print i; end if;
end do;

and

solve(ithprime(i)=prevprime(500),i)

??

What would be a better way?  I could do

for i from 1 to 200 while ithprime(i)<500 do end do; i;
                       ...

I am trying to construct a 7x7 matrix T such that T[i,j]=k, when multiply(monoid[i],monoid[j])=monoid[k]

I've defined monoid to be a boolean matrix such that there is at least one 1 in each row and column.

I am just not sure how to get my input to return the matrix. I keep getting a 7x7 zero matrix.

This is what I have:

T[i,j]:= Matrix(7);for i to 7  

do for j from 1 by 1 to 7    

do T[i,j]:=(multiply(monoid[i...