Hi everyone.

I have been experiencing a problem trying to solve a coupled system of 3 differencial equations

My problem is that a got a message back as I try to solve the system:

"Error, (in pdsolve/numeric/process_IBCs) improper op or subscript selector"

by apply this point that ''all dependent variables must be functions of the same independent variables''

i again accost with another error ''

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions can only contain derivatives which are normal to the boundary, got (D[1, 1](w))(x, -3/400000000)

''

please help me.....very very thanks

Here the potential of maple 2015 to the quantitative study of the decomposition of a vector table is shown in two dimensions. Application for the exclusive use of engineering students, which was implemented with embedded components.

Atte.

Lenin Araujo Castillo

Archivo Corregido:  Decomposición_Vectorial_Corregido.mw

I wish to display the following screen output from a procedure A(), in a textbox of mathcontainer or any other suitable component.

This output to screen listing is generated by print statements in the procedure A().

I have placed the procedure call A() in a textbox and then am using a button component to execute it with the objective of displaying the above output in another component display box; either a textbox or a MathContainer box...  I can see that A() executes, but do not know how to get the result to display in a component box.

Can anyone offer any help?

hi.may  help me for solve this nonlinear equations by numeric solver maple39.d39.pdfocx39.pdf

thanks alot

file format is pdf and word type

and show whether there is sink or source?

ode1:=a(t)*(diff(a(t), t))+c(t)*(diff(c(t), t))=3*t;
ode2:=a(t)*(diff(a(t), t))+a(t)*(diff(a(t), t))*b(t)*(diff(b(t), t))*c(t)*(diff(c(t), t))=3*t;
ode3:=a(t)*(diff(a(t), t))+a(t)*(diff(a(t), t))*c(t)*(diff(c(t), t))+a(t)*(diff(a(t), t))*b(t)*(diff(b(t), t))*c(t)*(diff(c(t), t))=2*t;

DEtools[DEplot3d]({ode1,ode2,ode3},[a(t),b(t),c(t)],t=0.5..1.4,a=0..1.4,b=0..1.4,c=0..1.4,[[a(1)=1,b(1)=2,c(1)=3]],scene=[a(t),b(t),c(t)], arrows=medium);

unable to convert to explicit first order system, do not understand this error message

Hi!!

I'm trying to define a function by a definite integral. For this, first we define the following  procedure

and the map

So, for a given m, my desired  "function" should be fun(curva(m,t)). The problem here, is that this function not work because, for instance, for m=2 the command fun(curva(2,t)) returns the below expession, which i think is wrong (where is the expression on "sin" or "cos(2*Pi*t)" ???)

Somebody can help me??

Thanks!!!

Hello!

Is it possible convert some equations in R² (ellipse, hyperbola) to polar coordinates using Maple? And a regular  object in R³ like a sphere, a cylinder or a cube could be converted from cartesian coordinates to polar coordinates?

Thank you so much.

Hi everybody,

i'm trying to do an elliptic regularization but i don't know how to proceed ?

Is someone know how to achieve to do that  with an example ?

thanks a lot !

PS: i know only how to do a linear regularisation.

I am considering the following PDE and I am getting an error, please suggest a better numerical method than the default one used in maple:

the PDE is:

u_{xx}u^3 - sin(xt)u_{tt} = u(x,t)

u(x, 0) = sin(x), (D[2](u))(x, 0) = cos(x), u(0, t) = cos(t), (D[1](u))(0, t) = sin(t)

Please suggest me a method that will also work for the following PDEs:

u^m* u_{xx} - sin(xt)u_{tt} = u^n

for m,n =0,1,2,3,... for the cases m=n and m not equal n

Here's the code:

Download nonlinear_hyperbolic_PDE.mw

Hi,

I'm trying to compile a worksheet and save it on Mac OSX. The worksheet ends with this code:

>

> save `E:\\class work\\THERMODYNAMICS\\ThermoSoft\\Compiled\\SteamTables.m`;

>

The guide tells me to change the path to a user specified folder where the .m file will be saved.

Since the path name already given is a Windows path, I have to change that to a Mac path.

My problem is, that i don't know what I should write?

I have tried writing:

save `\\Macintosh HD\\Users\\Name\\Folder\\SteamTables.m`;

But it gives an error, and can't save the file. 

LL_104)_NASDAQ.mw
Portfolio_Optimization.txt

Portfolio Optimization with Google Spreadsheet and Maple
 

I will in this post show how to manage data and do portfolio optimization in Maple by using google spreadsheet.

You can either use a direct link to the data:

https://docs.google.com/spreadsheets/d/1L5-yUB0EWeBdJNMdELKBRmBQ1JJ0QymrtDLkVhHCVn8/pub?gid=649021574&single=true&output=csv

or you can set up your own google spreadsheet. If you choice to set up your own spreedsheet follow the below road map:

1) select which market you want to follow:

NASDAQ

http://www.nasdaq.com/screening/companies-by-industry.aspx?exchange=NASDAQ&render=download

NYSE

http://www.nasdaq.com/screening/companies-by-industry.aspx?exchange=NYSE&render=download

AMEX

http://www.nasdaq.com/screening/companies-by-industry.aspx?exchange=AMEX&render=download

2) Create a new google spreadsheet and name two sheets Blad1 and Panel. In the first cell of Blad1 you put the formula:

=IMPORTDATA("http://www.nasdaq.com/screening/companies-by-industry.aspx?exchange=NASDAQ&render=download")

you need to change the url to match your selection in 1).

3) In the first cell of Panel you put the name "Ticker" and then you copy all the ticker names from Blad1.

4) In the script editor you put in the below java script code:

function PanelCreation_Stock() 

{
var ss = SpreadsheetApp.getActiveSpreadsheet();
var sourceSheet = ss.getSheetByName("Blad1");
var dstSheet = ss.getSheetByName("Panel");
var curDat = new Date();
var day1 = curDat.getDay();
if(day1 == 0 || day1 == 1)
{
return;
}
var lCol = dstSheet.getLastColumn();
var srcdate = dstSheet.getRange(1, 1, 1, lCol).getValues();

for(var k=1;k<=srcdate[0].length-1;k++)
{
if(Utilities.formatDate(srcdate[0][k],"GMT", "dd-MMM-yy") == Utilities.formatDate(curDat,"GMT", "dd-MMM-yy"))
{
return;
}
}
var snRows = sourceSheet.getLastRow();
var dnRows = dstSheet.getLastRow();

var srcStock = sourceSheet.getRange("A2:A" + snRows).getValues();
var srcLastSale = sourceSheet.getRange("C2:C" + snRows).getValues();

var dstStock = dstSheet.getRange("A2:A" + dnRows).getValues();
var dstLastSale = dstSheet.getRange("Z2:Z" + dnRows).getValues();

for(var j=0;j<dnRows-1;j++)
{
dstLastSale[j][0]="n/a";
}
var flag = "true";
var foundStock;
for(var i=0;i<snRows-1;i++) //snRows
{
var sStockVal = srcStock[i][0];

//var foundStock = ArrayLib.indexOf(dstStock,0, sStockVal);

flag="false";
for(var j=0;j<dnRows-1;j++)
{
if(dstStock[j][0].toString().toUpperCase() == srcStock[i][0].toString().toUpperCase())
{
flag = "true";
foundStock = j;
break;
}
}
if(flag=="true")
{
dstLastSale[foundStock][0] = srcLastSale[i][0];
}
else
{
var dnRows1 = dstSheet.getLastRow()+1;
dstSheet.getRange("A" + dnRows1).setValue(srcStock[i][0]);
dstSheet.getRange(dnRows1,lCol+1,1,1).setValue(srcLastSale[i][0]);
for(var k=2;k<=lCol;k++)
{
if(dstSheet.getRange(dnRows1, k).getValue()=="")
{
dstSheet.getRange(dnRows1, k).setValue("n/a");
}
}
}
}
dstSheet.getRange(1,lCol+1).setValue(curDat);
dstSheet.getRange(2, lCol+1, dstLastSale.length, 1).setValues(dstLastSale);
}

 
5) Set it to run each day at 12:00. The code will save the new last sale price for monday to friday with one days lag.

Now we can move on to Maple.

In Maple run the following code to load the data:

X := proc (Url) local theDLL, URLDownloadToFile, myDirectory, myFile, Destination, DL;

theDLL := "C:\\WINDOWS\\SYSTEM32\\urlmon.dll";

URLDownloadToFile := define_external('URLDownloadToFileA', pCaller::(integer[4]), szURL::string, szFileName::string, dwReserved::(integer[4]), lpfnCB::(integer[4]), 'RETURN'::(integer[4]), LIB = theDLL);

if FileTools[Exists]("C:\\mydir") = true then FileTools:-RemoveDirectory("C:\\mydir", recurse = true, forceremove = true) else end if;

FileTools:-MakeDirectory("C:\\mydir");
myDirectory := "C:\\mydir";
myFile := "data1.csv";
Destination := cat(myDirectory, "\\", myFile);

DL := proc () local M;

URLDownloadToFile(0, Url, Destination, 0, 0);
M := ImportMatrix("C:\\mydir\\data1.csv", delimiter = ",", datatype = string);
M := Matrix(M, datatype = anything)

end proc;

return DL()

end proc:

data := X("https://docs.google.com/spreadsheets/d/1L5-yUB0EWeBdJNMdELKBRmBQ1JJ0QymrtDLkVhHCVn8/pub?gid=649021574&single=true&output=csv");
L := LinearAlgebra:-Transpose(data);

If you use your own spreadsheet you need to change the url to match that spreadsheet.
Select File -> Publish to the web in google spreadsheet

We can now run the portfolio optimization in Maple:

with(Statistics):
with(ListTools):
with(LinearAlgebra):
with(Optimization):
with(plots):

Nr, Nc := ArrayTools:-Size(L):
symb := L[1 .. 1, 2 .. Nc]:
LL := L[2 .. Nr, 2 .. Nc]:
Nr, Nc := ArrayTools:-Size(LL):

# Removing stocks with missing observations
for i to Nc do if Occurrences("n/a", convert(Column(LL, i), list)) >= 1 then AA[i] := i else AA[i] := 0 end if
end do;

DD := RemoveInRange([seq(AA[i], i = 1 .. Nc)], 0 .. 1):
symbb := DeleteColumn(symb, DD):
LLL := map(parse, DeleteColumn(LL, DD)):
Nr, Nc := ArrayTools:-Size(LLL):

# Calculate Return
for j to Nc do
for i from 2 to Nr do

r[i, j] := (LLL[i, j]-LLL[i-1, j])/LLL[i-1, j]

end do
end do;

RR := Matrix([seq([seq(r[i, j], j = 1 .. Nc)], i = 2 .. Nr)], datatype = float[8]);
n, nstock := ArrayTools:-Size(RR):

# Portfolio Optimization
W := Vector(nstock, symbol = w):
y := Vector(n, fill = 2, datatype = float[8]):
s1 := Optimization[LSSolve]([y, RR])[2];
Nr, Nc := ArrayTools:-Size(s1):

j := 0:
for i to Nr do if s1[i] <> 0 then j := j+1; ss1[j] := symbb[1, i] = s1[i] end if end do;

Vector(j, proc (i) options operator, arrow; ss1[i] end proc);
LineChart(s1);

eq1 := a(t)*(diff(a(t), t))+a(t)*(diff(a(t), t))*b(t)*(diff(b(t), t))*c(t)*(diff(c(t), t));
eq2 := a(t)*(diff(a(t), t))+a(t)*(diff(a(t), t))*c(t)*(diff(c(t), t))+a(t)*(diff(a(t), t))*b(t)*(diff(b(t), t))*c(t)*(diff(c(t), t));

DEplot({eq1, eq2}, [b(t), c(t)], t = 0 .. 1, b = 0 .. 1, c = 0 .. 1, [[b(0) = 1, c(0) = 1]], arrows = large);

Error, (in DEtools/DEplot/CheckDE) only derivatives of dependent variables can be present

DEplot({eq1 = 3*t^2, eq2 = 2*t^3}, [b(t), c(t)], t = 0 .. 1, b = 0 .. 1, c = 0 .. 1, [[b(0) = 1, c(0) = 1]], arrows = large);
Error, (in DEtools/DEplot/CheckDE) only derivatives of dependent variables can be present

HI there,

I m getting an error message . Could someone help me

v1 := int(cos(tau)*g(tau), tau = t0 .. t);
int(cos(tau) g(tau), tau = t0 .. t)
v2 := int(-sin(tau)*g(tau), tau = t0 .. t);
int(-sin(tau) g(tau), tau = t0 .. t)
soln := C1*y1+C2*y2+v1*y1+v2*y2;
C1 y1 + C2 y2 + (int(cos(tau) g(tau), tau = t0 .. t)) y1

+ (int(-sin(tau) g(tau), tau = t0 .. t)) y2
soln := combine(soln);
(int(-y2 sin(tau) g(tau) + y1 cos(tau) g(tau), tau = t0 .. t))

+ C1 y1 + C2 y2
eval(soln, t = t0) = 0, eval(diff(soln, t), t = t0);
C1 y1 + C2 y2 = 0, -y2 sin(t0) g(t0) + y1 cos(t0) g(t0)
solve({%}, {C1, C2});
soln := eval(soln, %);
Error, invalid input: eval received (C1*y1+C2*y2 = 0, -y2*sin(t0)*g(t0)+y1*cos(t0)*g(t0)), which is not valid for its 2nd argument, eqns

Hi everyone,

Consider this Maple 18 doc:

FareySeq.mw

I am having the same problem as in:

http://www.mapleprimes.com/questions/203593-Maple-13-Fast-Maple-18-Crawling#comment213548

Namely, the produced animated .gif contains only one frame. Animation doesn't work, even though the button of the animation is active after I click on frame #1.

I have acer's solution at the end, but it doesn't work in this case either.

Breaks in both Maple 13 and 18.

TIA,

Yiannis

Hello!

With Maple 2015, license Student, is it possible calculate distance between point and plane in R³ ? Package geom3d is recommended for this or another? How can I calculate this?

Thank you so much!