1545 Reputation

19 years, 195 days

Complete the problem first......

First look at this:

```(() -> solve({args}, indets({args}, name) minus {t}))(eq || (`\$`(1 .. 4))); subs(t = 0, %); allvalues(%);```

What version of maple do You use?...

Could it be that you need a registered version of Maple?

Floor....

`plot(floor(x/3)*3, x=-5..5);`

You can see the same in all versions of ...

Simply compare these two codes:

```restart;
P1 := 1007;
P2 := 1014;
P3 := 1014.1;
evalf(P2 - P1, 2);
evalf(P3 - P1, 2);
```
```restart;
P1 := 1007.;
P2 := 1014.;
P3 := 1014.1;
evalf(P2 - P1, 2);
evalf(P3 - P1, 2);
```

The difference is obvious. No bugs here.

Thumb if You like.

BandMatrix....

For Your matrix you can use this:

```F := (n, f) -> LinearAlgebra:-BandMatrix(map(f, [\$1 .. n]), 0, n);
F(4, (n) -> d[n]);
```

Thumb if You like.

Valuation rules....

`for i from 1 to n do proc(i, a, b) a[i]:=b[i] end(i, a ,b) od;`

also

`for i from 1 to n do proc(i) a[i]:=b[i] end(i) od;`

provided a, b are not declared local.

Cancellation....

fgP,fgQ:=unapply(gP,x,y), unapply(gQ,x,y);
(fgP(QS1,QS2)*fgQ(S1,-S2)/(fgP(PS1,PS2)*fgQ(S1,-S2))) mod p;

Notice the cancellation.

Frontend may be enough for the task....

```e_n_1b := n_1 = (-w^(2*sigma)*tau + w^sigma)*s*nsp_1/(w^sigma*tau^2 - w^(2*sigma)*tau + w^sigma - tau) + tau*(w^(sigma - 1)*tau - w^(2*sigma - 1))*s*nsp_2/(w^sigma*tau^2 - w^(2*sigma)*tau + w^sigma - tau);

wpowers := (A, e) -> frontend(convert, [subsindets(A, dependent(w)^anything, expand), parfrac, e], [{Non}(identical(w^sigma)), {}]);
wpowers(e_n_1b, w);
wpowers(e_n_1b, w^sigma);
```

Thumb if You like.

Function parameters....

Edit: There was an error in my response, just like there was one in your code.
Correction:

```PlotGraph := proc(func::anything, opts::(seq({identical(zoom) = [range, range]})))
plots:-display(plot(func), ifelse(membertype(identical(zoom) = [range, range], [opts]), view = subs([opts], zoom), 'NULL'))
end proc;;
PlotGraph(x^2, zoom = [0 .. 1, -5 .. 5]);
```

Thumb if You like.

Solve twice....

One variable, constant bounds:

```a := solve({-1 < x, 0 < (2*x)/(x^2 - 1), x < 1}, {x});
solve(`and`(op(a)));
```

Thumb if You like.

Nullify unwanted elements....

See this example:

```A := map((x -> ifelse( x<1/3 or x>2/3, x, 0)), LinearAlgebra:-RandomMatrix(25, 25, generator = rand(0 .. 1.0))):
P := plots[matrixplot](A, heights = histogram, transparency = 0.3, style = patchnogrid, shading = zhue, orientation=[30, 30, 0]):
subsindets(P, [[anything, anything, 0.]\$4], (x -> NULL));
```

Thumb if You like.

cosh(Q)....

The point at K=4*I..1+5*I that appears a solution is at:

`{K = (1/2)*sqrt(nu*Pi*(-9*Pi*nu+(8*I)*h^2))/nu};`

The above is a solution to cosh(Q)=0. cosh(Q) appears in the denominator.

Thumb if You like.

Use algsubs....

```A := (4*x + 1)/(diff(f(y), y)) = a;
b := (3*(4*x + 1))/((diff(f(y), y))*(3*x + 1));
frontend(algsubs, [A, b, diff(f(y), y)], [{Non}(function), {}]);
```

Thumb if You like.

Numer....

Do You mean to do something like this:
"numer(L = 2/3)";
```f := e -> lcm((denom@lhs, denom@rhs)(e))*e; f(L = 2/3);```

Factor....

`collect(Your_stuff_goes_here, s, factor);` #normal might do also

Thumb if You like.

 5 6 7 8 9 10 11 Last Page 7 of 21
﻿