Items tagged with plot3d plot3d Tagged Items Feed

I've been trying to make a smooth plot of some ODEs. It should show C rapidly increasing at the innitiation, until they get into a quassi steady state,  and then all three variables increase much slower. This should look like a roughly straight line that elbows sharply into a smooth curve.

Any attempt to DEplot3d it i've made either just shows the time before the quassi steady state is reached, so shows the straight line; or smooths that time together with the next period, making the straight line look like a part of the smooth curve.


Model := [diff(B[1](t), t) = k[a1]*C(t)*(R-B[1](t)-B[2](t))-k[d1]*B[1](t), diff(B[2](t), t) = k[a2]*C(t)*(R-B[1](t)-B[2](t))-k[d2]*B[2](t), diff(C(t), t) = (-(k[a1]+k[a2])*C(t)*(R-B[1](t)-B[2](t))+k[d1]*B[1]+k[d2]*B[2](t)+k[m]*((I)(t)-C(t)))/h];
DissMod := subs((I)(t) = 0, Model);
AssMod := subs((I)(t) = C[T], Model);


Pars := [k[a1] = 6*10^(-4), k[d1] = 7*10^(-3), k[a2] = 5*10^(-4), k[d2] = 10^(-2), R = .5, k[m] = 10^(-4), C[T] = 100, h = 10^(-6)]

StateSol := DEplot3d(subs(Pars, AssMod), [B[1](t), B[2](t), C(t)], t = 0 .. 1000, number = 3, B[1] = 0 .. .5, B[2] = 0 .. .5, [[B[1](0) = 0, B[2](0) = 0, C(0) = 0]], scene = [B[1](t), B[2](t), C(t)], maxstep = .1, maxfun = 0, method = l)

I have two deformed planes, that i would like to draw with 3dplot, as well as drawing a curve marking their intersection.

the curves are given by the expressions:


C = -(k[d2]*B[2]+I*k[m]+k[d1]*B[1])/((B[1]+B[2]-R)*k[a1]+(B[1]+B[2]-R)*k[a2]-k[m]),

C = k[d1]*B[1]/(k[a1]*(R-B[1]-B[2]))

 

evaluated at

Pars := [k[a1] = 6*10^(-4), k[d1] = 7*10^(-3), k[a2] = 5*10^(-4), k[d2] = 10^(-2), R = .5, k[m] = 10^(-4), C[T] = 100, h = 10^(-6)]

with the variables B[1],B[2] and C within the bounds [0..0.5],[0..0.5],[0..100].

 

My method was to try and use solve to find a formula for the intersection curve- but i couldn't get 3dplot to plot it!

in Mathematica, there is the option called BoxRatios

"is an option for Graphics3D that gives the ratios of side lengths for the bounding box of the threedimensional picture."

It is sort-of like aspect ratio, but for 3D. It is set by default so make 3D plot looks "nice". I can't seem to find equivalent Maple option. The closest is the option "s=" for plot3d, but this just turns of/on "constrained scaling" and does not allow one to modify the "BoxRatios"

Let me give an example. Here is 3D plot in Mathematica and the same in Maple. I'd like to get the Maple 3D to look similar to Mathematica 3D in terms of the "aspect ratio". Maple on the z-axis is using the same size as in the x and y axis, and even though this is realistic, it does not make the plot as nice. I want to change this ratio.

T0[x_, y_, m_] :=20/Pi Sum[ (-1)^(n + 1)/n Exp[- (n Pi/10) y] Sin[ (n Pi/10) x], {n,1, m}]
Plot3D[T0[x, y, 70], {x, 0, 10}, {y, 0, 10}, PlotRange -> All,  AxesLabel -> {x, y, z}]

In Maple:

T0:= (x,y,m)-> 20/Pi*sum( (-1)^(n+1)/n*exp(- n*Pi/10*y)*sin(n*Pi/10*x),n=1..m);
plot3d(T0(x,y,50),x=0..10,y=0..10,scaling=unconstrained);

So Maple is using 1:1:1 box ratio. Mathematica default is 1:1:0.4, and I wanted to see if I can change Maple to be the same.

I get same plot in Maple using scaling=unconstrained or scaling=constrained. So this option is not very useful for what I want.

Is there a way to change the "BoxRatios" as defined above in Maple? There must be, right? Do I need to use different package?

 

Hi,

i would like to plot a graph in R^3 of a function f(n,t), where n is integer and t is real. For every t i would like to have a sequence of points. Is it possible?

thanks, anna rita 

I'm new here, so I'm not totally sure this is the right place to ask this. I apologize if it isn't, please let me know in that case.

 

My problem is that Maple won't recognize the built-in command 'complexplot3d'. For example, typing:

complexplot3d(z^2, z = -1-I .. 1+I)

doesn't do anything, and it's displayed again in blue as if it was not a command.

 

Any kind of help would be extremely helpful, as I have no clue of what's going on. Thanks in advance!

Im trying to draw a shpere but it always saying: 

Error, (in plot3d) unexpected option: z = -2 .. 2


this is the equation: x^2+y^2+z^2-4=0

i'm writing this way

plot3d(x^2+y^2+z^2-2^2, x = -2 .. 2, y = -2 .. 2, z = -2 .. 2)


what should I do? this is my first time with this software

 

best from Brazil,
Nina

Hello,

does anyone know a way to combine two plots in one where one is created with ScatterPlot3D (Package:Statistics) and the other one with plot3D (Package:plots)?

 

Normally you would write something like this:

but that only works if the plots are from the same package...

Hi,

 

I am having trouble in plotting the following surface (it is a quite complicated expression, but should be fine).

OwnSurface := [-Re(arctan(exp(I*Pi*(1/4))*(u^2+v^2)^(1/2)))-(1/2)*ln((1+(u^2+v^2)^2)^(1/2)), -arctan(1/((2*(u^2+v^2))^(1/2)-1))+7*Pi*(1/4), (1/8)*ln(u^2+v^2-(2*(u^2+v^2))^(1/2)+1)-(1/8)*ln(u^2+v^2+(2*(u^2+v^2))^(1/2)+1)-(1/4)*arctan(1/((2*(u^2+v^2))^(1/2)+1))+(1/4)*arctan(1/((2*(u^2+v^2))^(1/2)-1))];

plot3d(OwnSurface, u = -.4 .. .4, v = -.4 .. .4, labels = [x1, x2, x3]);

The only thing maple does is plotting a box with a diagonal line. How can I fix this?

 

Is there an elegant way to plot the region between the surfaces z=-y^2 and z=x^2, only on the domain of the XY-plane bounded by the triangle with vertices (0,0), (1,0) and (1,1)?

I am trying to have the output of DETOOLS as 3dpolarplot. As in the following example:

 

EF := {2*(diff(w[2](t), t)) = 10, diff(w[1](t), t) = sqrt(2/w[1](t)), diff(w[3](t), t) = 0}; with(DEtools); DEplot3d(EF, {w[1](t), w[2](t), w[3](t)}, t = 0 .. 100, [[w[1](0) = 1, w[2](0) = 0, w[3](0) = 0]], scene = [w[1](t), w[2](t), w[3](t)], stepsize = .1, orientation = [139, -106])

 

how can I get the output as a polarplot in 3d where, w[2] and w[3] have range 0..2*pi.

Please help in this respect asap.

 

Hello everyone!

Is this a bug that the following two commands work differently?:

densityplot(sin(x*y), x = -5 .. 5, y = -5 .. 5, colorscheme = ["zgradient", ["blue", "green", "yellow", "red"], zrange = -5 .. 5], style = surface)

plot3d(sin(x*y), x = -5 .. 5, y = -5 .. 5, view = -10 .. 10, colorscheme = ["zgradient", ["blue", "green", "yellow", "red"], zrange = -5 .. 5], style = surface)

The second one works fine in that if you increase the magnitude of sin(x*y) (e.g. 3sin(x*y)) the coloring changes accordingly. But the first one plots sin(x*y) or 5sin(x*y), etc. just the same!

Many thanks for you comments in advance!

So, I have the following inputted into Maple 2015, but the plots are not displaying.  I've double checked the syntax, but I don't see anything wrong.  Any thoughts?  Thank you!

 

plot3d(cos(theta)^2, phi = 0 .. 2*Pi, theta = 0 .. Pi, axes = BOX, style = line, grid = [50, 50], view = [-1 .. 1, -1 .. 1, -1 .. 1], coords = spherical);

 

Trying to plot that here actually displays the plot, but it isn't showing up in my Maple worksheet.

how draw a line x+y = 1 on a sphere x^2+y^2+z^2 = 1 in maple

Let us consider the maximum value of the polynomial

x^4+c*x^2+x^3+d*x-c-1

on the interval x=-1..1 as a function g of the parameters c and d. General considerations suggest its continuity. However, a plot3d of g does not  confirm it.  Also the question arises "Is the function g(c,d) bounded from below?". Here is my try with the DirectSearch and NLPSolve:

 

restart;
``

g(10, -10)

9.

(1)

plot(x^4+x^3+10*x^2-10*x-10-1, x = -1 .. 1)

 

plot3d(g, -5 .. 5, -5 .. 5, grid = [100, 100], style = surface, color = "DarkOliveGreen")

 

DirectSearch:-GlobalOptima(proc (a, b) options operator, arrow; g(a, b) end proc, {a = -1000 .. 1000, b = -1000 .. 1000}, variables = [a, b])

[-167.208333252089, Vector(2, {(1) = 999.9999999975528, (2) = 166.20833325208952}, datatype = float[8]), 815]

(2)

DirectSearch:-GlobalOptima( (a, b) -> g(a, b), variables = [a, b])

DirectSearch:-GlobalOptima(proc (a, b) options operator, arrow; g(a, b) end proc, variables = [a, b])

Error, (in Optimization:-NLPSolve) invalid input: PolynomialTools:-CoefficientVector expects its 1st argument, poly, to be of type polynom(anything, x), but received HFloat(HFloat(undefined))*x^4+HFloat(HFloat(undefined))*x^3+HFloat(HFloat(undefined))*x^2+HFloat(HFloat(undefined))*x+HFloat(HFloat(undefined))

 

``

 

Download bound.mw

 

I have an equation for r(t) that involves 3 (slidable) constants; an equation for phi(t) that involves the same three constants and is written in terms of arctan; theta is a slidable constant. How do I plot this on an x,y,z plot? I want an animation in terms of t.

1 2 3 4 5 6 7 Last Page 1 of 9