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?

How I can highlight p(1,1) and Q(2,4) on the curve f(x)=x^2?

Hello,

I´ve got 2-dimensional array (100x100) of values, which I'd like to plot as 3d graph, where x,y axis would be coords in the array and z-axis would be the value there.

Values are formed by procedure, not by function, so I can't use standard 3d plot, that plots a 3d surface. I think the best is to use the pointplot3d function, but you have to assign each point separately. 

I tried to create a list of pointplot3d plots, one for each row of the array, and then display these plots all in one, but I just received GLException error.

Do you have any ideas?

Social Golfer Problem

So we have 8 golfers : 2 golfers per group and 4 groups.Over several golfing days, how can they be arranged so no golfer plays with the same player twice.

turns out they can play for 7 days.

So my attempt in maple....golfing.mw

For 3 players per group and 4 groups this is the arrangement over 4 days:

(these was done in wolfram demonstrations project but the code appears to be obfuscated)

have eroror when run command below

animate(dualaxisplot(plot(sin(x+A), x=0..5), plot(cos(x+A), x=0..5), style = line, gridlines = false), A = 0 .. 5);

Dear collagues

Hi,

I write a code to solve a system of ODE. It solve the ODES in a wide range of parameters but as I decrease NBT below 0.5, it doesnt converge. I do my best but I couldn't find the answer. Would you please help me? Thank you

Here is my code and it should be run for 0.1<NBT<10. the value of NBT is input directly in res1.

restart:
EPSILONE:=1000:
Digits:=15:

a[mu1]:=5.45:
b[mu1]:=108.2:
a[k1]:=1.292:
b[k1]:=-11.99:




rhop:=4175:
rhobf:=998.2:
mu1[bf]:=9.93/10000:
k1[bf]:=0.597:

rhost(eta):=1-phi(eta)+phi(eta)*rhop/rhobf;
k:=unapply(k1[bf]*(1+a[k1]*phi(eta)+b[k1]*phi(eta)^2),eta);


eq1:=(diff(u(eta), eta))*a[mu1]*(diff(phi(eta), eta))+2*(diff(u(eta), eta))*b[mu1]*phi(eta)*(diff(phi(eta), eta))+((diff(u(eta), eta))+(diff(u(eta), eta))*a[mu1]*phi(eta)+(diff(u(eta), eta))*b[mu1]*phi(eta)^2)/(eta+EPSILONE)+diff(u(eta), eta, eta)+(diff(u(eta), eta, eta))*a[mu1]*phi(eta)+(diff(u(eta), eta, eta))*b[mu1]*phi(eta)^2+1-phi(eta)+phi(eta)*rhop/rhobf:
eq2:=(1+a[k1]*phi(eta)+b[k1]*phi(eta)^2)*(diff(T(eta), eta))/(eta+EPSILONE) + (a[k1]*(diff(phi(eta), eta))+2*b[k1]*phi(eta)*(diff(phi(eta), eta)))*(diff(T(eta), eta))+(1+a[k1]*phi(eta)+b[k1]*phi(eta)^2)*(diff(T(eta), eta, eta));
eq3:=diff(phi(eta),eta)+phi(eta)/(N_bt*(1+gama1*T(eta))^2)*diff(T(eta),eta):

eq2:=subs(phi(0)=phi0,eq2):
eq3:=subs(phi(0)=phi0,eq3):


eval(dsolve({eq3,phi(1)=phiv},phi(eta)),(T)(1)=1):
Phi:=normal(combine(%)):
Teq:=isolate(eval(eq2,Phi),diff(T(eta),eta)):
ueq1:=eval(eq1,Phi)=0:
ueq2:=subs(Teq,ueq1):


lambda:=0;
Ha:=0;
N_bt:=cc*NBT+(1-cc)*0.8;
kratio:=k1[p]/k1[bf]:






GUESS:=[T(eta) =0.0001*eta, u(eta) =0.1*eta, phi(eta) = 0.3*(eta-1)^4];
res1 := dsolve(subs(NBT=0.48,gama1=0.2,phiv=0.06,{eq1,eq2,eq3,u(0)=lambda*D(u)(0),D(u)(1)=0,T(0)=0,phi(1)=phiv,T(1)=1}), numeric,method=bvp[midrich],maxmesh=4000,approxsoln=GUESS, output=listprocedure,continuation=cc):
G0,G1,G2:=op(subs(subs(res1),[phi(eta),u(eta),diff(T(eta),eta)])):

masst:=evalf(int((1-G0(eta)+G0(eta)*rhop/rhobf)*G1(eta), eta = 0..1));
heatt:=(1+a[k1]*G0(0)+b[k1]*G0(0)^2)*G2(0);

plots:-odeplot(res1,[[eta,T(eta)]],0..1,legend=[T],color=["Black"],linestyle=Solid,axes=boxed,thickness=3);
plots:-odeplot(res1,[[eta,u(eta)]],0..1,legend=[u],color=["Black"],linestyle=Solid,axes=boxed,thickness=3);
plots:-odeplot(res1,[[eta,phi(eta)]],0..1,legend=[phi],color=["Black"],linestyle=Solid,axes=boxed,thickness=3);

>
>

Thank you

Amir

I know this has been dealt with before here, but I have forgotten the proper way to trim a small imaginary round-off from a result. I cannot locate the proper answer in MaplePrimes; and the Mapledocs are either quiet about it or it is hidden in a difficult-to-find place.

What I want is to ignore ny imaginary part below a certain threshold. Ideally, it takes the threshold in relation to the real part but I am not particular about this.

Thanks in advance,

Mac Dude

If for example I type :

It gives out this :

I want dsolve to know I'm using real numbers so it gives out something like :

I tried some assumptions and stuff like assume(a,real), but I didn't manage to figure it out.

If that matters, I'm using Maple 18 student.

EDIT : I know dsolve is not necessary for this particular example, but I want to know if it's possible with dsolve or maybe an other tool that can handle ODE.

Thank you in advance !

 Maple code when open saved maple sheet?

and keep it can run code and show latest result after open

i mean user can not extend or find code or display code in maple sheet , just can run result

Hello all

I am trying to write  a tutorial about systems of linear equations, and I want to demonstrate the idea that when you have a system of 3 euqtions with 3 unknowns, the solution is the intersection point between these planes. Plotting 3 planes in Maple 2015 is fairly easy (you plot one and just drag the others in), but I don't know how to plot the intersection point. Can you help please ?

My equations are:

x-2y+z=0

2y-8z=8

-4x+5y+9z=-9

The intersection point is (29,16,3)

Thank you !

Heun.mw

Download Heun.mw

Hi all

I want to know if there's any plot for the function:

Dear all,

I'm having a bit of trouble figuring out some error messages in the Maple debugger. When I debug some procedures I'm getting error messages like "Error, invalid input: %1 expects its %-2 argument, %3, to be of type %4, but received %5". Since the procedure calls a lot of other procedures I have to check every single one of them, to see, which one gives me the error message. It would be easier if I would know, what procedure/type hides behind the variables %1,... Is there a way to check this from within the debugger?

Thanks a lot.

I want to write a procedure which recieves a natural N and print all the possible 2ˆN sequences of zeros and ones with size N.

 For instance, if I set N = 3; my procedure returns a set with 8 lists of zeros and ones: [[0,0,0],[1,0,0],[0,1,0],[0,0,1][1,0,1],[0,1,1],[1,1,1]] ;

Many thanks!

For those who are running Maple and/or MapleSim on the Mac, and who may have missed an earlier question on this site, we wanted to let you know that there are some problems with running Maple and MapleSim on the new Mac OS, Mac OS X 10.11 (El Capitan). We’re working on a solution, which we expect will be ready in a few weeks. We’ll keep you posted, but in the meantime, please delay updating your Mac OS for now to avoid problems.

eithne

Im working on an assignment involving a rocket being launched vertically from rest. 
I have a eq. of motion that looks like this: v(t)=vs*ln(-m0/(r*t-m0))-g*t

Now Im supposed to implement the air resistance, and write an expression for the total force on the rocket. After that creating a plot displaying the speed v(t) as a function of time for both, with and without air resistance. 

All constants are known including CD, rho and A for the dragforce. 

I would be grateful if anyone could help me with this!

Thanks