with plots

Eq[1] := 5*(diff(f(eta), `$`(eta, 3)))+(m+3)*f(eta)*(diff(f(eta), `$`(eta, 2)))-(2*m+1)*(diff(f(eta), eta))^2-(2*(2*m+1))*h(eta)-(m-2)*eta*(diff(h(eta), eta))

Eq[2] := g(eta)-(diff(h(eta), eta))

Eq[3] := 5*(diff(g(eta), `$`(eta, 2)))/Pr+(m+3)*f(eta)*(diff(g(eta), eta))-5*m*g(eta)*(diff(f(eta), eta))

BCs := [f(0), (D(f))(0), (D(f))(18), g(0)-1, g(18), h(18)]

pars := {Pr = .1, m = 0}

for i to 3 do eq[i] := subs(pars, Eq[i]) end do

hi. i write a program consists a for loop in which cycle 'i', x[i] and hz[i] is calculated.

now i want to plot hz[i] correspond to each x[i].

this code is working but how to brief them?

hhz := piecewise(x < X[1], hz[1], X[1] < x and x < X[2], hz[2], X[2] < x and x < X[3], hz[3], X[3] < x and x < X[4], hz[4], X[4] < x and x < X[5], hz[5], X[5] < x and x < X[6], hz[6], X[6] < x and x < X[7], hz[7], X[7] < x and x < X[8...

So, I've defined procedures in a .mw file and now I want to somehow read this file into another worksheet, in order use the aforementioned procedures. Is this doable? I've unsuccessfully tried using the READ command. The error messages I get are :

Error, on line 1, syntax error, character `?` unexpected:

<?xml version="1.0" encoding="UTF-8"?>

Error, while reading ``file.mw``

Not sure where to go from here. Any help would be great!

Hi,

I am trying to graph some curve C and the vectors from the vector field corresponding to t=1 and t=-1.

My Vector Field is defined by some F(x,y,z) and my Curve is some r(t) [didn't bother putting the actual values here since I only want a general expression...]

In other words I want to have some kind of a curve in 3D with vectors coming out at points t=+/-1 along the curve itself. Is there any way to do such thing?

Thanks!

Hi,

I would like the title of my y axis to say "SaO_² (%)".  When I put in "SaO[2] (%)" the output is SaO[2] (%); However, when I put in SaO[2] (%) the output is "SaO₂ PLOT (...) "  How can I get both parts correct (i.e., the subscript of the 2 and the displaying of the percent sign)?

Thanks!

Hello, I have tried to solve a system of ODE about the gyroscope for a week. I have entered the following commands : >eq1:=(l^2+r^2/2)*(sin(theta(t))^2*diff(diff(psi(t),t),t)+2*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t)))-r^2*(cos(theta(t))*diff(psi(t),t)+omega(t))*diff(theta(t),t)=0;>eq2:=(l^2+r^2/2)*(diff(diff(theta(t),t),t)-(diff(psi(t),t))^2*sin(theta(t))*cos(theta(t)))+r^2*(diff(psi(t),t)*cos(theta(t))+omega(t))*diff(psi(t),t)*sin(theta(t))=g*l*sin(theta(t)); >eq3:=diff(diff(psi(t...

I want to know for what combinations of parameters (x,y) the following holds: Prob(X(x) + Y(y) < z) = alpha. X and Y are random variables (I want to look at different distributions), z and alpha are constants. I tried the following (simplified with X=Y):

Dear Sir / Madam,

Would you recommend me how to isolate the real and imaginary part and simplify the equation A? Maple cannot evaluate this equation to numeric value...

A := evalf(sqrt((Re(B))^2 + (Im(B))^2));
A := evalf(abs(B));

are not working.

I have attached my equation function.txt

f(x)=(sin@@3)(0.5143600172e-2+1.000000039*arcsin(arcsin(sin(sin(-0.1543083839e-1+.9999999610*RootOf((sin@@3)(_Z)-(sin@@3)(0.5143619410e-2+1.000000000*RootOf(-(sin@@3)(_Z)+x))))))))

i want to find inverse of this funtion i am getting nothing

>s := solve(w = f(x), x)

    s:=

i am getting like above.

i want to find inverse of above function plz give me reply to me.

f :=  (0.1953125000e-2+.3750000000*(RootOf(-512*x-7+192*_Z^(2/3)+24*_Z^(1/3)+512*_Z-6*(512*_Z+192*_Z^(2/3)+24*_Z^(1/3)+1)^(2/3)+12*(512*_Z+192*_Z^(2/3)+24*_Z^(1/3)+1)^(1/3))^(2/3))+0.4687500000e-1*(RootOf(-512*x-7+192*_Z^(2/3)+24*_Z^(1/3)+512*_Z-6*(512*_Z+192*_Z^(2/3)+24*_Z^(1/3)+1)^(2/3)+12*(512*_Z+192*_Z^(2/3)+24*_Z^(1/3)+1)^(1/3))^(1/3))+RootOf(-512*x-7+192*_Z^(2/3)+24*_Z^(1/3)+512*_Z-6*(512*_Z+192*_Z^(2/3)+24*_Z^(1/3)+1)^(2/3)+12*(512*_Z+192*_Z^(2/3)+24*_Z^(1/3)+1)^(1/3)))^3

Hi all

I have used pointplot3d to plot a set of points using three vecotrs v1, v2 and v3 which all contain many elements and the same number of elements.

I would like to know how to plot a new set of elements for new vectors v4,v5, and v6 on to the same plot as the previous set of values.

Could I then give each of the two sets a different color?

I am to solve this differential equation (nonlinear):

C'(t) = G[0]+0.1*sqrt(C(t))-C(t)*n+C[ude]*n

G[0] = 1.7*1.2*83

n = 6

C[ude] = 417

i need to use programming since the DE is nonlinear.

Does anybody have the solution to that?

Thanks in advance.

Best regards,

Marie

f:=x->x^3:
y[1]:=0.5:
y[3]:=subs(x=y[1],f(x)):
i := 2:
y[i] := y[i-1]-(y[1]-y[2+1])*(1/2):
i:=1:
g[i] := (y[i+1]-y[i+2])*x/(y[i]-y[i+1])+(y[i+2]*y[i]-y[i+1]*y[i+1])/(y[i]-y[i+1]):
g[i] := unapply(g[i], x):
h[i] := unapply(solve(y = g[i](x), x), y):
for i from 2 by 1 to 7 do
g[i]:=(f@h[i-1])(x);
g[i]:=unapply(g[i],x);
y[i+2]:=subs(x=y[i+1],g[i](x));
s:=solve(z=g[i](x),x);
h[i]:=unapply(s[1],z);
end do:

for i from -1 by -1 to -8 do

f:=x->x^3:
y[1]:=0.5:
y[3]:=subs(x=0.5,f(x)):
i := 2:
y[i] := y[i-1]-(y[1]-y[2+1])*(1/2):
i:=1:
g[i] := (y[i+1]-y[i+2])*x/(y[i]-y[i+1])+(y[i+2]*y[i]-y[i+1]*y[i+1])/(y[i]-y[i+1]):
g[i] := unapply(g[i], x):
h[i] := unapply(solve(y = g[i](x), x), y):
for i from 2 by 1 to 5 do
g[i]:=(f@h[i-1])(x);
g[i]:=unapply(g[i],x):
s:=solve(y=g[i](x),x):
h[i]:=unapply(s[1],y);
end do;

 for this code i got error "Error, (in kernels)...

I have a system of equations  AX=B where the matrix A from dimension 307x200, X from 200x1, and B from 307x1
How can I use Gauss elimination to solve these system of equations.