Maple 13 Questions and Posts Maple 13 Questions and Posts Feed

These are Posts and Questions associated with the product, Maple 13

Dear friend! Hope you will be fine, I want to general matrix for m by m order. The file attached. I am waiting kind response.

Help.mw

PhD (Scholar)
Department of Mathematics

 

can anybody help me? i want to check the consistency of my scheme. My equation is too long if i check manually, so i used maple 13 to simplify my equation. But it cannot simplify it because of length of output exceed limit 1000000

restart

eqn1 := u+(1-exp(-m))*u[t]+(1-exp(-m))^2*u[tt]/factorial(2)+(u-(1-exp(-m))*u[t]+(1-exp(-m))^2*u[tt]/factorial(2))-u-(1-exp(-m))*u[x]-(1-exp(-m))^2*u[xx]/factorial(2)-u+(1-exp(-m))*u[x]-(1-exp(-m))^2*u[xx]/factorial(2)+(1-exp(-m))^2*u+(1-exp(-m))^2*u^3-(1-exp(-m))^2*(4*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3))*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)/((((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3+(x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)+(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3+(x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3));

(1-exp(-m))^2*u[tt]-(1-exp(-m))^2*u[xx]+(1-exp(-m))^2*u+(1-exp(-m))^2*u^3-(1-exp(-m))^2*(4*((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-16*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+4*(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)/((((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3+(x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)+(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3+(x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3))

(1)

a := simplify(eqn1);

`[Length of output exceeds limit of 1000000]`

(2)

``


Download consistency_expmle_4.mw

.

can anybody help me? i want to check the consistency of my scheme. My equation is too long if i check manually, so i used maple 13 to simplify my equation. But it cannot simplify it because of length of output exceed limit 1000000

restart

eqn1 := u+(1-exp(-m))*u[t]+(1-exp(-m))^2*u[tt]/factorial(2)+(u-(1-exp(-m))*u[t]+(1-exp(-m))^2*u[tt]/factorial(2))-u-(1-exp(-m))*u[x]-(1-exp(-m))^2*u[xx]/factorial(2)-u+(1-exp(-m))*u[x]-(1-exp(-m))^2*u[xx]/factorial(2)+(1-exp(-m))^2*u+(1-exp(-m))^2*u^3-(1-exp(-m))^2*(4*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3))*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)/((((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3+(x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)+(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+(t+1-exp(-m)))-(4*(x+1-exp(-m)))*sinh(x+1-exp(-m)+(t+1-exp(-m)))+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+(t+1-exp(-m)))^3+(x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3));

(1-exp(-m))^2*u[tt]-(1-exp(-m))^2*u[xx]+(1-exp(-m))^2*u+(1-exp(-m))^2*u^3-(1-exp(-m))^2*(4*((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-16*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+4*(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)/((((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3)*((x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3+(x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)+(((x+1-exp(-m))^2-2)*cosh(x+1-exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+1-exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+1-exp(-m)+t)^3)*((x^2-2)*cosh(x+1-exp(-m)+t)-4*x*sinh(x+1-exp(-m)+t)+x^6*cosh(x+1-exp(-m)+t)^3)*(((x+1-exp(-m))^2-2)*cosh(x+2-2*exp(-m)+t)-4*(x+1-exp(-m))*sinh(x+2-2*exp(-m)+t)+(x+1-exp(-m))^6*cosh(x+2-2*exp(-m)+t)^3+(x^2-2)*cosh(x+t)-4*x*sinh(x+t)+x^6*cosh(x+t)^3))

(1)

a := simplify(eqn1);

`[Length of output exceeds limit of 1000000]`

(2)

``


Download consistency_expmle_4.mw

.


can anybody help me..? why my graph not come out? Is that any mistake in my coding?

restart

y := x^2-x*(exp(I*k*`Δx`)+exp(-I*k*`Δx`)-m^2+m^2*((4*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3))*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3)*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3))*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3)/(epsilon*((((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3)*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3))*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3+((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3))+(((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3)*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3))*((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3+((1^2-2)*cosh(1+1)-4*sinh(1+1)+1^6*cosh(1+1)^3)))))+1;

x^2-x*(exp(I*k*`Δx`)+exp(-I*k*`Δx`)-m^2+2*m^2*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(epsilon*(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3)))+1

(1)

subs(m = 1-exp(-m), %);

x^2-x*(exp(I*k*`Δx`)+exp(-I*k*`Δx`)-(1-exp(-m))^2+2*(1-exp(-m))^2*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(epsilon*(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3)))+1

(2)

subs(epsilon = .17882484, %);

x^2-x*(exp(I*k*`Δx`)+exp(-I*k*`Δx`)-(1-exp(-m))^2+11.18412856*(1-exp(-m))^2*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3))+1

(3)

subs(k = n*Pi, %);

x^2-x*(exp(I*n*Pi*`Δx`)+exp(-I*n*Pi*`Δx`)-(1-exp(-m))^2+11.18412856*(1-exp(-m))^2*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3))+1

(4)

subs(`Δx` = m, %);

x^2-x*(exp(I*n*Pi*m)+exp(-I*n*Pi*m)-(1-exp(-m))^2+11.18412856*(1-exp(-m))^2*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3))+1

(5)

subs(m = 0.1e-2, %);

x^2-x*(exp((0.1e-2*I)*n*Pi)+exp(-(0.1e-2*I)*n*Pi)-(1-exp(-0.1e-2))^2+11.18412856*(1-exp(-0.1e-2))^2*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3))+1

(6)

j := subs(n = 1, %);

x^2-x*(exp((0.1e-2*I)*Pi)+exp(-(0.1e-2*I)*Pi)-0.9990006498e-6+0.1117295170e-4*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3))+1

(7)

complexplot3d(x*j, x = -2-I .. 2+I);````

complexplot3d(x*(x^2-x*(exp((0.1e-2*I)*Pi)+exp(-(0.1e-2*I)*Pi)-0.9990006498e-6+0.1117295170e-4*(-cosh(2)-4*sinh(2)+cosh(2)^3)^2/(-2*cosh(2)-8*sinh(2)+2*cosh(2)^3))+1), x = -2-I .. 2+I)

(8)

a := fsolve(x*j, x);

0., .9865070072, 1.013677544

(9)

b := fsolve(x*j, x = 1);

0., .9865070072, 1.013677544

(10)

with(plots):

complexplot({a, b}, numpoints = 100, color = green, filled = true, title = "Stability Region");

 

``

``

``


Download Stability_1.mw

i want to find the stability of this equation, but there is seem to have some problems..can somebody help me..

 

y := A*(1/x+x*exp(-2*sqrt(-1)*b))+4*(exp(h)-1)^2*(2*exp(-sqrt(-1)*b)-3*(exp(h)-1)^2*x^(-exp(h)+1)*exp(sqrt(-1)*b*(-exp(h)+1-1))+3*x^(-exp(h)+1)*exp(sqrt(-1)*b*(-exp(h)+1-1)))/(3*(1-r))-exp(-2*sqrt(-1)*b)/x-x;

subs(A = (1+r)/(1-r), %);

subs(r = (1/3)*(exp(h)-1)^2, %);

subs(b = m*Pi*(exp(h)-1), %);

subs(m = 1, %);

subs(h = 0.5e-1, %);

> ans := solve(%, x);
Warning, solutions may have been lost
ans:=

> r1 := ans[1];
Error, invalid subscript selector
> r2 := ans[2];
Error, invalid subscript selector

Hye, have a good day everyone..

can anyone know why I got this error? :(second_grade.mw

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system


Here I have attached the file. thank you in advance :)

This might be a ridiculous question, but what is the easiest way to check the monotonicity of a function? I didn´t find anything helpful in the Maple Help, so maybe one of you could give me a nudge here :)

I find mapleform software   but I dont know how to use , and does anyone know another method to convert.  I know use 

> with(MmaTranslator);
> MmaToMaple();  

I can  automatically translated .nb files.  but this have  

Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist

another Error, missing operator or `;`

also I try to use 

> with(MmaTranslator);
> FromMmaNotebook(Mma_notebook_filename, options);

I still dont know how is works?  can you explains for me and show me some example  , here is my example files.plot.zip

 

Hello everyone!

I'm trying to understand how to work with PDEs in Maple by solving a basic PDE such as the wave equation, along with conditions at t = 0 in terms of arbitrary functions. I stumbled upon the very same example, included in the Maple Online help at this page under the part "Three Textbook Examples":

I tried to copy and paste the code and execute it, but to no avail: it seems like pdsolve doesn't give me any other result than a blank space, unlike in the example reported in the Online help. Can anyone tell me if (and where) I am doing wrong? Is it possible that my version of Maple (Maple 13) may not be compatible with the procedures required by this example?

 A warm thank you to whoever will help me!

Hello,

I'm trying to calculate and export a huge "Strategy-Matrix". The calculation works in a smaller version with less possible strategies. For this version I use a specialised simulation-PC with enough memory (that was my problem using my own PC). But now I've got this error massage ("Kernel connection has been lost. You should save this worksheet and restart Maple. Executing commands in Maple requires a connection to the Maple kernel. Firewalls have been known to cause problems with kernel connections in Maple. Please ensure that any firewall software is configured to allow Maple to create connections to the kernel. Consult the FAQ for more information.")

The firewall should not be the problem. Below the code. num_Strat is 63001. 

Thank you in advance for your help.

 

m_Gi:= Matrix(num_Strat, num_Strat):

for t from 1 by 1 to num_Strat
do
Digits:= 5:
v_h2:= evalf((-(24/17)* v_pEi[t]^2)*v_h1 + (v_pEi[t]*(-(128/51)*v_pKi[t]+980/51))*v_h1+ v_pEi[t]*((40/51)*v_pKj+(15/17)*v_pEj)-((110/51)*v_pKi[t]^2)*v_h1 + v_pKi[t]*((55/51)*v_pKj+1220/51*v_h1+(10/17)*v_pEj)):
m_Gi(..,t):= v_h2:
end do:
ExportMatrix(TestMatrixNK, m_Gi):

How do we print "Hello" in green in a .mw file in Maple 13?

Thanks!

mapleatha

I cannot make Maple 13 read the maple.ini file no matter where it is located.

I am using Windows 7, 64 bit.

I thought that ...maple 13\bin.win would do it.

The file contains a few commands of mine. Three of these commands

"read" three procedure files A.m, B.m C.m. respectively, and the other two

are

with(LinearAlgebra):
with(DEtools):

Could somebody up here direct me properly about this?

I would appreciate it.

Thank you!

mapleatha

 

There is no menu item called "Startup Code"

in the "Edit" menu item of my Maple 13.

Please help!

 

Thanks!

I have a surfdata-plot which - by interpolation - goes through different circles with different inclinations with respect to each other. The surface is color-coded with zhue from violet to red. In order to see the valleys behind peaks I would like to have a filled 2D-plot of my zhue-surfdata colors with contours wherever the surface passes z-coordinate integers located somewhere below the surface (e.g. at z=-7). Unfortunately, I haven't found examples of contourplots for surfdata plots so far.

restart

with(plots);

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, graphplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

(1)

with(CurveFitting);

[ArrayInterpolation, BSpline, BSplineCurve, Interactive, LeastSquares, PolynomialInterpolation, RationalInterpolation, Spline, ThieleInterpolation]

(2)

R := [30, 32.5, 37.5, 42.5, 47.5, 50]:

incl := [0, .5, 1, -2, 5, 0]:

phases := [0, (1/12)*Pi, (1/8)*Pi, (1/4)*Pi, (1/2)*Pi, 0]:

colors := [grey, black, blue, red, green, grey]:

orbit := [R[j]*cos(t), R[j]*sin(t), incl[j]*cos(t-phases[j])]:

display(seq(spacecurve(orbit, t = 0 .. 2*Pi, color = colors[j], view = [-56 .. 56, -56 .. 56, -15 .. 15], labels = ["x [AU]", "y [AU]", "z [AU]"]), j = 1 .. 6));

 

NULL

R := [30, 32.5, 37.5, 42.5, 47.5, 50]:

incl := [0, .5, 1, -2, 2*evalf(Pi)-1.75, 0]:

pointplot(R, incl, color = [grey, red, blue, green, black, grey], labels = ["radius", "incl"]);

 

NewR := [seq(30+.4*i, i = 0 .. 50)]:

Newincl := ArrayInterpolation(R, incl, NewR, method = spline):

pointplot(NewR, Newincl, labels = ["radius", "incl"]);

 

phases := [0, (1/12)*evalf(Pi)+.3, 3*evalf(Pi)*(1/4), (1/4)*evalf(Pi)+.2, evalf(Pi)/(2.5)+.5, 0]:

pointplot(R, phases, labels = ["radius", "phase"]);

 

Newphases := ArrayInterpolation(R, phases, NewR, method = spline):

pointplot(NewR, Newphases, labels = ["radius", "phase"]);

 

t := [seq(0+i*(2*Pi*(1/50)), i = 0 .. 50)]:

f := proc (i, j) options operator, arrow; [NewR[i]*cos(t[j]), NewR[i]*sin(t[j]), Newincl[i]*cos(t[j]-Newphases[i])] end proc;

proc (i, j) options operator, arrow; [NewR[i]*cos(t[j]), NewR[i]*sin(t[j]), Newincl[i]*cos(t[j]-Newphases[i])] end proc

(3)

Surface := [seq([seq([NewR[i]*cos(t[j]), NewR[i]*sin(t[j]), Newincl[i]*cos(t[j]-Newphases[i])], i = 1 .. 51)], j = 1 .. 51)]:

plots[surfdata](Surface, labels = ["x [AU]", "y [AU]", "z [AU]"]);

 

``


Download Worksheet_Surfdata-contourplot.mwWorksheet_Surfdata-contourplot.mw

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