Hi, I've been introudeced to maple today and we were asked to think about how we'd write a two argument Ackermann code. I appreciate that computing power can lead to problems when workin gout this kind of recursive problem. I've done a bit of reading around but no joy on the maple front to be honest.

I thought the code may look like something like this but it's not right. Any thelp would be greatly appreciated!

# Ackermann-Peter function

I Follow the instructions to connect with matlab but when i type for example ( Matlab[setvar]("x",3.14);) it doesn't work.

i don't know why? my windows is 7 .

In Maple 14, some functions involving coth do not plot correctly when using "plot".

Example:

> restart;
rS := 10000; cS := 1; 
ZWarbSolid := rS*coth((((1/4)*rS*cS*I)*exp(2*Pi*freq))^.5)/(((1/4)*rS*cS*I)*exp(2*Pi*freq))^.5;
> with(plots); 
p1 := plot([Re(ZWarbSolid), -Im(ZWarbSolid), freq = .4 .. 3.4], color = red, scaling = constrained, style = point);p2 := plot([freq, -Im(ZWarbSolid), freq = .4 .. 3.4], color = blue, style = point...

s:=1/(y+a)^2;

s2:=1/(y+a)^(2/3);

type(s, smth rational)  should return true and type(s2, smth rational) -false. How to do that?

Hello,

I am doing some complicated computations in Maple, and I noticed something was wrong. For some reason, Maple makes the sum of 1 a zero, although he knows that it is a 1.

this is the code:

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?