Hi,

according to my previous question

http://www.mapleprimes.com/questions/201435-ODE-With-Constraint

I wrote the following code. at first, the code solve the equation for f and when it slves that, I want to solve Theta in such a way that use the values of f in previous calculation. I use the command 'known' but i couldnt find thesolution

I would be most grateful if you could help me in this problem

Thanks for your attentions in advance

restart; # Notice that Restart (capital R) has no effect (to catch that use semicolon, not colon)
a:=0.13:
b:=0.41:
reynolds:=1.125*10^7;  
Eq1:=diff(f(x),x$3)+diff(f(x),x$2)*f(x)+b^2*sqrt(2*reynolds)*diff(diff(f(x),x$2)^2*x^2,x$1);
Eq2:=diff(g(x),x$3)+diff(g(x),x$2)*g(x)+c*a^2*sqrt(2*reynolds)*diff(diff(g(x),x$2)^2*x,x$1);
eq1:=isolate(Eq1,diff(f(x),x,x,x));
eq2:=subs(g=f,isolate(Eq2,diff(g(x),x,x,x)));
EQ:=diff(f(x),x,x,x)=piecewise(x<c*0.1,rhs(eq1),rhs(eq2));


c:=75:
;
Q:=proc(pp2) local res,F0,F1,F2;
print(pp2);
if not type(pp2,numeric) then return 'procname(_passed)' end if:
res:=dsolve({EQ,f(0)=0,D(f)(0)=0,(D@@2)(f)(0)=pp2},numeric,output=listprocedure);
F0,F1,F2:=op(subs(subs(res),[f(x),diff(f(x),x),diff(f(x),x,x)])):
F1(c)-1;
end proc;


fsolve(Q(pp2)=0,pp2=(0..102));
se:=%;
res2:=dsolve({EQ,f(0)=0,D(f)(0)=0,(D@@2)(f)(0)=se},numeric,output=listprocedure);
G0,G1,G2:=op(subs(subs(res2),[f(x),diff(f(x),x),diff(f(x),x,x)])):

plots:-odeplot(res2,[seq([x,diff(f(x),[x$i])],i=0..2)],0..2); #This plots from and past 0.1*c
pr:=1;
prt:=0.89;

Eq11:=diff(theta(x),x$2)+pr*diff(theta(x),x$1)*f(x)+pr/prt*b^2*sqrt(2*reynolds)*diff(diff(f(x),x$2)*diff(theta(x),x$1)*x^2,x$1);
Eq22:=diff(g(x),x$2)+pr*diff(g(x),x$1)*f(x)+pr/prt*a^2*c*sqrt(2*reynolds)*diff(diff(f(x),x$2)*diff(g(x),x$1)*x^1,x$1);
eq11:=isolate(Eq11,diff(f(x),x,x));
eq22:=subs(g=theta,isolate(Eq22,diff(g(x),x,x)));
EQT:=diff(theta(x),x,x)=piecewise(x<c*0.1,rhs(eq11),rhs(eq22));


QT:=proc(pp3) local res3,theta0,theta1;
print(pp3);
if not type(pp3,numeric) then return 'procname(_passed)' end if:
res3:=dsolve({EQT,theta(0)=1,D(theta)(0)=pp3,known=f},numeric,output=listprocedure);
theta0,theta1:=op(subs(subs(res),[theta(x),diff(theta(x),x)])):
theta0(c);
end proc;

fsolve(QT(pp3)=0,pp3=(0..200));
res3(0);

Amir

in maple 15

https://drive.google.com/file/d/0B8F2D27rfQWgVXE1alN0V3JWU1U/edit?usp=sharing

there are 3 equation to be minimized

and i limit x between x + 5 and x - 5 as constraints

though f1 got a error in first line of command,

later i type a correct command for f1 in later part of script

Pokaż pełną treść

What is the difference between the "old" Maplet format and the "new" Möbius Project App format? Is the Möbius Project just a new name for writing Maplets? Is the programming different; are the mechanics of putting the file in a form to be used by students different; does the file containing the program that can be directly run by students have a different format and extension; are the kinds of things students see different? Is the Maplets function on the way to being discontinued in favor of the Möbius Apps?

I cannot find anything in the documentation for either format that relates it to the other format, although the general description of the output seems entirely the same.

Hi:

how find the coefficients C1,C2?
q(0)=0,D(q)(0)=0

q(T) := sin(T)*_C2+cos(T)*_C1-(1/3000)*Pi*(4012562293500*Pi^3*cos(T)^3-32100498340000*Pi^3*cos(T)^2-3009421720125*Pi^3*cos(T)+16050249170000*Pi^3+435778855000*Pi*cos(T)^2-217889427500*Pi-3539762622):

We know that the set f:={(1,2),(2,a),(3,b)} can introduce a function from {1,2,3} to {2,a,b}. I want Maple to know f as a function. Is this job possible at Maple? I thought to find the Cartesian product of above latter sets and then try to select "f" as one of its subset but this did not help me to force Maple to take "f" as a function. Indeed, I want to work with this kind of function (like plotting and doing f+g, f-g, fog for two functions for example).

Thanks for your time and help.

how to plot tail of normal distribution

would like to see the shape of tail of distribution

 Hi all,

i want to solve a system of differential equations in maple with "dsolve, numeric",

i got th error message: "Error, (in f) unable to store in Matrix ..."

what is the problem?

thanks a lot

Suppose (1 + 2x)^n = a0 + a1*x + a2*x^2+...+an*x^n.

I want to find value of n so that max(a0, a1, ..., an) is a8. 

I tried directly. 

With n = 12. 

restart:
A:=expand((1+2*x)^12,x);
max(coeffs(A));

And with n = 11. 

B:=expand((1+2*x)^11,x);
max(coeffs(B));

Therefore,  n = 12 or  n = 11. 

How can I solve the problem with Maple?

What is the command to get a sequence of the first twenty prime numbers

why system inconsistent 

when only want to obtain a equation H

in terms of a1 , a3, b2, b3, c1, c2, x1, x2, lambda

pdsolve([
Diff(H(x1,x2,lambda), x1) = a1*x1 + a3*lambda,
Diff(H(x1,x2,lambda), x2) = b2*x2 + b3*lambda,
Diff(H(x1,x2,lambda), lambda) = c1*x1 + c2*x2
], H(x1,x2,lambda));

Hi all. I am using Maple 18.

1. Does Maple only output an answer, or can it be made to display the steps it took to achieve the result, in a way that is logical and understandable by the user? For example, I would like to see it display steps when integrating or differentiating, or dividing a polynomial by another polynomial. Basically any process that involve simplifyinf or factorisation.

2. I am trying to divide x^2 + 5x + 9 by x+2, but all Maple displays is the expression

How can I make it do the polynomial division and output the steps it took? I tried using simplify and evaluate from the context menu but the output is still just the expression.

3. Has something happened to this site? It was working fine in Opera until two days ago, and has sinmce started displaying "Error generating page" where the replies to questions should be. No problems with the "posts" section, only in the "questions" section.

Thanks in advance.

Is there a document with all elementairy plottools commands to find?

i must construct a numberline and lines to e4xplain with a arrow and text under it.

with(Finance):
with(Statistics):
X:=RandomVariable(Binomial(nn, pp));
Drift(X);

How to listplot two vectors on the same graph?

That means having two lines.