genDE := Diff(g(t),t)^2 - 1 - x*(g(t))^2 - (y-x^2)/12*(g(t))^4 = 0;gensol := dsolve(genDE,g(t));genfun := rhs(gensol[1]);

from above, i got generating function and then continue to derive got the following error

u = ln(-sqrt((-z^2+y)*(-3*z+sqrt(12*z^2-3*y)))*z^2/(sqrt(-z^3*(-3+2*sqrt(3)*csgn(z)))*(-z^2+y)))/x;

tau := subs(u=z,solve(q,z));

it said solution may have lost

if this can not be direct to...

restart;with(Statistics):

MartinPoisson := (p/((p+2*q*T)/(1+p*T+q*T^2)-(p*T+q*T^2)*(p+2*q*T)/(1+p*T+q*T^2)^2))^r*r^2*T^n/factorial(n);

Dist := subs(T=t,MartinPoisson) assuming t > 0;

MartinDist := Distribution(CDF = unapply(piecewise(t>1,Dist,0),t));

MartinDist := Distribution(CDF = unapply(Dist, t)) assuming t > 0; #Change 2

X:=RandomVariable(MartinDist);

MartinDensity := PDF(X,t);

MartinDensity := subs(t=x, MartinDensity);

SURF is the surface z=2x^2+y^2 for x=0 .. 2, y=0 ..4 and view= -1 ..16.

L is a straight line laser beam x=7-2t, y=t, z=4t-1 which hits SURF at the point P=(1, 3, 11).

Graph SURF, the laser path in red (t=0 to 3), the normal vector at P in black, and the reflected path in green.

Include the equations of the normal vector and the reflected path.

The surface z=f(x,y)=5/(1+x^2+y^2) is a hill. A bug walks along the path with x(t)=2+3cos(t) and y(t)=-1+2sin(t)

and z(t) is on the surface (0≤t≤2pi).

Plot the surface z for -4≤x≤7, -4≤y≤4 and includethe path (in thick black) of the bug on the surface.

(this will require 2 graphing commands plus a display 3d. Think spacecurve for the bug's path.)

plot part of the surface f(x,y)=10-x^2-3y^2:

plot3d(10-x^2-3y^2, x=-1 .. 4, y=1 .. 2, view=0 .. 11, axes= normal); (press enter)

Determine the normal vector and the equation of the tangent plane at the point (2,1,3) on the surface.

Display the surface , a black vector and the tangent plane, and give the equation for each of them

label the normal vector and the tangent plane, and give the equation for each of them

LL_61)_Poker.mw

I am having some problem with this procedure located in the "experiment" sub section.
I am trying to calculate the win probability in a poker game pre flop ie each player
get two cards. I am trying to find out which two card combination will result in the
highest probability of winning in the end.

In the first example I simply assume that player 1 gets [DA, HA] ie dimond ace

how to find eigen value this 4x4 matrix:

i copy from notepad, it do not have line return, so i copy to MS WORD and copy here from MS WORD

but some characters missing, if there is some missing please tell me, it seems old web page is better

code and error here

https://skydrive.live.com/redir?resid=E0ED7271C68BE47C!351

restart;with(LinearAlgebra):

f:=[x3-x2^3,-x2,x1^2-x3]:

with(Physics):

u := subs(x3=x3(t),subs(x2=x2(t),subs(x1=x1(t),u)));u := -Physics[`*`](Physics[`*`](Physics[`*`](3, Physics[`^`](conjugate(x2(t)), 2)), x2(t))+Physics[`^`](x1(t), 2)-x3(t), Physics[`^`](1+Physics[`*`](3, Physics[`^`](conjugate(x2(t)), 2)), -1));

u = subs(conjugate(x3(t))=x3(t),subs(conjugate(x2(t))=x2(t),subs(conjugate(x1(t))=x1(t),u)));

would like to remove all conjugate back to normal variable

have a matrix An = [alpha(i,j)]
alpha(i,j) = (j-i+1)*(j-1)!/(i-1)!*a(j-i+1)*x    for j-i >= 0, j - i even

                   0, j-i >= 0, j-i odd,           

                   -1, i-j = 1,             

                    0, otherwise

DOES ANYONES HAS THE SLIGHTEST IDEA OF HOW I COULD MODEL AND SUMILATE THE WAVE PATTERN OF A REGISTERED EARTHQUAKE??? FOR INSTANCE I AM STUDYING THE EARTHQUAKE EFFECTS OF THE NOVEMBER QUAKE IN GUATEMALA OF LAST YEAR. I WOULD LIKE TO MODEL THE MOVEMENT OF THAT QUAKE OVER A GIVEN AREA NEAR THE EPICENTER.

ANY HINTS OF HOW I COULD EVEN START TO MAKE THE MODEL AND SIMULATE??

Working.mw

This one is wokring fine.

However,

See the new worksheet: testP.mw

with the source code: tmP.c

######################################################################...

Hello, 

I would like to solve this differential equation: Eq2 := diff(z(x), x, x) = p*(1+(diff(z(x), x))^2)^.5 

I do not want an analytical solution, but I would like to obtain solutions on a graph, and I want to be able to change the value of the parametre p according to the value of x. 

For example, for x from 1 to 10 I want p=2, then for x from 10 to 11 I want p=10, then for x from 11 to 20 p=2...... 

Does anybody have any ideas please? :) 

Hi there, i want to write a procedure which can initiat the sexond procedure accodingly. For example, f1(a,b,c,d,e,f) where a,b,c,d,e,f is set of different atoms. say, Here I only need to deal with a and d type atoms to calculate energy values through f2 procedure.    

f1:=proc(a,b,c,d,e,f)   
.....
.....
end proc;

How i can write a procedure to assign values, like f1(yes,no,no,yes,no,no) or with number as f1(1,0,0,1,0,0) to show i only need a and d results. 

I try to make the inverse of a one to one piecewise function, but the solve command gives a function of lists. How can I avoid this?

restart;
f := x -> piecewise( x<=1, (x-1)/2, x<2, -3*(x-1)+3, 2*(x-2)+3 ):
solve( f(x)=y, x ): g := unapply(%,y): lprint(g(x));
piecewise(x <= 0, [1+2*x], x < 3, [2-(1/3)*x], 3 <= x, [1/2+(1/2)*x])