Hi, I'm starting my program using Maple, then will bring some data to calculate in matlab, later on bring the data back to be calculated in Maple.

My equations have lots of trigonometry multiplication of matrix. When I run it in matlab, it runs fine. But when I try to use Matlab[FromMFile]("test.m",'evaluate') it gives me this statement:

"Error, (in Matlab:-FromMFile) unable to strore '(3.500000000*I)*2^(1/2)' when datatype=float(8)"

Could someone please...

Hello.

I took the battery out of my laptop by a mistake while it was in sleep mode. Now maple can't open the document even though i did save it before. It opens a file where 80% of my work is gone.

Download BET.mw

Maple says Incomplete worksheet when I open the file.  

Can someone please help me?

Hi all,

I am trying to get the following matrix:

restart:
int00:=Matrix(3,3);
int00[2,2]:=a;
int00[2,3]:=b;
int00[3,2]:=c;
int00[3,3]:=d;
int00;
cols:=MTM:-sum(int00,1);
rows:=MTM:-sum(int00,2);
int00[1,1..3]:=-cols:
int00[1..3,1]:=-rows:
int00;

int00[1,1]:=a+b+c+d:
int00;
MTM:-sum(int00,1);
MTM:-sum(int00,2);

How do I get the sum of all elemnets in the 'original' matrix? Bascially, the int00[1,1] term in the toy example.

Z := proc (theta) options operator, arrow; sin(theta+2.3)+sin(2*theta+2*(-1.5))+sin(3*theta+3*(-0.94e-1)) end proc

Zbar := proc (theta, phi) options operator, arrow; Z(theta+phi)-Z(theta)end proc

plot3d(int(signum(Zbar(theta, phi)+l)*abs(Zbar(theta, phi)+l)^(1/4), theta = 0 .. 2*Pi, phi = 0 .. 1, l = 0 .. 1)

Hello, I am having some trouble plotting 3D of these function. There is no error message or anything, but when I enter plot3d command it takes too...

I want to find all the max or min values in a range, from a numerical function [t,x(t)], numerical solution of a differential equation.

For example if the solution was periodic in time, I want to find the max values of x in a certain intervalot t

Hi Guys. I´ve got a problem with the mod (mods,modp) function of maple.

I want to calculate the following congruences:

mod(4,-5) = -1 BUT maple returns 4

mod(-3,-4)=-3 but maple returns 1

and so on....

Does anybody of you know how to fix this? thanks!

Hello,

I have a problem. I want to make plot 3d for pressure drop in diffuser by posuille formula,

but with increasing of "l"(length of channel) graph shows wrong pressure drop. Normally with i increasing of "L" we should have...

> restart; R := 8.314;print(??); # input placeholder> G1[1] := piecewise(`and`(T > 298, T < 700), -7976.15+137.093038*T-24.3671976*T*log(T)-0.1884662e-2*T^2-0.877664e-6*T^3+74092/T, `and`(T >= 700, T < 933), -11276.24+223.048446*T-38.5844296*T*log(T)+0.18531982e-1*T^2-0.5764227e-5*T^3+74092/T, `and`(T >= 933, T < 2900), -11278.361+188.684136*T-31.748192*T*log(T)-0.1230622e29/T^9, 0);print(??); # input placeholder

> G1[2] := piecewise(`and`(T > 298, T < 700...

The system consists of 3 equationsHow can you single out one of the solutions in obtaining module and plot for it?
a := 1; m := 1/64 
l1 := 2*(diff(ur(e, t), t))+(1+tanh((1/20)*e))*(diff(ui(e, t), `$`(e, 2)))+2*a*(ui(e, t)*ur(e, t)*ur(e, t)+(ui(e, t)*ui(e, t))*ui(e, t))+m*ui(e, t)*(diff(ur(e, t)^2, e))+m*ui(e, t)*(diff(ui(e, t)*ui(e, t), e)) = 0 
l2 := -2*(diff(ui(e, t), t))+(1+tanh((1/20)*e))*(diff(ur(e, t), `$`(e, 2)))+2*a*(ur(e, t...

Let say that you have two datasets:

restart:
with(MyFinance):
z1 := GetCSV("http://www.quandl.com/api/v1/datasets/USER_18Q/18R.csv?&auth_token=BXckVSG3twWkV272c9Q4&trim_start=1968-04-01&trim_end=2013-04-12&collapse=daily&transformation=rdiff&sort_order=desc");
z2 := GetCSV("http://www.quandl.com/api/v1/datasets/USER_18Q/1EV.csv?&auth_token=BXckVSG3twWkV272c9Q4&trim_start=1994-01-03&trim_end=2013-04-12&collapse=daily&...

** comparison link has been updated for Maple 2015 **  All three M's have identical speed performances

Someone has created a comparison of Mathermatica 9, Maple 17 and Matlab 2013a here http://12000.org/my_notes/rankTest/test.htm

The timings and codes used were as follows:

2.2.1 Mathematica

Remove["Global‘*"];  
Share[];  
n = 7000;  
m = Table[Random[], {i, 1, n}, {j, 1, n}];  
AbsoluteTiming[MatrixRank[m];]


2.2.2 Maple

restart;  
kernelopts(gcfreq= 2^22):  
UseHardwareFloats:= true:  
gc():  
 
n:=7000:  
 
M:= LinearAlgebra:-RandomMatrix(  
     n  
    ,n  
    ,generator=0.0 .. 1  
    ,outputoptions=[datatype=float[8]]  
):  
 
time[real](LinearAlgebra:-LA_Main:-Rank(M));



2.2.3 Matlab

clear all;  
n=2500;  
A=rand(n,n);  
tic();  
rank(A);  
toc()

Certainly that must be the slowest possible route for Maple ever.  Anyone care to post their code and timings?

Maple gives the fourier transformation:

fourier(int(y^k*exp(-k*x*y/phi+I*y*omega), y = 0 .. infinity), omega, t) as:

2*Pi*t^k*exp(-t*k*x/phi)*piecewise(t < 0, 0, 1)+piecewise(0 < t, -2*Pi*t^k*exp(-t*k*x/phi), 0)

but it is identically 0.

if change the order of integration and integate exp(I*y*omega) first, the transform is equal to

int(Dirac(y-t)*y^k*exp(-k*x*y/phi), y = 0 .. infinity)

which is Heaviside(t)*t^k*exp(-t*k*x/phi)

When I solve a systems of equations, fsolve gives the solution

.

I now want want to define the variables x, y and z as respectively 5, 6 and 1 - but without copying the solution. Because this is a solution that I get in a loop, so it changes in every step. Maybe there is a easy way to put it in a vector?

I often have large matrices and have to grab only every nth column or row (for plotting

for example). So far I've been doing it this way:

nthColumns:=proc(v_m, v_n)  

v_m[..,[seq(i, i=1..rtable_size(v_m)[2], v_n)]]

end;

nthRows:=proc(v_m, v_n) 

 v_m[[seq(i, i=1..rtable_size(v_m)[1], v_n)],..]

end;

A := LinearAlgebra:-RandomMatrix(20, 20, outputoptions = [datatype = float[8]]);

I only want to display 'position' in animation while tend <= 12:
This is the trial:

restart:with(plots):
for tend from 0 to 50
while (tend <= 12) do
pos:=animate(textplot3d,[[-1.8,0,alpha*tend, 'position']], tend=0..50,frames=150);
od;
display(pos);

I receive error messages:
error, (in plots/animate) 0=0..50 is an invalid animation range

Then I try:

restart:with(plots):
while (tend <= 12) do