Can anyone help me with the question? Keep getting the error. I don't know what is going on. Thank you in advance.

QUESTION_6.docx
 

Download Question_6_maple.mw

Hello all

I need to calculate the distance L in a curvefitting. The calculation has been made in MathCAD and I need to copy it in Maple, but I cant seem to make it solve.

Here are the two pictures of my problem:

In MathCAD:

In Maple:

Here is the copy of Maple:

V[Ed13] := [257.184, 230.4, 184.3, 138.2, 92.2, 46.1, 0];
 = 
    V[Ed13] := [257.184, 230.4, 184.3, 138.2, 92.2, 46.1, 0]

x[12];
 = 
            [-149, 208, 567, 925, 1283, 1642, 2000]

L[1] := x[13];
 = 
         L[1] := [0, 358, 717, 1075, 1433, 1792, 2150]

solve(V[Rdc] = CurveFitting:-ArrayInterpolation(x[12], V[Ed13], L[1]), L[1]);
 = 
Error, invalid input: solve expects its 1st argument, eqs, to be of type {`and`, `not`, `or`, algebraic, relation(algebraic), ({list, set})({`and`, `not`, `or`, algebraic, relation(algebraic)})}, but received 153.25 = [HFloat(246.00524369747902), HFloat(211.1381615598886), HFloat(164.98435754189944), HFloat(118.92625698324022), HFloat(72.93816155988858), HFloat(26.784357541899443), HFloat(-19.315642458100548)]

Find 
L[1];

Answer should be 808
 

Hope someone can help as I honestly have no idea how to do this :-)

I want to calculate the below triple integral numerically. I tried different methods and reduced accuracy but nothing works.

Anyone capable?

Thanks

restart;
Digits := 15;
with(VectorCalculus);
v1 := <0, 0, 1>;
v2 := <sin(theta2), 0, cos(theta2)>;
v3 := <VectorCalculus:-`*`(sin(theta3), cos(phi3)), VectorCalculus:-`*`(sin(theta3), sin(phi3)), cos(theta3)>;
v1v2 := CrossProduct(v1, v2);
v2v3 := CrossProduct(v3, v2);
DotProduct(v1v2, v2v3);
(simplify(VectorCalculus:-`*`(%, VectorCalculus:-`*`(Norm(v1v2), Norm(v2v3))^VectorCalculus:-`-`(1))) assuming (0 < theta2, theta2 < Pi, 0 < theta3, theta3 < Pi));
evalf(VectorCalculus:-`*`(Int(VectorCalculus:-`*`(arccos(%), sin(theta3)), [theta2 = 0 .. Pi, theta3 = 0 .. Pi, phi3 = 0 .. VectorCalculus:-`*`(2, Pi)], epsilon = 0.001, method = _d01akc), VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, Pi), 2), Pi), 4), Pi)^VectorCalculus:-`-`(1)));

Hey there

I have noticed that after some time, Maple freezes and says in the bottom left corner that it is "evaluting" when I  copy something. This is becoming unbearable, as it makes every assignment take a lot longer. Why is this happening, and how can I prevent it?

Hi

I need your help with this problem. The date protocol here in the antipodes is day/month/year.

I get this error...Maybe there is an alternative.

does_not_evaluate.mw

Hello, I have a question about converting an output of a expression into another.

solve({-1 < x, 0 < 2*x/(x^2 - 1), x < 1}, {x})

This is my expression which evaluates into: . Which is correct but I would like it to be a classic range, more precisely (-1,0).

Is there a command to do that?

Thanks,

David.

I compared the numerical value of an integral expression in Gradshteyn and Rhyzhik (8.431.5) for the BesselI(v,z) function to the BesselI(v,z) function in Maple 2019.2 and got different answers.  The GR expression is the difference of two integrals and it looks like Maple only evaluates the first integral, which is correct only if the order v is an integer.  Please see the attached.

BesselI.mw

Hello,

I am trying to pass the graphics to .txt format so that I can plot them in Origin (OriginLab).

npontos := 10;
dt := 0.02;
fd := fopen("out.txt", WRITE);
for i from 0 to npontos do
    t[i] := i*dt;
    PA[i] := x[i](t);
    printf("%e %e\n", t[i], PA[i]);
    fprintf(fd, "%e %e\n", t[i], PA[i]);
end do;
0.000000e+00 

Error, (in fprintf) number expected for floating point format
fclose(fd);
npontos := 10;
dt := 0.02;
fd := fopen("out.txt", WRITE);
for i from 0 to npontos do
    t[i] := i*dt;
    PB[i] := x[i + n](t);
    printf("%e %e\n", t[i], PB[i]);
    fprintf(fd, "%e %e\n", t[i], PB[i]);
end do;
fclose(fd);
0.000000e+00 

Error, (in fprintf) number expected for floating point format
npontos := 10;
dt := 0.02;
fd := fopen("out.txt", WRITE);
for i from 0 to npontos do
    t[i] := i*dt;
    qA[i] := x[i + 2*n](t);
    printf("%e %e\n", t[i], qA[i]);
    fprintf(fd, "%e %e\n", t[i], qA[i]);
end do;
fclose(fd);
0.000000e+00 

Error, (in fprintf) number expected for floating point format
npontos := 10;
dt := 0.02;
fd := fopen("out.txt", WRITE);
for i from 0 to npontos do
    t[i] := i*dt;
    qB[i] := x[i + 3*n](t);
    printf("%e %e\n", t[i], qB[i]);
    fprintf(fd, "%e %e\n", t[i], qB[i]);
end do;
fclose(fd);
0.000000e+00 

Error, (in fprintf) number expected for floating point format
npontos := 10;
dt := 0.02;
fd := fopen("out.txt", WRITE);
for i from 0 to npontos do
    t[i] := i*dt;
    Tg[i] := x[i + 4*n](t);
    printf("%e %e\n", t[i], Tg[i]);
    fprintf(fd, "%e %e\n", t[i], Tg[i]);
end do;
fclose(fd);
0.000000e+00 

Error, (in fprintf) number expected for floating point format
npontos := 10;
dt := 0.02;
fd := fopen("out.txt", WRITE);
for i from 0 to npontos do
    t[i] := i*dt;
    Ts[i] := x[i + 5*n](t);
    printf("%e %e\n", t[i], Ts[i]);
    fprintf(fd, "%e %e\n", t[i], Ts[i]);
end do;
fclose(fd);
0.000000e+00 

Error, (in fprintf) number expected for floating point format

trabalho_final_2019.mw

I would like to know how to create a group with some elements (numbers) {1,2,4,...} with a defined operation on it like X_42 (modular  multiplication by 42).

What is the command to yield the expansion of Z^N in terms of x and y.

Z complex = x+iy,

N integer >1

e.g. Z^2 = x^2 - y^2 + 2.i.x.y

could you just delete the post !!!!!!!!!!!!!!!!!!!!!!!!!!!x!xxxx!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Hi,

Let's take the last example (Maple 2019) given in help[sparsematrixplot] (representation of the adjacency matrix of a graph).
Vertices of this graph are labelled 1, 2, ...20.
Suppose I change these names as a, b, ...t.
I would like the tickmarks of the sparsematrixplot output match the names of the vertices of the graph, and not the integers 1, 2, ..20

I tried this:
S := [$1..20] =~ StringTools:-Char~(96 +~  [$1..20]);
plots:-sparsematrixplot(..., tickmarks=[S, S])

But the only the tickmars of the columns are changed, not those of the rows.

Is it possible to change the names of the tickmarks ?


Thanks in advance.

Hello,

Please, what is going wrong that it is not graphing the ODE system solution?

eqs := seq(eq[i], i = 1 .. 6*n);
cis := seq(ci[i], i = 1 .. 6*n);
sol := dsolve([eqs, cis], numeric, stiff = true, output = listprocedure);
for i to n do
    graf1[i] := odeplot(sol, [t, x[i](t)], 0 .. 5, color = black);
end do;
for i from 11 to 2*n do
    graf2[i] := odeplot(sol, [t, x[i](t)], 0 .. 5, color = blue);
end do;
for i from 21 to 3*n do
    graf3[i] := odeplot(sol, [t, x[i](t)], 0 .. 5, color = green);
end do;
for i from 31 to 4*n do
    graf4[i] := odeplot(sol, [t, x[i](t)], 0 .. 5, color = red);
end do;
for i from 41 to 5*n do
    graf5[i] := odeplot(sol, [t, x[i](t)], 0 .. 5, color = pink);
end do;
for i from 51 to 6*n do
    graf6[i] := odeplot(sol, [t, x[i](t)], 0 .. 5, color = orange);
end do;

display(seq(graf1[i], i = 1 .. n));
display(seq(graf2[i], i = 11 .. 2*n));
display(seq(graf3[i], i = 21 .. 3*n));
display(seq(graf4[i], i = 31 .. 4*n));
display(seq(graf5[i], i = 41 .. 5*n));
display(seq(graf6[i], i = 51 .. 6*n));

trabalho_final_2019.mw

BIOMAT_project.mw

I am looking at an SIR model and an adapted model that has a vaccine for measles. I would like to use Maple to plot the difference between the two models. 

All the work that's been done so far is in the attached file. I would like to have multiple graphs, changing the starting values but I just need the general idea of how to do this first. 

Hi

I'm trying to rearrange a transfer function to get my s^2 term with other specific terms . Below is the denominator of the transfer function. I know that my goal is to get (Izz*u)/m for the s^2 term, so I divide by Df*Dr*(a+b)*u*m, as shown below.

However, if we look at the s^2 terms for now, the equation should be further simplified leaving me with simply (Izz*u)/m, but Maple doesn't want to do this for some reason, and if I use simplify on the above expression it tries to expand everything again, and forgetting about my collect(fun, s).

Is there any way to get around this?

Many thanks