fcns:=[u(x,y,t),v(x,y,t),h(x,y,t),c(x,y,t)]:

gr:=0.5:pr:=0.71:sc:=0.7:m:=1.0:k:=0.3:
  fcns:=[u(x,y,t),v(x,y,t),h(x,y,t),c(x,y,t)]:
  IC := [u(x,y,0)=0,v(x,y,0)=0,h(x,y,0)=0,c(x,y,0)=0]:
  BC:=[u(0,y,t)=0,h(0,y,t)=0,c(0,y,t)=0,u(x,0,t)=1,v(x,0,t)=0,h(x,0,t)=1,c(x,0,t)=1,u(x,10,t)=0,h(x,10,t)=0,c(x,10,t)=0];
  eq1:={diff(u(x,y,t),t)+u(x,y,t)*diff(u(x,y,t),x)+v(x,y,t)*diff(u(x,y,t),y)=diff(u(x,y,t),y$2)+gr*h(x,y,t)+gr*c(x,y,t)-m*u(x,y,t)
,diff(h(x,y,t),t)+u(x,y,t)*diff(h(x,y,t),x)+v(x,y,t)*diff(h(x,y,t),y)=1/pr*diff(h(x,y,t),y$2),diff(c(x,y,t),t)+u(x,y,t)*diff(c(x,y,t),x)+v(x,y,t)*diff(c(x,y,t),y)=1/sc*diff(h(x,y,t),y$2)-k*c(x,y,t)}:
  pds:= pdsolve(eq1,IC,BC,fcns,numeric):
  pds:= pdsolve(eq1,IC,BC,fcns,numeric,spacestep = 1/100):

for the above problem i made this code.

Hi MaplePrimes Web managers

I suspect you have forgotten to update the lists for MapleSim "2017" you seem to stop at "MapleSim 2016", while I hear from your local office that MapleSim 2018 is "just around the corner" ...
It would make it easier to sort the issues for the last version, no ?

Thanks for your support
Sincerely
Ivar

Hi All!

Is there a way to

1) sechedule the job

2) make it parallel

for computing Ax=B where A is a matrix (n*n) and B is a vector (n).

thanks.

Lamba is meant to be reserved for liouville function plotting some data i just realized it isn't that, so if the people that make the numtheory package please unprotec lamba in the future, or protec it but changed its proc assignment to the liouville function
function that would really sweet thanks

I wanted to have 3 graphs in 1 coordinate system. The Maple command I used was,

restart;
g := x -> (10*x)/(3+x);

f := x -> (10*x)/(5+x);

h := x -> (8*x)/(5+x);

plots:-display(
  plot( [g(x),f(x), h(x)], x=0..30, color=[red,blue, green], legend = ["k=3,c=10", "k=5, c=10","k=5,c=8"]));

Warning, unable to evaluate 1 of the 3 functions to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

The error comes up. only f(x) and h(x) were shown. How do I fix this? Thank you.

I 'm solving asystem of equations and one of the result are Hfloat what does it means and how i convert it to number???

how to find skin friction value below code

restart

PDEtools[declare]((U, W, T, C)(y), prime = y):

R1 := .1; R0 := .1; m := .1; a := .1; Ha := .1; Nt := .1; Nb := .1; Pr := 6.2; Le := .6; Bi := 1; Ec := .1; k := 1; r := .1; A := 1;

sys := diff(U(y), `$`(y, 2))+(R1*(diff(U(y), y))-2*R0*W(y))*exp(a*T(y))-a*(diff(U(y), y))*(diff(T(y), y))-Ha = 0, diff(W(y), `$`(y, 2))+(R1*(diff(W(y), y))+2*R0*U(y))*exp(a*T(y))-a*(diff(W(y), y)) = 0, diff(T(y), `$`(y, 2))+R1*Pr*(diff(T(y), y))+Pr*Ec*exp(-a*T(y))*((diff(U(y), y))*(diff(U(y), y))+(diff(W(y), y))*(diff(W(y), y)))+Pr*Ha*Ec*((U(y)+m*W(y))*(U(y)+m*W(y))+(W(y)-m*U(y))*(W(y)-m*U(y)))/(m^2+1)^2+Nb*(diff(T(y), y))*(diff(C(y), y))+Nt*(diff(T(y), y))*(diff(T(y), y)) = 0, diff(C(y), `$`(y, 2))+Pr*Le*R1*(diff(C(y), y))+Nt*(diff(C(y), `$`(y, 2)))/Nb = 0:

ba := {sys, C(0) = 0, C(1) = 1, T(1) = 0, U(0) = 0, U(1) = 0, W(0) = 0, W(1) = 0, (D(T))(0) = Bi*(T(0)-1)}:

r1 := dsolve(ba, numeric, output = Array([0., 0.5e-1, .10, .15, .20, .25, .30, .35, .40, .45, .50, .55, .60, .65, .70, .75, .80, .85, .90, .95, 1.00])):

with(plots); 

p1u := odeplot(r1, [y, U(y)], 0 .. 1, numpoints = 100, labels = ["y", "U"], style = line, color = green); 

plots[display]({p1u})

I have modeled a 3D flexible manipulator in Maplesim. I intend to measure the torsional strain of the flexible beams. How can I do this??

restart; with(plots);
G := 1; M := .1; R := 1; P := .72; alpha := .1; phi := 1; K := 1; n := 2; beta := 1;
                               1
                              0.1
                               1
                              0.72
                              0.1
                               1
                               1
                               2
                               1
ode1 := {(1+(4/3)*R)*(diff(theta(x), x, x))+(1/2)*P*f(x)*(diff(theta(x), x))+alpha*theta(x) = 0, n*(diff(f(x), x, x))^(n-1)*(diff(f(x), x, x, x))+f(x)*(diff(f(x), x, x))/(n+1)+G*theta(x)-M*(diff(f(x), x)) = 0, f(0) = 0, theta(10) = 0, (D(f))(0) = beta*K*((D@@2)(f))(0), (D(f))(10) = 1, (D(theta))(0) = -phi*(1-theta(0))};
 /7  d  / d          \                     / d          \
{ - --- |--- theta(x)| + 0.3600000000 f(x) |--- theta(x)|
 \3  dx \ dx         /                     \ dx         /

                         / d  / d      \\ / d  / d  / d      \\\
   + 0.1 theta(x) = 0, 2 |--- |--- f(x)|| |--- |--- |--- f(x)|||
                         \ dx \ dx     // \ dx \ dx \ dx     ///

     1      / d  / d      \\                  / d      \      
   + - f(x) |--- |--- f(x)|| + theta(x) - 0.1 |--- f(x)| = 0,
     3      \ dx \ dx     //                  \ dx     /      

  f(0) = 0, theta(10) = 0, D(f)(0) = @@(D, 2)(f)(0),

                                           \
  D(f)(10) = 1, D(theta)(0) = -1 + theta(0) }
                                           /
dsol := dsolve(ode1, numeric, method = bvp[midrich], range = 0 .. 10);
Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

Hello,

I need a little help...

what I want to do, if i have 3 equations  with multiple variables and i want to get derivative of some of these variables wrt other variables while setting others as constants

here in the example

1st i thought the last equation should have partial m / partial f ... how can i get that

2nd if ,another problem, i wanted to set z as a constant so it won't consider getting its derivative, how do i get that

diff_test.mw

Download diff_test.m

hope i was clear in my explanation. thank you

Hello,

So my questions is fairly simple but I am not sure if Maple can do this. I want to create a Matrix full of partial derivatives (Jacobian), but computing this Jacobian is maxing out my RAM. So I would like to create each row (or column) of the Jacobian seperately and save each row (or column) in a Matrix that has already been saved to a file. Basically I would like to store values in a file (in a loop) without overwriting the file.

Hi every body:

I have to obtain A_max where is the maximum value of the f(t) function at steady state but at first i don't know steady sate what's the mean and the second i don't know how to obtain it with maple. help me plz. tnx. 

f:=(t)->A*sin(x*t)

A := -6000.*x^3*((5.925123867*10^87*I)*x+5.830850605*10^88*x^2+1.752916719*10^90)/((4.241561702*10^89*I)*x^4-2.669934485*10^90*x^5-(3.812926555*10^90*I)*x^2+7.038145625*10^91*x^3+3.174301659*10^89*I-8.376628660*10^90*x)

Sin.mw

I am looking for advice as to avoiding or identifying what may elude to the cause of the error "Error, (in simplify/do) invalid simplification command" in circumstances where i have simply enclosed an expression with simplify(expression), and in most cases an error is not returned.

Thankyou in advance