Bug in my Division by Zero Calculator...

Basically it spits out the subset of values for which a division by zero error will occur for the function you specify on  range you specify for each of it's arguments, but I get an ambigous error when ever exponentiation features in the function I specify, which of course dramatically reduces the application of the calculator. Division,addition,substraction and multiplication are currently the only available arithmetic operators availble for the function window that I know the error will not occur.

If some one can help it is much appriciated

DIVISION_BY_ZERO_CALCULATOR.mw

problem with piecewise and differentiate piecewise...

Dear all

I have an first-order ode

y'(x) =0, for x in [0,1] but x different to 1/2

and we define a function y(x)= 2   if 0<=x<0.5  and y(x)=0 for 0.5<=x <=1

The attached code define these function ( maybe piecewise function is not welle defined) and how can we show that y is solution or not of the differential equation

ode.mw

many thanks

Plotting depth more real on 3d plots...

I was playing around with some of the maple examples(specifically the knots) and all the plots seem difficult to visualize in 3D... well, more difficult than they should.

In some cases the plots look more like projections and 2D because the graphical representation does seem to project in to the depth of the screen in a natural way.

Mainly it happens when rotating when some of the curve's appear to stay in front due to the highlighting and size issues.

I'm not sure exactly what the problem is but I imagine it definitely can be improved. It may be simply that there are no other depth cues in the plot to make things more real.

1. I'd like anti-aliased edges. The edges of the plot are hard and not correct coloring wise(do not fade like in real life).

2. Some type of depth perspective where the size changes with depth in to the screen in a meaningful way.

3. Proper coloring of this. Maple seems to handle the lighting correctly but it is far too weak to give the impresson of depth

I think maple essentially does all this(except probably AA, although there is probably a setting somewhere) but it seems they are not properly tweaked for maximum realism. I played with the lighting and projection and it seems to help a little but not much.

Any ideas on how to improve this? Ideally something that doesn't have to be set every plot.

how I solve a Hard Triple Integral?...

hello everyone,
INGT.mw

 (1)

 (2)

 (3)

 (4)

 (5)

 (6)

 (7)

 (8)

 (9)

 (10)

 (11)

 (12)

 (13)

 (14)

 (15)

 (16)

 (1)

 (2)

 (3)

 (4)

 (5)

 (6)

 (7)

 (8)

 (9)

 (10)

 (11)

 (12)

 (13)

 (14)

 (15)

 (16)

I'm trying to calculate a triple integral complicated by a procedure that changes each time a variable xi, while the program takes a lot of time and it gives me the message "Warning, computation interrupted". If anyone can help me I will be very happy

Simplify command on Binomial expressions...

Is there something I should be doing whenever I use simplify to avoid things like this, or should I stop using the "is" function all together?

 >

 >
 (1)
 > #And we have:
 (2)
 > #And since:
 (3)
 >

"How do I solve Fredholm or Volterra integral equa...

Can I know how to write code about to solve the second kind of linear Fredholm and also Volterra approximation integral by Simpson's method and Galerkin method, and I need a book about advance math to the integral equation, thanks.

Some Questions about Maple Symbolic Toolbox in Mat...

1. Does it have any detail manual for Matlab programming using Maple symbolic engine? (Except help documentation of Maple Symbolic Toolbox in Matlab).
2. How to maximize the running speed of the codes programmed by Maple symbolic toolbox? Any Notes?
3. How to use the function 'GenerateMatrix' of Maple in Matlab? What's the syntax to use this function in Matlab?

Best Regards.

Si

How to make Maple to integrate term by term?...

Hello

I wonder how integrate can be applied term by term to the following nonlinear differential equation

`(diff(y(t), t))*(diff(y(t), t, t, t))-(diff(y(t), t))*y(t)^2-(diff(y(t), t))*(diff(y(t), t, t))-(diff(y(t), t))*A*y(t)`

The expected output will be something like

`(diff(y(t), t))*(diff(y(t), t, t))-(1/3)*y(t)^3-(1/2)*A*y(t)^2-(1/2)*(diff(y(t), t))^2+C+int((diff(y(s), s, s))^2, s = 0 .. t)`
`map(x-> integrate(x,t),(diff(y(t), t))*(diff(y(t), t, t, t))-(diff(y(t), t))*y(t)^2-(diff(y(t), t))*(diff(y(t), t, t))-(diff(y(t), t))*A*y(t))`

solves most of it but not the first part.

Many thanks

Ed

Problem solving log equation...

I was attempting to construct  a log question and noticed a problem with the answers that are being resolved by maple.

restart:
with(StringTools):
a := 2:
b := 10:
c := -4:
d := 5:
rs := 4:
#notice I use the word filler below, in order to randomize later and replace "filler" with various "log[base]" or "ln".
equation1 := convert(filler(a*x+b)-filler(c*x+d) = rs, string):
equationAct := parse(SubstituteAll(equation1, filler, cat(ln))):

#it appears maple is applying the log properties and suggesting solutions outside the real domain, for the individual log expressions.  The plot of the equation seems to confirm this.
plotA := plot(lhs(equationAct), x = -10 .. 10):
plotB := plot(rhs(equationAct), x = -10 .. 10):
plots[display](plotA, plotB);

Newton metho for system...

restart:
NewtonRaphsonSYS:=proc(X,F,X0,TOL,itmax)
local Xn,H,FX,J,i,m,n,Err;
m:=LinearAlgebra[Dimension](Vector(X));
Xn:=Vector(m);
H:=Vector(m);
FX:=Vector(m);
J:=Matrix(m,m);
Xn:=X0;

for n to itmax do
FX:=eval(F,[seq(X[i]=Xn[i],i=1..m)]);
J:=evalf(VectorCalculus[Jacobian](F,X=convert(Xn,list)));
H:=-Inverse(J).Vector(FX);
Xn:=Xn+H;

printf("%3d %.8f",n,Xn[1]);

for i from 2 to m do;
printf("%.8f",Xn[i]);

end do;
for i from 1 to m do;
printf("%.3g",H[i]);

end do;
printf("\n");
if (LinearAlgebra[VectorNorm](H,2)<TOL) then break end if;
end do;
end proc:
F:=[xy-z^2-1,xyz+y^2-x^2-2,e^x+z-e^y-3];
X:=[x,y,z];
X0:= <1,1,1>;
TOL:=0.00000001;
itmax:=10;

[  2             2    2             x        y    ]
[-z  + xy - 1, -x  + y  + xyz - 2, e  + z - e  - 3]
[x, y, z]
Vector[column](%id = 4442374402)
-8
1 10
10
NewtonRaphsonSYS(X,F,X0,TOL,itmax);

Maple is not able to simplify x^3+a*x+b) -(y^3+a*y...

Maple is not able to simplify x^3+a*x+b) -(y^3+a*y+b) )^2/(x-y)^4   into (x^2 + x*y + y^2 + a)^2/(x - y)^2

wtf

p.s.  how do I get maple to show steps in simplifications

Error, when trying to remove a task......

Hello people :)

As the captian says, im trying to remove an old task i've made.
But i get this:

Error in Get, invalid object [_XML_reply_data_get("reference" =
"_Maplets_reference_12","parameter" =
"#xa6","tning"))]

And i have no idea what it is, but it won't erase my task :'D

Have a great weekend you all
Best regards Lucas :)

How do i work with multivariate functions?...

Hi

I am a student of economics writing my thesis. I have an inverse demand function defined:

p[i,j]=1-q[i,j]-αq[i,k]-β(q[h,j]+αq[h,k])

where α,β are scalars, i{1,2,}, j {1,2,3}, h/=i and k/=j.

I wanted to know how do i input this kind of function in the program and how do i calculate the function of q[i,j] (the inverse)

Thank you!

How to extract x and y values from a list of order...

I would like to extract the x and y values from a list containing ordered pairs, [x,y], and display as a list.

For example, given

A:= [ [4,2], [9,1], [6,8] ]

how can I obtain

X:=[4,9,6] and Y:=[2,1,8]

using a simple command?

Thank you!

numtheory package lamba wtf...

Why is lamba protected as something it traditionally isnt for this package?

 >

 >
 (1)
 >
 >
 (2)
 >