hello everyone, I have problem in my project written by matlab to solving my equations with rang-kuta 4 order method. I send my project file.  please help me to run it correctly. thank you.

input Q_j n j V Z F_ve a_1 a_2 p_1 p_2 L A f m k c h t i;

for jj=1:n;

            s=Q_j(jj);

end

b=(-V)*Z;

o=(F_ve*A)*(((a_1^2)*p_1)-((a_2^2)*p_2))/L;

% RK4 method;

% md2z/dt2+cdz/dt+kz=f;

% dz_1/dt=z_2;

% mdz_2/dt+cdz_1/dt+z_1k=f;

% dz_2/dt=1/m(f-cz_2-kz_1);

f=0;

m=1;

k=s+b+o;

c=0;

h=0.5;

t=0:0.5:t;

z1=zeros (1,n);

z2=zeros (1,n);

z1 (1)=0;

z2 (1)=1;

% N=length (t);

f1=@(t,z1,z2)z2;

f2=@(t,z1,z2)f-(c*z2)-(k*z1);

for i=1:n-1;

    k1=f1 (t (i),z1 (i),z2 (i));

    m1=f2 (t (i),z1 (i),z2 (i));

    k2=f1 (t (i)+h/2,z1 (i)+0.5*k1*h,z2 (i)+0.5*m1*h);

    m2=f2 (t (i)+h/2,z1 (i)+0.5*k1*h,z2 (i)+0.5*m1*h);

    k3=f1 (t (i)+h/2,z1 (i)+0.5*k2*h,z2 (i)+0.5*m2*h);

    m3=f2 (t (i)+h/2,z1 (i)+0.5*k2*h,z2 (i)+0.5*m2*h);

    k4=f1 (t (i)+h,z1 (i)*k3*h,z2 (i)*m3*h);

    m4=f2 (t (i)+h,z1 (i)*k3*h,z2 (i)*m3*h);

    z1 (i+1)=z1 (i)+(h/6)*(k1+(2*k2)+(2*k3)+k4);

    z2 (i+1)=z2 (i)+(h/6)*(m1+(2*m2)+(2*m3)+m4);

end

plot (t,z1)

plot (t,z2)

plot (z,t)

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

Download Stability_1.mw

Hi everyone!

I use two "ArrayTools:-Copy" commands to copy elements from Vector A to Vector B.

restart;
A := LinearAlgebra:-RandomVector(10); A := convert(A, Vector[row]);
B := Vector[row](10); ArrayTools:-Copy(2, A, 1, B, 1); ArrayTools:-Copy(2, A, 7, B, 7); B;

My question is: Could I obtain the same result using only one command?

P.S. I tried to use "ArrayTools:-BlockCopy", but didn't get correct result.

I want to retrieve all the commands input by user in a maple document saved on disk. How this can be done? DocumentTools:-Retrieve permits to receive only the labeled expressions but not the input commands.

The following error occurred when I simulate a build-in model, anyone could help me to solve this problem? Thanks first

Hi
how can we write a same result n-times near each other


x:=3:
for i from 1 to 4 do       # here the 4 can be a variable and can change...
    x;
    od;
                                  3
                                  3
                                  3
                                  3

I need the result to come in that form  :    3  ,  3  ,  3  ,  3
and because the number of times can change i can't just easily write x,x,x,x;

thanks in advance

guys, i put a test solution in my differential equation and when i excuted it, maple calculated the solution in term of functions itself. 

i mean, i have functions A(r),B(r),C(r),phi(r) in a differential equation and then i put A(r)=(1-a/r)^alpha and B(r)=(1-a/r)^beta and C(r)=(1-a/r)^gamma and phi(r)=Q*(1-a/r)^lambda , maple give me result with A and B and ... . i want to maple give the result with out any function and just in term of (1-a/r)^... and alpha , betta , .... and c

and is there any way to arrange in term of power of (1-a/r)^... or something like that ?

what is your opnion ?

and also i obtained the differential equation by maple and it is not simplify completely and some terms can cancel each other how can i simplify it completely?

thanks in advance

de.mw

worksheet/expressions/copypasteMaple

Gerschgorin := proc (A::Matrix) local Delta, m, n, AA, R, C, i, c, eig, P, Plt; Delta := proc (i, j) if i = j then 0 else 1 end if end proc; m, n := LinearAlgebra[Dimension](A); AA := Matrix(m, n, proc (i, j) options operator, arrow; Delta(i, j)*abs(A[i, j]) end proc); R := evalm(`&*`(AA, Vector(m, 1))); C := {seq(('plottools[circle]')([Re(A[i, i]), Im(A[i, i])], R[i], color = violet), i = 1 .. m)}; c := {seq(('plottools[point]')([Re(A[i, i]), Im(A[i, i])], color = blue, symbol = diamond), i = 1 .. m)}; eig := evalf(LinearAlgebra[Eigenvalues](A)); P := {seq(('plottools[point]')([Re(eig[i]), Im(eig[i])], color = red, symbol = box), i = 1 .. m)}; Plt := `union`(`union`(C, c), P); plots[display](eval(Plt), scaling = constrained) end proc

A := Matrix([[5, 8, 4, -3], [8, -9, 7, 5], [0, 4, 4, 2], [5, -5, 9, -9]]); evalf(LinearAlgebra[Eigenvalues](A), 3); Gerschgorin(A)

worksheet/expressions/pasteMathML

F := Matrix([[2, -1/2, -1/3, 0], [0, 6, 1, 0], [1/3, -1/3, 5, 1/3], [-1/2, 1/4, -1/4, 4]]); evalf(LinearAlgebra[Eigenvalues](F)); Gerschgorin(F)

Could you print A & F ?

regards

Hello i would love some help with this code i need Maple to make an exercise guide in word automaticly the thing is that i don't know how to automaticly insert the plot in the word, i think there is no fprintf format to add images or plots maybe with exportplot? i really don't know; here is the Code is Basic for a starter student i just need the comand to insert the plot just that.
-------------------------------------------------------------------
fd:fopen("Graphs.doc", WRITE);
 for i from 1 by 5 to 11do
  if i=1 then plotsetup(bmp,plotoutput=plot1);
 a:= -50;
 b:=30;
 elif i=6 then
 plotsetup(bmp,plotoutput=plot2);
 a:= 50;
 b:=-30;
  elif i=11 then  plotsetup(bmp,plotoutput=plot3);
  a:= 2;
 b:=-3;
 end if;
 h:=(i*x^(2)-((3*a)/(i-42))*x-(i*b)/(-i+42));
  fprintf(fd,"the quadratic equation is: ")
fprintf(fd,"%q \\n",h);
 #`HERE IS WHERE I NEED THE PLOT TO BE`
  plot(i*x^(2)-((3*a)/(i-42))*x-(i*b)/(-i+42),x)
   end do; 
 fclose(fd);
--------------------------------------------------------------------
pd. i can't do it manually i need it to be automatic later i will add some inputs to make it even faster.

Hi, as I can't manage to copy and paste on mapleprimes, I would be glad to get a hint ...

Hello,
I am looking here for a tutor, that can help me doing some Maple V (5) programming.
I am a mathematic student and we use Maple programming.
I will be happy to pay a small amount of money for each exercice you help me doing,

if anyone is interested, please contact me here.

(The procedures that we usually have to write are for example:
 Newton-Raphson Method, Chebyshev Polynomial,...  I don't think it is hard for you.
Thank you very time for your time and your help.

How to find the nth derivative of (logx)/x  and (e^x)logx by using leibenitz theorem....?

Hello,

I would like to export maple inputs as xml code. The aim is to recover the semantics of the inputs.

I have tried without success with the Export Function.

Have you some ideas to export maple inputs as xml and with preserving the semantics of the input ?

Thanks a lot for your help

Maple function plot allows to build a linear plot for two vectors, when one vector contains X coordinates, another -- Y coordinates (similar to plot below):

Is it possible to build such a plot for several Y-vectors like this:

I want to solve the following differential equation

solve(y''(x)=(λ*x* y[x])/Sqrt(-1+ x), y(x),x)

But do not know how to actually solve it. Any suggestion?