# Items tagged with studentstudent Tagged Items Feed

### How can error in fsolve be corrected?...

December 17 2014
1 2

Could anyone assist in rectifying this error ''Error, (in fsolve) {f[1], f[2], f[3], f[4], f[5], f[6], f[7], f[8], f[9], f[10], f[11], theta[11]} are in the equation, and are not solved for''. Here is the worksheet FDM_Revisit_1.mw

### how to plot inequality with maple 12...

December 14 2014
0 2

I'd like to plot the following inequalities:

sqrt(x)<=1/sqrt(2)

1/sqrt(2)<sqrt(x)<=1/sqrt(2)

### can maple give answer for int(cos(x)^n,x) ?...

December 13 2014
2 11

Was trying to see if I can get the reduction formulas for int(cos(x)^n,x) in maple. But it seems no assumption used can make Maple give any result for this.  Mathematica gives a result using Hypergeometric2F1 (even with no assumption on n, which I am not sure about now), but was wondering why maple can't do this one:

restart;
int( (cos(x))^n,x) assuming n::integer;

int( (cos(x))^n,x) assuming n::posint;
same

In Mathematica, I get:

I am newbie in Maple, so may be I am missing some command or doing something wrong.

ps. I was trying to obtain

But this is lost case now. I just need to find out first why int(cos(x)^n,x) does not evaluate to anything in Maple.

fyi, the Hypergeometric result for $\int cos^n(x) \,dx$ can be seen in this reference (half way down the page):

http://www.integraltec.com/math/math.php?f=cosPower.html#cos

ps. can't one enter Latex in this forum like at stack exchange?

### Error, (in isolate) cannot isolate for a function ...

December 10 2014
0 9

Hello,

Im solving 4 ODE equations with BC. im trying to shoot the initial value but im having this error:

""Error, (in isolate) cannot isolate for a function when it appears with different arguments""

anyone could help me???

 >
 >
 >
 >
 >
 >
 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 >
 Error, (in isolate) cannot isolate for a function when it appears with different arguments
 >

Download shooting92.mw

### How to plot the Taylor Polynomial approximation...

December 08 2014
0 2

Hi everyone, I have been trying to plot the Taylor Polynomial approximation with the following code. However, my maple crushes everytime I run it. I indexed some of the variables to get the plot. The code works fine without the index. What did I do wrong?

y := array(1 .. 2);

Digits := 10;

n := 30;

h := .1;

T := 0;

X := 1;

f := (x, t) -> 1/(3*x(t)-t-2);

one := 1/(3*x(t)-t-2);

two := diff(f(x, t), t);

first := diff(x(t), t)

for k to n do

y[1] := subs(t = T(k), x(T(k)) = X(k), one);

y[2] := subs(first = y[1], t = T, x(T(k)) = X(k), two);

X[k+1] := X+sum(y[i]*h^i/factorial(i), i = 1 .. 2);

T[k+1] := T+h

end do;

X[n];

data := [seq([T[n], X[n]], n = 0 .. 30)];

p[2] := plot(data, style = point, color = blue);

p[3] := plot(data, style = line, color = blue);
display(p[2],  p[3])

The code without Index (which works fine)

y := array(1 .. 2);

Digits := 10;

n := 30;

h := .1;

T := 1;

X := .1547196278;

f := (x, t) -> 1/(3*x(t)-t-2);

one := x(t)^4*e^t-(1/3)*x(t);

two := diff(f(x, t), t);

first := diff(x(t), t);

for k to n do

y[1] := subs(t = T, x(T) = X, one);

y[2] := subs(first = y[1], t = T, x(T) = X, two);

X := X+sum(y[i]*h^i/factorial(i), i = 1 .. 2);

T := T+h

end do

### BodePlot function doesn't like inputs with units...

December 08 2014
0 5

I have been having problems with using the BodePlot function with units:

R1 := 18.2*10^3*Unit('Omega');

R2 := 10^3*Unit('Omega');

C1 := 470*10^(-12)*Unit('F');

C2 := 4.7*10^(-9)*Unit('F');

# wo is in hertz

wo := 1/sqrt(R1*R2*C1*C2);

# Q is unitless

Q := wo*R1*R2*C2/(R1+R2)

with(DynamicSystems);

sys := TransferFunction(wo^2/(s^2+wo*s/Q+wo^2));

This is the error message I got:

Error, (in Units:-Standard:-+) the units 1 and Hz have incompatible dimensions

I think the problem is that the BodePlot function doesn't expect 'wo' to have units.

So I tried to work around the issue by using the loglogplot but it doesn't seem to like

complex function even when I used abs to find the magnitude (with or without units).

Any workaround is appreciated.

### cluster analysis sing maple syntax... problem...

December 07 2014
1 4

Strings := {"buffalo", "zebra", "aardvark", "aardwolf", "anteater", "antelope", "bumblebee"};
{"buffalo", "zebra", "aardvark", "aardwolf", "anteater",

"antelope", "bumblebee"}
CP := {seq(seq(StringTools:-Take(x, StringTools:-CommonPrefix(x, y)), x in Strings), y in Strings)};
P := table(map(proc (x) options operator, arrow; x = "" end proc, Strings));
for x in CP do for y in minus(CP, {x}) do if and(StringTools*IsPrefix(y, x), implies(assigned(P[x]), StringTools*IsPrefix(P[x], y))) then P[x] := y end if end do end do;
Error, cannot determine if this expression is true or false: StringTools*IsPrefix("a", "")

Hi I tried to apply these command in maple to see hoe cluster graphic after that can work .. but i face this problem ,

I am a begginer in maple working on it for a term paper to demonstrate it is ability to view the cluster analysis .

if u have any benefcial resources I ll be so happy.

thanks very much

### modification of expression...

December 07 2014
0 2

Hello people in mapleprimes,

I want to ask a question about modification of expression.

Basically, what I want to do is to change the expression of sqrt((-a)^2) to sqrt(a^2).

The expression I want to modify is this:

aa:=phi[n, j] = a[j, D]-a-sqrt((-a[j, D]+a+tau[n, j])^2+4*gamma*f[c, n, j]/L[j])

I want to change this to

phi[n, j] = a[j, D]-a-sqrt((a[j, D]-a-tau[n, j])^2+4*gamma*f[c, n, j]/L[j])

To do this, I wrote as

subs(sqrt((-a[j, D]+a+tau[n, j])^2+4*gamma*f[c, n, j]/L[j]) = sqrt((a[j, D]-a-tau[n, j])^2+4*gamma*f[c, n, j]/L[j]), phi[n, j] = a[j, D]-a-((-a[j, D]+a+tau[n, j])^2+4*gamma*f[c, n, j]/L[j])^(1/2))

Or, I wrote as

subs(op([2,3,2,1,1],aa)=(a[j, D]-a-tau[n, j])^2,aa)

If there is better ways, please tell me them.

Best wishes.

taro

### I want to know some constant in Rutherford scatter...

December 02 2014
1 2

I don't know what mean is T and V0x?

Please.

### Quartic equation Solution in Maple...

December 01 2014
0 2

Dear all,

I want to know how can I solve a fourth or higher order equation to find out its roots, actually coefficients of each power terms itself are functins of other parameters not constant.

Please reply as soon as possible, as this is very urgent,

Regards.

### How to simulate elastic collision on maple?...

December 01 2014
2 0

As title , there has 2 balls(A and B) which I can change Quality 、Radius and Position for two dimensions(x and y axis),

A has a velocity and B quiescence , let A strike to B .

As shown below

I did a sample ,but i can't make it accurately .

test.mw

I want to let my sample like the below sample

http://www.mapleprimes.com/questions/202467-How-Does-Maple-Make-A-Rutherford-Scattering-

### How to point on specific summand in expression?...

November 29 2014
1 2

Hello,

Is it possible in Maple to select summand with some command?

For example:

eq := a+b = 0;
a + b = 0
eq-b;
a = -b

How to generalize second line?(in this case I explicitly told to subtract b, but I want to pick b by some command)

And I am not talking about solving, isolation or anything like that in this case. I just want to transfer something from lhs to rhs by selecting it beforehand.

Thank you.

### How to solve a nonlinear system of equations, if I...

November 27 2014
1 1

Hello,

I would like to solve a system of 9 nonlinear equations, with the constraints on all 9 variables to be that they are nonnegative. How can I do this?

My code is below - I am trying NLPSolve and have tried solve, but am getting stuck.

with(Optimization);

restart; eq1 := 531062-S/(70*365)-(.187*(1/365))*(H+C+C1+C2)*S/N = 0;eq2 := (4/365*(T+C))*S/N-(.187*(1/365))*(H+C+C1+C2)*T/N-(1/(70*365)+1/(5*365))*T = 0; eq3 := (.187*(1/365))*(H+C+C1+C2)*S/N-(4/365)(T+C)*H/N-(1/(70*365)+1/(4*365))*H = 0; eq4 := (.187*(1/365))*(H+C+C1+C2)*T/N+(4/365*(T+C))*H/N-(1/(70*365)+3/(8*365)+.2*(1/365)+.1)*C = 0; eq5 := .1*C-(1/(70*365)+1/(4*365)+1/60+.5)*C1 = 0; eq6 := (1/60)*C1-(1/(70*365)+1/(4*365)+1/210+.5)*C2 = 0; eq7 := .5*C1-(1/(70*365)+1/60+0.1e-2)*CT1 = 0; eq8 := .5*C2-(1/(70*365)+1/210+(1/9)*(0.1e-2*7))*CT2+(1/60)*CT1 = 0; eq9 := N-S-T-H-C-C1-C2-CT1-CT2 = 0; soln := NLPSolve({eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {C, C1, C2, CT1, CT2, H, N, S, T}, assume = nonnegative);

### Number Theory-FactorSet...

November 22 2014
1 1

Nm= p1. p2 ...pm + 1, for m more than or equals 1.

So N1 = p1 + 1 = 2 + 1 = 3, N2 = p1 p2 + 1 = 2  3 + 1 = 7, etc.

We prove that Nm is not divisible by any of p1, p2, . . . , pm, so that Nm is either a prime or it is divisible by a prime larger than pm.

(c) Use Maple to find out which of these numbers Nm, for m = 1, 2, . . . , 15, is actually prime.

Use Maple to compare pm with the smallest prime number that divides Nm, for m =1, 2, . . . , 15.

### How to simulate "Inelastic collision" on maple?...

November 22 2014
0 0

I want to simulate Inelastic collision

There has 2 ball which I can change  Quality and Radius.

one ball  move to another stirless ball with diferent angle

 1 2 3 4 5 6 7 Last Page 1 of 134
﻿