Can someone help me whith this plot? Somehow Maple always says unable to evaluate the function to numeric values in the region.

File is attached.

Plots.mw

The equaltion given below has 4 roots. I only want to see the 2nd solution for each alpha and it will come in a list. I don't want to see all the solutions. What to do

for alpha from .4 by .1 to 5 do
S := solve(x^4-alpha, x);
S[2] ;
end do;

Dif_Out_plot.mwWhy this difference in the results of the graphics?There are differences between the command smartplot and the command smartplot implicitplot, to display the graph of the function.Following the line of site: web

Hi all,

I need to solve this linear system:

 

R11+R12+fx=1

R21+R22+fy=1

R31+R32+fz=1

R11+R21+R31=1

R12+R22+R32=1

R11*R22*fz+R12*fy*R31+fx*R21*R32-R31*R22*fx-R32*fy*R11-fz*R21*R12=1

 

It is a solution for a 3D rotation matrix in which I know the last coloumn (fx,fy,fz) and I need to find the other six elements of the matrix R11,R12...Obvioulsy there are infinite solution, but in particular...

I am currently doing a project for mass transfer class. Here is a differential equation that I need to solve:

For the solution, and plug in the value of C_AS and C_A0, I finally got the concentration profile of:

Then I write the commands and Maple won't give me a plot (no warnings appeared). Can someone help me with this situation?

Recently, I encountered the problem in which one's asked to investigate if the following fraction is an integer or not. Of course the problem has an theoritical approach, but I wanted to make the problem dynamic. The fraction is p=(k^2-87)/(3k+117) and we want to find k in which p is an integer. What I did is to use the function msolve and a loop for but it gives my just the empty set which is incorrect. Any hint? Thanks

I have a command in Maple that looks like

Hey,

i got some linear optimization problem which i solved via maximize from the simplex library and i'd like to know if there's some way to generate only the final simplex tableau.

After typing

with(LinearAlgebra):
avector:=Vector(3,symbol=a):
bvector:=Vector(3,symbol=b):
'a x b'=CrossProduct(avector,bvector)

Maple writes out

Hi,

iam trying to define a commutator rule (using the Physics package) but iam getting an error i do  not understand.

 

I have (quantum) operators like this: A_ij  i and j are numbers from 0 to 3

The commutator relation i want to define is: [A_ij,A_kl] = A_ilδ_jk - A_kjδ_il

 

My worksheet looks like this:

 

what am i doing wrong???

> restart;
> with(PDEtools);
> a := 3;
> U := u(r, t);
> wave := a^2*(diff(U, r$2))+(diff(U, r))/r) = diff(U,t$2));
> ics := u(1, t) = 0;
> bcs := diff(U, t) = piecewise(0<= r <b,-2, b<= r<1,0), u(r, 0) = 0;
> s := pdsolve(wave, {bcs, ics});

Error, (in pdsolve/sys) too many arguments; some or all of the following are wrong: [{u(r, t)}, {u(1, t) = 0, u(r, 0) = 0, diff(u(r, t), t) = -2}]

HI,is it possible to solve the ODE system (that I have uploaded) by adding events at a certain timestep, instead of reformulating the system for each step.The ODE system expresses two coupled movements y(t) and X(t). The current position of y(t) is controlled by a parameter which depends on values calculated in the previous step, while changing the sign for each step.Thank you! 

We have the following sequence of natural numbers
5, 7, 11, 13, 15, 19, 21, 29, 31, 33, 37, 39, 45, 55, 57, 63, 83, 85, 87, 91, 93, 99, 109,
 111, 117, 135, 163, 165, 171, 189, 245, 247, 249, 253, 255, 261, 271, 273, 279, 297, 325,
 327, 333, 351, 405, 487, 489, 495, 513, 567, 731, 733, 735, 739, 741, 747, 757, 759, 765,
 783, 811, 813, 819, 837, 891, 973, 975, 981, 999, 1053, 1215, 1459, 1461, 1467, 1485, 1539,
 1701, 2189, 2191, 2193, 2197,...

Hello:
    I  want to get the variational value of omega when KR[N]>0.2, but  the omega of my procedure remain the same,  how do I change my program to get the variational value of omega when KR[N]>0.2.

Thank you !

My procedure:
restart:
> Digits:=30;
>
> h0:=0.156;
> d:=0.32*h0;
> l:=2;
> h1:=h0-d;
> h2:=h0+d;
> h3:=0.6*h0;
> g:=9.8;
> d1:=1;

rsolve({g(m, n) = g(m − 1, n) + g(m, n − 1)},g,'genfunc'(a,b));

why no solution