solve Derivative in maple

see in photo

What do people find a good screen size is? 

Personally I find two monitors would be ideal for using Maple and especially Maplesim.  I'm going to be adding a second monitor, probably a 20".  Maybe 24" is better?  I think 20" would be enough.  Just wondering how many monitors people out there work with and their monitor sizes.

Dear all

I have n, p two integers greater than one

I would like to minimize the following real number by a fixed positive constant if its possible

(n/(n+p))^(n/p)-(n/(n+p))^((n+p)/p)

Many thanks

Dear all,

I am creating an animation, and I was wondering if I can add a multiplier to my equation in a specific range.  So starting from z:=0.2 add a multiplier (z+1). The code I have so far is added. Does anyone know a code for this?

Kind regards

restart; 
with(plots); 
a := -1/2; b := 1/2; c := -2; d := 2; n := 20; 
g := proc (x) options operator, arrow; value(Int(sigma(t), t = 0 .. x)) end proc; 
sigma := proc (z) options operator, arrow; 2*sqrt(2*h^2-4*z^2)*z/h^2 end proc; 
h := i/n; 
for i to n do 
an2[i] := plot(sigma(z), z = -(1/2)*h .. (1/2)*h, view = [a .. b, c .. d], color = AQUAMARINE); 
an3[i] := plot(2*g(x), x = 0 .. (1/2)*h, view = [a .. b, c .. d], color = RED) 
end do; 
p := plots[display]([seq(an2[i], i = 1 .. n)], insequence = true); 
q := plots[display]([seq(an3[i], i = 1 .. n)], insequence = true); display(p, q)

equation 1 : x_i+1=x_i− (f·g_y−f_y·g)/(f_x ·g_y −f_y ·g_x)
equation 2: y_i+1=y_i− (fx·g−f·g_x)/(f_x·g_y-f_y·g_x)

My quesiton are, deriving equations (1) and (2) above and constructing a single Maple function called newt2d that implements both of these recurrence relation.

I apolgize in advance if I don't write my question correctly.  This is my first time posting a question. 

> restart;

with(plots);

pr := .72; p := 0; n := [.5, 1, 1.5]; s := 0; a := .2; b := 0; L := [red, blue, green]; l := 0; k := 1;

for j to nops(n) do R1 := 2*n[j]/(1+n[j]); R2 := 2*p/(1+n);

sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2) = 0, diff(diff(theta(eta), eta), eta)+pr*k*f(eta)*(diff(theta(eta), eta))+R2*pr*k*(diff(f(eta), eta))*theta(eta)+(2*(a*(diff(f(eta), eta))+b*theta(eta)))/(1+n[j]) = 0, f(0) = 1, (D(f))(0) = b*((D@@2)(f))(0), (D(f))(1.8) = 0, theta(0) = 1+s*(D(theta))(0), theta(1.8) = 1], numeric, method = bvp);

fplt[j] := plots[odeplot](sol1, [eta, diff(diff(f(eta), eta), eta)], color = L[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = L[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

plots:-display([seq(tplt[j], j = 1 .. nops(n))]);

staganation_point1.mw
 

can we chage the axis sir ?? like  f'' vs eta to f'' vs lambda.

Download staganation_point1.mw

Hello every one,

I'm using Maple to create a compiled program of my model Maplesim. I had a error with "getcompiledProc" command, which is :

Error, (in GetCompiledProc) non-numeric initial condition for `Main.RAB.value`(t): Float(undefined)

the "Main.RAB.value(t)" represent a probe used to visualize an output of my system.

Does anybody had a such error ?
Thank you for your responses.

Larbi ANIBA

I am writing a package to define a few Records and collect some quantities and formulae to do some arithmetic with.

At some point I started running into a severe problem: Each time I want to use one of the exports form the module, the Maple kernel dies with a "Stack limit exceeded" message. The issue is: This seems to happen only on Maple 2016. I had some issues in Maple 2015 but that got sorted by rewriting some stuff. I never had any issue like that with Maple 15. I went as far as commenting out most of my .mapleinit code to make sure that did not cause issues. The code works as intended in Maple 15 and in Maple 2015.

I am posting this here and am asking others to see if they can reproduce this. Needless to say I will SCR this. But I actually need to get this to work.

The error can be triggered by simply entering

Booster;

after loading the package. Needless to say, it should just echo the name "Booster", to get more one would need to enter eval(Booster) (due to the last-name evaluation rules).

Thanks,

M.D.

BooMBA.mw

Hello everybody,

I want to solve this pde.the desire solution is V(r,z). three boundar conditions are written that two of them are related to rhe radial and one is related to the longitudinal coordinate.

I attached the solution for you. but this solution is derived by Matlab. Now, I just want to resolve it by Maple, but I couldn't reach it. Please let me know the correct way asap.

Thanks a lot.PDE.mw
 

Download PDE.mwPDE.mw

Howdy!!

Just thought I'd put a quick one in. I'm trying to solve the Schrodinger equation given as:

under the constraint

I'm trying to plot this for a range of values of the parameters a,b and fixed p=0,1.

I'm aiming to get a solution of the form

and plotting it.

Thanks in advance

Recently I have used a FilterForge filter to generate hyperbolic tiles from images (Lua code with Möbius transform etc. available). If an input image (1 see below) is selected a tile (2) is generated if the hyperbolic parameters (5,3,2) were chosen.  How could this be made with Maple?

> restart; for j to nops(n) do sys := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))+1-(diff(f(eta), eta))^2 = 0, (diff(diff(theta(eta), eta), eta))/pr+f(eta)*(diff(theta(eta), eta))-(diff(f(eta), eta))*theta(eta) = 0; bcs := f(0) = 0, (D(f))(0) = l+b*((D@@2)(f))(0), (D(f))(-.5) = 1, theta(0) = 1+s*(D(theta))(-.5), theta(2) = 0; n := [1, 2, 3, 4, 5, 6]; pr := .71; p := 0; q := 0; b := 0; l := 0; s := 0; L := [red, blue, orange]; R1 := 2*n[j]/(1+n[j]); R2 := 2*p/(1+n); p := proc (f1, th1, { output::name := 'number' }) local res1, fvals, thvals, res2; option remember; res1 := dsolve({sys, f(1) = 0, theta(0) = 1+th1, (D(f))(0) = f1, (D(theta))(0) = th1, ((D@@2)(f))(0) = f1-1}, numeric, :-output = listprocedure); fvals := (subs(res1, [seq(diff(f(eta), [`$`(eta, i)]), i = 0 .. 2)]))(0); thvals := (subs(res1, [seq(diff(theta(eta), [`$`(eta, i)]), i = 0 .. 1)]))(0); res2 := dsolve({sys, f(0) = fvals[1], theta(0) = thvals[1], theta(5) = 0, (D(f))(0) = fvals[2], (D(f))(5) = 1}, numeric, :-output = listprocedure); if output = 'number' then [fvals[3]-(subs(res2, diff(f(eta), `$`(eta, 2))))(0), thvals[2]-(subs(res2, diff(theta(eta), eta)))(0)] else res1, res2 end if end proc; p1 := proc (f1, th1) p(args)[1] end proc; p2 := proc (f1, th1) p(args)[2] end proc; p(.3, -.2); par := fsolve([p1, p2], [.3, -.2]); res1, res2 := p(op(par), output = xxx); plots:-display(plots:-odeplot(res1, [[eta, f(eta)], [eta, theta(eta)]]), plots:-odeplot(res2, [[eta, f(eta)], [eta, theta(eta)]])); plots:-display(plots:-odeplot(res1, [[eta, diff(f(eta), eta)], [eta, diff(theta(eta), eta)]]), plots:-odeplot(res2, [[eta, diff(f(eta), eta)], [eta, diff(theta(eta), eta)]])); plots:-display(plots:-odeplot(res1, [[eta, diff(f(eta), eta, eta)]]), plots:-odeplot(res2, [[eta, diff(f(eta), eta, eta)]])); plots:-display(plots:-odeplot(res1, [[eta, diff(f(eta), eta)]])); fplt[j] := plots[odeplot](res1, [eta, diff(diff(f(eta), eta), eta)], color = L[j], axes = boxed); tplt[j] := plots[odeplot](res1, [[eta, theta(eta)]], color = L[j], axes = boxed) end do;
> plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

Dear Sir

In my above problem i trying to plot for set of values of n but in plot command it not executing , can you do this why it is not executing ??

ha := (diff(c(t), t))/(c(t)*(diff(c(t), t))-c(t));
solve(subs(m=ha,f(m))*subs(m=subs(c(t)=a(t),ha)), f(m)) = subs(m=ha+subs(c(t)=a(t),ha), f(m), f);

just expect to find a function ?
 

1. How to generate times series data that have some logic law

but do not include distributive law?

Hi guys! 

I have a PDE system. The mayority of the equations are equal to zero, but two of them are:

where a, b, c, d are CONSTANT parameters. I know that if a=b=d=c=1 the system is inconsistent. But I also know that if a=-1, b=d=0 and c=1 the system is consistent and exist the solution. I wanna know if there's a way to ask maple to find another selections of my parameter that make my PDE consistent and what it's the solution for that selection of a,b,c,d. 

Here's my PDE system (sys2). 

test.mw

thank you so much for your time!