# Items tagged with lineline Tagged Items Feed

### Plot of line and point in 3d...

April 06 2016
1 2

hello,
I try to draw together point (0,0,12) and the line  x=4t, y=-2t, z=2t on one graph 3D. Somebody has some pointers? Thank you
Agnieszka

### How to sketch a straight line in three-dimensional...

March 17 2016 Maple
1
2

A wealth of knowledge is on display in MaplePrimes as our contributors share their expertise and step up to answer others’ queries. This post picks out one such response and further elucidates the answers to the posted question. I hope these explanations appeal to those of our readers who might not be familiar with the techniques embedded in the original responses.

The Question: How to sketch a line in space?

vahid65 wanted to know how to sketch the line with equation (x-2)/3 = (y-1)/4 = (z-3)/3 in a three-dimensional space.

This question was answered using two different approaches that we will discuss.

The first approach, given by Preben Alsholm, suggested using these commands: solve, subs, and plots:-spacecurve.

Preben provided the following lines of code:

1. {(x-2)/3 , (y-1)/4 , (z-3)/3} =~t;
2. solve(%,{x,y,z});
3. L:=subs(%,[x,y,z]);
4. plots:-spacecurve(L,t=-5..5,thickness=3,color=red);

The first line sets the three expressions equal to t using the element wise operator =~.  This distributes the operation of equality over the elements of the set of expressions, forming a set of equations.

Result-1: {(x-2)/3=t, (y-1)/4=t, (z-3)/3=t}

The second line invokes the solve command. The solve command, solve(equations, variables), solves one or more equations or inequalities for their unknowns. So in this line, this command was used to solve each expression for its corresponding unknown.

Result-2: {(x= 2+3t, y=1+4t, z= 3+3t}

You may have noticed that the % symbol is used within the command. This symbol is referring to the equation label that was created chronologically last.

The third line uses the subs command. The subs command, subs(x=a,expr), substitutes “a” for “x” in the expression expr. In this case, “a” is the reference to the equation label referenced by (%) and expr is the sequence of variables x,y,z. The square brackets around expr forms a list, and the replacement imposed by the subs command replaces each of the three names x, y, and z with their equivalents in the set of equations returned by the solve command.

Result-3: [2+3t,1+4t, 3+3t ]

Finally the last line uses the plots:-spacecurve function.  This function, spacecurve(sc,r,opts),  graphs a parametrically defined curve in three-dimensional Cartesian space.

In this example the first parameter, sc, is replaced by L since L has the properties of a list. The parameter  r is an equation containing the parameter name and the parameter range; here it is the equation t=-5..5. Last but not least, the opts parameter (which is optional) gives the user the opportunity to specify the optional properties of the graph, properties such as thickness and color of the space curve.

Result 4:

Another contributer, Carl Love, suggested that these commands could be combined using the zip function, zip( f, u, v), which is an unusual application of it. The zip command applies elementwise to corresponding members of the two lists u and v, the binary function f, creating a new list r, whose members are f(uk,vk).

Carl provided the following code:

plots:-spacecurve(  zip(solve, [(x-2)/3 , (y-1)/4 , (z-3)/3] =~ t, [x,y,z]),  t= -5..5, thickness= 3, color= red, axes= normal);

In this case zip is applying solve to u, a list of equations, and v, a list of variables. The equations in u are those in Result-1. The zip function returns the solution of each equation for the indicated variable, that is, the new list r is the list in Result-3. This list is the first argument of plots:-spacecurve. Finally, the parameter range and opts are up to the user’s preferences.

A second approach was suggested by Kitonum. The difference between this approach and the first one is that here we are not using the solve command. Instead, we are treating the line as an intersection of two planes.

Kitonum's code is as follows:

1. L := [(x-2)*(1/3), (y-1)*(1/4), (z-3)*(1/3)]:
2. plots[intersectplot](L[1]-L[2], L[2]-L[3], x = -4 .. 4, y = -4 .. 4, z = -4 .. 4, linestyle = 1, thickness = 3, axes = normal, orientation = [60, 75], view = [-1 .. 3, -4 .. 2, -1 .. 4]);

The first statement initializes a list L that contains all the desired expressions in the given symmetric form of the line.

In the second statement Kionum used the plots[intersectplot] command. This command graphs, in three-dimensional Cartesian space, the intersection of a pair of surfaces. This example uses the calling sequence: intersectplot(expr1, expr2, x=a..b, y=c..d, z= e..f, options)

What are expr1 and expr2 used in Kitonum’s solution? The pairwise differences of L[1] and L[2], and L[2] and L[3] are planes, the intersection of which is the desired line in space.

expr1 = L[1] – L[2] = (x-2)*(1/3)- (y-1)*(1/4)

expr2= L[2] - L[3] = (y-1)*(1/4)-(z-3)*(1/3)

The variables x,y, and z are the coordinate names and define the axes, whereas a, b, c, d, e, f define the ranges. In this case they were taken as -4 or 4.

In the options parameter Kitonum used linestyle, thickness, axes, orientation, and  view. The help page “plot3d,option” details these and other options that apply.

This blog was written by our intern Pia under the supervision of Dr. Robert Lopez. We both hope that you find this useful. If there is a particular question on MaplePrimes that you would like further explained, please let us know.

### How to sketch a line in space?...

February 22 2016
0 4

How I can sketch a line with equation (x-2)/3 = (y-1)/4 = (z-3)/3 in three dimentional space?

### Break to form new lines...

November 11 2015
1 5

Working on a code to create random lines. Here is the code:

M:=1; N:=1500; R3:=rand(1..3): # M=lines, N=steps

for i from 1 to M do
X[i,0]:=0;
Y[i,0]:=0;
for j from 0 to N do
if j=1 then X[i,1]:=1; Y[i,1]:=0;
elif j>1 then
r:=R3();
if r=1 then X[i,j]:=2*X[i,j-1]-X[i,j-2]; Y[i,j]:=2*Y[i,j-1]-Y[i,j-2]
elif r=2 then X[i,j]:=X[i,j-1]+Y[i,j-1]-Y[i,j-2]; Y[i,j]:=Y[i,j-1]-X[i,j-1]+X[i,j-2];
elif r=3 then X[i,j]:=X[i,j-1]-Y[i,j-1]+Y[i,j-2]; Y[i,j]:=Y[i,j-1]+X[i,j-1]-X[i,j-2];
end if;
end if;
R[i,j]:=[X[i,j],Y[i,j]];
K[i,j]:=[X[i,j],Y[i,j]];
end do:
end do:

Now the code works fine. But I don't want the lines to cross its own path. I want it to stop if R[i,j]=R[i,k] for any 0<k<i-2 for j=0..N. What I have so far is:

for i from 1 to M do
for j from 1 to 40 do
R[i,j];
for k from 1 to j-2 do
if R[i,j]=K[i,k] then x:=j; print(b[i]=j); break; end if;
L[i,j]:=R[i,j];
end do;
end do;
end do;

Now it works aswell, but it doesn't stop when it hits the crossing point. I want it to stop so i get one b[i]=j per M. So I need it to break two for loops so that it goes through i M times and stops everytime it reaches a crossing point

### Secant line equation...

October 25 2015
0 2

How I can write a procedure in Maple that return a equation of line that passes of two points (a,f(a)) and (b,f(b)) on the curve y=f(x) and the final plot y=f(x) and secant line in the same screen?

### Line appears like a plane - geom3d...

September 05 2015
0 4

Hello!

I defined a line normally, but when I try plot this line in R3, appears a shape of a plane.

How can I plot only a line in a 3D window?

Line has coordinates x,y, e z.

### Method of line(MOL)?...

May 01 2015
0 5

Hi

Do the following equation solve by method of line(MOL) with this software?

du/dt=0.01 d^2u/dx^2

conditions:

u(x,0)=10

u(0,t)=0

u(1,t)=0

t=1,delta t=0.1,x=1,Nx=10

thank you

### How to plot a vertical line in Maple 15?...

November 29 2014
0 6

How can I plot a vertical line say x = a in Maple 15 worksheet?

Thanks

### Line style and legend for graph...

November 15 2014
1 3

AoA... I want to plot the following functions

f(x)=sin(x)

g(x)=sin(x)-x+x^2

h(x)=cos(x)-sin(x)

I(x)=exp(x)-x

J(x)=sin(x)+exp(x)

in one coordinate having the line style and legend like attached file.

1-s2.0-S0377042714003331-main.pdf

PhD (Scholar)
Department of Mathematics

### how to add another line in this template of formul...

October 19 2014
1 4

how to add another line in this template of formula

use the expression template

how to add another line like

### Numerical Integral by trapizidal rule...

September 24 2014
0 9

Hello.

I have a line integral and i want to calculate and graph this line integral.

The line integral comes from a physical problem, the integrand is good in the interval [0, 2*Pi] and the value of the integrand at zero = the value of the integrand at 2*Pi.

The problem is when i use the trapizidal rule to calculate the line integral.

The results is so bad, i.e.  the value of the line integral at zero not equal to the value of the line integral at 2* Pi.

where is the problem.

Any one could help me.

### why F3 (split execution group) not smart enough to...

August 21 2014
0 3

One would expect that hitting F3 will split from _end_ of the current line (where the cursor is at).

But what Maple does is actually split everything from the current cursor location. Which means if the cursor happened not to be exactly at the end of the line, the current line itself will also be split and broken.

It is much more logical to split starting from end of current _line_ (where the cursor is at), not current character, because that is what normally one would want to do. One will have large block of code, and want to split it from one line down to the end.  It is a simple usability issue, which Maple UI seems to suffer allot from.

Is there a way to modify this behavior? I keep hitting F3 and forget to move the cursor to the end of the line before, and end up wasting time having to fix things afterwords since the line itself is split.

### Error, (in f) unable to store .......................

May 27 2014
0 2

Hello Hello everybody
I have to solve the following differential equation numerically

 > restart:with(plots):
 > mb:=765 : mp:=587 : Ib:=76.3*10^3 : Ip:=7.3*10^3 : l:=0.92 : d:=10: F:=490: omega:=0.43 :
 > eq1:=(mp+mb)*diff(x(t),t\$2)+mp*(d*cos(theta(t))+l*cos(alpha(t)+theta(t)))*diff(theta(t),t\$2)+mp*l*cos(alpha(t)+theta(t))*diff(alpha(t),t\$2)+mp*(d*diff(theta(t),t)^2*sin(theta(t))+l*(diff(theta(t),t)+diff(alpha(t),t))^2*sin(alpha(t)+theta(t)))-F*sin(omega*t)=0;
 (1)
 > eq2:=(mp+mb)*diff(z(t),t\$2)-mp*d*(sin(theta(t)+alpha(t))+sin(theta(t)))*diff(theta(t),t\$2)-mp*l*sin(alpha(t)+theta(t))*diff(alpha(t),t\$2)+mp*(d*diff(theta(t),t)^2*cos(theta(t))+l*(diff(theta(t),t)+diff(alpha(t),t))^2*cos(alpha(t)+theta(t)))+9.81*(mp+mb)-F*sin(omega*t)=0;
 (2)
 > eq3:=mp*(d*cos(theta(t))+l*cos(alpha(t)+theta(t)))*diff(x(t),t\$2)-mp*(l*sin(theta(t)+alpha(t))+d*sin(theta(t)))*diff(z(t),t\$2)+(Ip+Ib+mp*(d^2+l^2)+2*mp*d*l*cos(alpha(t)))*diff(theta(t),t\$2)+[Ip+mp*l^2+mp*d*l*cos(alpha(t))]*diff(alpha(t),t\$2)-mp*sin(alpha(t))*(l*d*diff(alpha(t),t)^2-l*d*(diff(alpha(t),t)+diff(theta(t),t))^2)+mp*9.81*l*sin(alpha(t)+theta(t))+mp*9.81*d*sin(theta(t))=0;
 (3)
 > eq4:=mp*l*cos(alpha(t)+theta(t))*diff(x(t),t\$2)-mp*l*sin(alpha(t)+theta(t))*diff(z(t),t\$2)+(Ip+mp*l^2+mp*d*l*cos(alpha(t)))*diff(theta(t),t\$2)+(Ip+mp*l^2)*diff(alpha(t),t\$2)-mp*9.81*l*sin(alpha(t)+theta(t))+l*d*mp*diff(theta(t),t\$1)^2*sin(alpha(t))=0;
 (4)
 > CI:= x(0)=0,z(0)=0,theta(0)=0,alpha(0)=0,D(x)(0)=0,D(alpha)(0)=0,D(z)(0)=0,D(theta)(0)=0;
 (5)
 > solution:=dsolve([eq1,eq2,eq3,eq4, CI],numeric);
 Error, (in f) unable to store '[0.]/(0.17571268341557e16+[-0.25659510610770e15])' when datatype=float[8]

I don't know why it says : Error, (in f) unable to store '[0.]/(0.17571268341557e16+[-0.25659510610770e15])' when datatype=float[8]

Help pleaase!

thank you !!!

### printf and print output on the same line...

March 20 2014
0 3

Hi, i'm wondering if theres anything i can do printf and print and have the outputs on the same line. I want the string wihtout any quotations so i think i can only do printf to get this, but then i'm putting an entry out of an array after this. This is a small example out of my code, i would like the output to be all on one line:

The synchronising vector is (S[j-1]) The synchronising word is (X[j-1]) The length of the synchronising word is ...

if (sync=true) then
printf("The Synchronising vector is"); print(S[j-1]);
printf("The sychronising word is");
print(X[j-1]);
printf("The length of the synchronising word is");
print(length(X[j-1]));
else printf("There is no synchronising word");
end if;

At the moment its outputing each print on a new line.

Any help would be much appreciated, thanks.

### getting parametric equations...

February 23 2014
0 4

Hello Maple users friends,

I have two lines in the space (x,y,z) described by the equations in L1 and L2:

L1:= {4*x + 3*y + z = 0, x + y - z - 15 = 0}:

L2:={12*x + 5*y + 7*z -13 = 0, 9*x + y -3*z - 5 = 0}:

I would like the get the parametric (with z=t) equations P1 and P2 of the two lines..

I see the "form" of such parametric equations P1 and P2 using "solve"

solve(L1, {x, y}); solve(L2, {x,y});

but I do not know how to use those values to get my parametric equations P1 and P2 to continue with additional computation (area, volume etc).

Thanks for your attention and help.

JJ

 1 2 3 Page 1 of 3
﻿