I ‘d like to compute something as follows at certain r and z, but first I must compute the integral and then substituting the r and z values.

A(r, z) = int( f(r,z, teta), teta=0..0.78)

But the problem: f is a complicated function of the 3 variables (including BesselJ(0,..), sin & root square functions, ..!!) thus a simple int command...

I hope to know

what the meaning of the follwing message

illegal use of a former parameter

Hi, I'm trying to fit an equation:

starting_salary= (intercept+SEX*B0) +B1*GPA  +B2* METRICS

(metrics,sex are  dummy variables, either 0 or 1) to see how sex coefficient affects the intercept ( how it affects starting_salary)

(later on I  want to see how metrics affects starting_salary for either sex=1 and sex =0 .. anyways I disgress) 

STARTING_SALARY :=:  

GPA:=:      

METRICS: =;  

restart;
lucoef:=Array(0..4, (i->i^2+3));
j:=1;
lucoefnew:=sum (lucoef[r]*lucoef[j-r],r=0..j);

Error, bad index into Array

===========
What it wants from me?!

This is the Pde I am trying to solve, with the initial profile function at t=0 called Piinitial, and IC=initialconditions,boundary conditions,

Piinitial := piecewise(r <= 1, (1/9)*r^2+(1/12150)*(31*(8807640902+1170*sqrt(24475122685023))^(2/3)+109499254+22553*(8807640902+1170*sqrt(24475122685023))^(1/3))/(8807640902+1170*sqrt(24475122685023))^(1/3), r <= 2 and 1 < r, (4/405)*(8807640902+1170*sqrt(24475122685023))^(1/3)*r+(14128936/405)*r/(8807640902+1170*sqrt(24475122685023...

This is the eq I am trying to solve

pde1:= diff(Phi(r,t),t,t) = (-(diff(Phi(r, t), r, r))*r^2-2*(diff(Phi(r, t), r))*r+2*(diff(Phi(r, t), r))^2+4*(diff(Phi(r, t), r))*(diff(Phi(r, t), r, r))*r+2*(diff(Phi(r, t), t, r))^2*r^2)/(r*(-r+2*r*(diff(Phi(r, t), r, r))+4*(diff(Phi(r, t), r))));

with initial conditions..

z=1000;

IC:={Phi(r,0)=Piinitial,D[2](Phi)(r,0)=0,D[1](Phi)(0,t)=0,D[1](Phi)(z,t)=0}:

pdsolve runs and returns the module..

Hi there

I have, for a while, been trying to find a way to force maple to display the unit im inputting. 

Im studying to be an engineer, and for an example the unit i often would like displayed is kilonewtons per square metre (kN/m2). 

Once i hit enter, maple converts it to Pa, which isnt really usable for me ... 

The question is, is there any way to stop maple from doing the conversion in output......

Does anybody have such package that converts function in series of powseries package i.e. addable, substructuble and so on. And mayby with functionality of series itself i.e. remember asymptotic point for which O(x^n) was present as simple series returns. I think with such pakage limit in attach limit_test.mw

plot3d(x^2+y^2+2*x, x = 0 .. 3, y = sqrt(4-x^2) .. sqrt(9-x^2))

It will only plot as far as x = 2, except for one line.  Maybe because sqrt(4-x^2) is undefined for x > 2?

How can I plot over an annular region?

p := proc() a:=1; b:=2; a/b end proc;

subs(1=1.0, eval(p));

proc () local a, b; a := 1.0; b := 2; a^1.0/b end proc;

why it generates a^1.0? It's rather nonsens. How to avoid such behavior?

I have a procedure that gets two lists A and B, I want to restore the name of the entry lists as string. How can I do that?

f(A,B) [procedure]

e.g. a(MyList,LL)

I want to restore a:="MyList" and b="LL"

TIA

Hi,

I would like to animate the solution to

x'(t)=y(t)

y'(t)= -4x(t)

that passes through x(0)=2, y(0)=1 as it orbits (0,0). I need to include the phaseportrait as the background. Any ideas on how to accomplish this?

Thanks!

Hey guys,

Really struggling with computing an algorithm on maple for the Stable Marriage Problem? Would be so appreciative for any help. Thanks Claudio

Does anyone know how to create one in maple? please see the example attached..

Compute the coefficients of
> a:=`@@`(D,6)(u)^3; 
#and
> b:=`@@`(D,5)(u)*`@@`(D,6)(u)*`@@`(D,7)(u);
#in the differential polynomial
> c:=`@@`(D,18)(u^3);.

Helppppppp :) TIA