I have looked at the maple help for the ErrorPlot function, but i can't make it work for my plot. 

I have the following procedure:

K := 1.; Xrng := [seq(dt*i, i = 0 .. ntime)];

value := array(1 .. ntime+1);

ave := array(1 .. ntime+1);

for j to ntime+1 do:

discountfactor := exp(-r*j*dt);

Stot := 0.;

Stot1 := 0.;

for i to npath do:

Spath := max(Spot[i][j]-K, 0);

Stot := Stot+Spath;

When I use pdsolve to solve a system of PDEs, and tell it to use HINT='*', does the solution I get ONLY contain those obtainable by separation of variables, or does it use separation of variables as a starting point and deliver the general solution? And how can I tell?

Dear all,

I wrote some simple functions in Maple 13, called "myproc" file,

and used startup code (edit-->startup code) to initialize it in a worksheet by: read "myproc";

Now I am using Maple 14, it does not work, says, "unable to read", sort of thing.

I have also tried to rename the file to myproc.m and read "myproc.m";

does not work either.

Is there any compatibility problem?

Thanks.

casper

I was wondering if there's a simple way to add something like an interpolation loop to a procedure. For example, say I have a function of 2 vars which is expensive to compute, f(x,y). I can use option cache to speed up points which are already computed, but is there an easy way to interpolate between them using ArrayInterpolation? That is, if f varies on scales of order 1, if I have calculated it at 0.1, 0.2, then It can be guessed with a spline at 0.15 accurately enough. What...

I have the following procedure that generates Sierpinski's triangle:

seirpinski := proc (n)
> local Mat, Vector1, Vector2, Vector3,
> Prob1, Prob2, Q, prob, counter, fractalplot;
> Mat := linalg[matrix]([[.5, 0.], [0., .5]]);
> Vector1 := linalg[vector]([0, 0]);
> Vector2 := linalg[vector]([0, 0.5]);
> Vector3 := linalg[vector]([0.5, 0]);
> Prob1 := .333; Prob2 := .333;
> Q := linalg[vector]([0, 0]);
> writedata("fractaldata", [[Q[1...

I guess someone is kidding with me...

I want unsubscribing of this s... place!

I have this non-linear ODE in q(p):

ode := diff(q(p), p) = (q(p)-a*p)/(p-a*q(p))*(p^(-b)-1)/q(p)^(-b);

for some real parameters a and b, it may be solved in exact form for the special value a=1.

Is there any hope of getting a closed-form solution or of expressing the solution in a neat way for the case 0<a<1? even some specific value of a would be of interest, e.g. a=1/2.

dsolve(eval(ode,{a=1}));

let us consider the series

y = a_0 - sum(a_i,i=1..n)

here the terms a_i are defined recurrsively as follows

a_i = (a_(i-1) -4*ln(x)/ln(10) + 4/10 + 4*ln(a_(i-1)) )/ ( 1+ 4/(a_(i-1) * ln(10)))

how can i program this in maple. so i can...

what is the general sum formula for the following seires (we may notice from these few terms it has a pattern)

f:=x->x^3:
y[1]:=0.5:
y[3]:=subs(x=y[1],f(x)):
i := 2:
y[i] := y[i-1]-(y[1]-y[2+1])*(1/2):
i:=1:
g[i] := (y[i+1]-y[i+2])*x/(y[i]-y[i+1])+(y[i+2]*y[i]-y[i+1]*y[i+1])/(y[i]-y[i+1]):
g[i] := unapply(g[i], x):
h[i] := unapply(solve(y = g[i](x), x), y):
for i from 2 by 1 to 7 do
g[i]:=(f@h[i-1])(x);
g[i]:=unapply(g[i],x);
y[i+2]:=subs(x=y[i+1],g[i](x));
s:=solve(z=g[i](x),x);
h[i]:=unapply(s[1],z);
end do:

for i from -1 by -1 to -8 do

 Hello,

I am trying to use a worksheet developed by Professor Batista [ritz06.mws] for calculating plate vibration frequencies using the Ritz method, which is available from the Applications directory - Mechanical Engineering folder.  The worksheet was developed using Maple 12 and I am using Maple 11; a careful reading of the file didn't show any commands that were incompatible across the two versions and yet I am not able to solve the worksheet using the default inputs provided in the worksheet.  I have attached a PDF output of the file with the error message.  Any assistance that could be offered would be greatly appreciated.

Please send your response to wayne.bell@fpl.com.

Kind regards,

Symbolic_integral.mw

Thanks in advance

Friends,

I am very beginner of Maple softs, in order to proceed with my project i need to model a simple Human Knee joint and  to do its simulation, Somebody please help on this, otherwise i cant able to proceed with my project.

Waiting for a solution

Vinil