f[2](x) = -(3*R^2*x^11)/30800 + R^2*x^9/420 - (177*R^2*x^7)/9800 + M1*M2*x^5/40 + (17*R^2*x^5)/700 + ((-(3*M1*M2)/10 - (443*R^2)/21560)*x^3)/6 + (M1*M2/40 - (137*R^2)/26950)*x

Hi,

I am writing the following code and MAPLE is giving me operator error. Please help (file: doubt_6.mw)

d2:=100000000

for m in set_m do
    for n in set_n do
        SOL1 := fsolve({ODE11, ODE12}, {N, t__2});
        N1:=eval(N,SOL1);
        t_2_1 :=eval(t__2,SOL1);
        T_1:= eval(T, [lambda = 3, a = 300, b = .15, c = .25, A__m = 300, A__d = 150, A__r = 50,C__m = 4, P__m = 8, P__d = 10, 
        P__r = 12, theta__m = .15, theta__d = .12, theta__r = 0.5e-1, h__m = .2, h__d = .3, h__r = .5, i__m = .1, i__d = .1, 
        i__r = .1, i__om = .1, i__OD = .15, i__c = .3, i__e = .2, M = 2, alpha = 0.2e-1, t__2 = t_2_1]):
        t_31:= T_1 /m ;
        t_41:= T_1 /(m*n) ;
        if (N1<=t_41 and 2>=t_31) then
            d1:= eval(TCS__1, [lambda = 3, a = 300, b = .15, c = .25, A__m = 300, A__d = 150, A__r = 50, C__m = 4, P__m = 8, 
            P__d = 10, P__r = 12, theta__m = .15, theta__d = .12, theta__r = 0.5e-1, h__m = .2, h__d = .3, h__r = .5, i__m = .1, 
            i__d = .1, i__r = .1, i__om = .1, i__OD = .15, i__c = .3, i__e = .2, M = 2, alpha = 0.2e-1, t__2=t_2_1, N=N1]):
            if (d1<= d2) then
                d2:= d1;
                print("Value is updated",d2,N1,t_2_1,"for",m,n)
            end if
        end if    
        print(N1,t_2_1,t_31,t_41,d1,m,n)
    end do
end do

Thanks in advance

Any idea why Maple can simplify this expression only when multiplied by -1 but not otherwise?

expr:=sqrt(-2 + 2*I*sqrt(3));
simplify(expr)

But now

simplify(-1*expr)

Multiplying by -1 should not have mattered. Becuase -2 + 2*I*sqrt(3)=(1+I*sqrt(3))^2  Therefore sqrt(-2 + 2*I*sqrt(3)) is the same as 1+I*sqrt(3) and Maple knows this:

simplify((1 + sqrt(3)*I)-expr) gives zero.

But Maple only simplifies -expr  and not +expr.

Just wondering is someone can suggest why that is and if this is considered OK?

Maple 2020.1, Physics 724

Hello

I need to check if the solution (sols) of a (nonlinear) polynomial system of equations (the coefficients are not numeric) has only one solution for, let's say, y and z (one solution for y and one solution for z).  I cannot use allvalues (not threadsafe according to CodeTools:-ThreadSafetyCheck) and then count the number of solutions. 

I have removed all solutions that fit the command 

ormap(x->x=true,map(has,rhs~(op~(sols)),_Z))

(RootOf works too).  (Please tell me if I am using ormap correctly).  

Even using the above command and then checking if nops(sols)=1, not all one solutions are caught.  

Many thanks

Ed

Any idea what this could mean?

I'm comparing 2 matrices where one of them is empty, and the other one is unitfree.

Hello,

I want to ask, I have a system of equations... 

eq1:=K_1=(((n_Cl-x)*u_Cl)*(n_H*u_H))/(n_HCl*m);
eq2:=K_2=(((n_Na-x)*u_Na)*(n_OH*u_OH))/(n_NaOH*m);
eq3:=K_w=(n_H*u_H/m)*(n_OH*u_OH/m);
eq4:=K_NaCl=((n_Na-x)*u_Na*(n_Cl-x)*u_Cl)/m^2;
eq5:=(n_NaCl-x)=(n_Na-x)+n_NaOH;
eq6:=(n_NaCl-x)=(n_Cl-x)+n_HCl;
eq7:=(n_Na-x)+n_H=(n_Cl-x)+n_OH;
eq8:=2*ionic=(n_H/m)+((n_Cl-x)/m)+((n_Na-x)/m)+(n_OH/m);
eq9:=u_H=0.4077*ionic^2-0.3152*ionic+0.9213;
eq10:=u_Na=0.0615*ionic^2-0.2196*ionic+0.8627;
eq11:=u_OH=0.1948*ionic^2-0.1803*ionic+0.8887;
eq12:=m=r*V;
eq13:=u_Cl=(1.417625986641341e-01)*exp^(-ionic/2.199955601666953e-02)+2.369460669647978e-01*exp^(-ionic/3.756472377688394e-01)+5.859738096037875e-01;

but when I use solve command, Maple starts computing, but suddenly gives message: Solutions may have been lost. 

I guess problem is the last exponential equation. Is there any way how to make Maple to be able to compute this system?

What is the best way to get some one on one training in maple , has anyone done this?

Area of shadow : i have studied one way of solution by hand, but now in Maple ?
Note: for math illustrations : geogebra seems to be wel suited

I would like to plot a specific vector with an initial point of <4,3,-5> and a terminal point of <3,-1,4>. I have been searching through Maple Primes and Maple Soft.

restart;
M1 := 12.3:
M2 := 12.4:
M3 := 12.5:
R := 50:

EQ:={(diff(F(x), x $ 4)) - M1*diff(G(x),x$2) -2*R*F(x)*diff(F(x),x$3)=0, diff(G(x),x$2)+ M2*(diff(F(x),x$2)-2*G(x)) +M3*(diff(F(x),x)*G(x)-2*F(x)*diff(G(x),x))=0}:

IC:={D(F)(-1)=0, D(F)(1)=0,F(-1)=-1,F(1)=1,G(-1)=0,  G(1)=0}:

sol:= dsolve(EQ union IC,numeric,output=Array([-1,-0.9,-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1]));
Error, (in dsolve/numeric/bvp) uanble to achieve continuous solution with requested accuracy of 0.1e-5 with maximum 128 point mesh (was able to get 0.47e-5), consider increasing `maxmesh` or using larger `abserr`
plots:-odeplot(sol,[x,F(x)],x=-1..1,color=red,axes=box)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution
plots:-odeplot(sol,  [x, diff(F(x), x)], x = -1 .. 1, color =  green, axes = box)
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution
plots:-odeplot(sol,[[x,F(x)],[x,diff(F(x),x)],[x,G(x)]],x=-1..1,color=[red,green,blue],axes=box):
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution
plots:-odeplot(sol, [x, G(x)], x = -1 .. 1, color = blue, axes = box);
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

solve( sum((1/(1+x+x^2))*y^n, n=1..infinity) = 1, parametric=true, real);

How to solve two variables summation equation with parametric true or real triangularize?

how to solve for x and y?

I just wonder how it is possible to store data between sessions. What I found out is that if a matrix is connected with a datatable, the matrix data is stored. If it is not connected with a datatable, it is not.

In the attached example there are 2 matrices. If you run the worksheet, both variables are defined as visible in the variable overview.

After reopening without running the worksheet, just the variable connected with a datatable is.

Any ideas how to also get the second variable stored between sessions (obviously not with the same method as used for a)?

storetable.mw

Download Analisa_Dinamik_Limb_v1_(30).mwAnalisa_Dinamik_Limb_v1_(30).mw

So I have an equation that basically takes the component of vectors to be used as an equation. The variables that I after are FB1z, FB2x, and FB3y For example here is my equation: 

EOM1:=(AFB1[1]+AFB2[1]+AFB3[1])=TEOM[1]

EOM2:=(AFB1[2]+AFB2[2]+AFB3[2])=TEOM[2]

EOM3:=(AFB1[3]+AFB2[3]+AFB3[3])-TEOM[3]:

FBBp1:=FBPP1=(EulP1[1]+EulP2[1]+EulP3[1]):
FBBp2:=FBPP2=(EulP1[2]+EulP2[2]+EulP3[2]):
FBBp3:=FBPP3=(EulP1[3]+EulP2[3]+EulP3[3]):

However there are unknown variable in AFB2[1] named FB2x and AFB3[1] named FB3y. Then AFB1[2] has unknown equation named FB1z and AFB3[2] has FB3y and so on. While in my FBBp1,FBBp2,and FBBp3 holds all of the variable of FB1z, FB2x, and FB3x
I have tried to use 'solve' command to find the variable but my computer won't stop processing it:
sls:=solve({EOM1,EOM2,EOM3,FBBp1,FBBp2,FBBp3},{FB1z,FB2x,FB3y}):

I tried to use the Gauss-Elimination by forming a matrix but it doesn't work as well since I am really confused how to take out the variables out of the vector component.

zzz:=Matrix([0,AFB2[1],AFB3[1],jjj[1]],[AFB1[2],0,AFB3[2],jjj[2]],[AFB1[3],AFB2[3],0,jjj[3]],[FBP1[1],FBP2[1],FBP3[1],EulP[1]],[FBP1[2],FBP2[2],FBP3[2],EulP[2]],[FBP1[3],FBP2[3],FBP3[3],EulP[3]]):
GaussElimination:=(zzz)

I would be very grateful If someone could help me. Thankyou

Edit: here are the .txt files and .mpl files that required to run the program 

Inverse_Kinematics_ADRIAN2.mw

RotInertiax0_ADRIAN.txt

Download DisplacementXYZ.txt

inersia_platfrom.txt

There is an .mpl file that I couldn't upload so I will upload it in the comments

I have to find the distance between a point and a plot. The point is (2,-3,4) and the plane is x+2y+2z=13. How do I plot this?

i need to find a y(2) and than work with it in a loop

but i cant do it because i solved my equation with "dsolve" and "sol(2)" and its solution is list

really need a help