Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

i have already release memory with gc() in the most inner for loop ]

,but it still use up my memory

once it printed two lines each time, actually it can release memory

but why it still use up memory

read "D:\\hello.m":
with(combinat):
appendto("D:\\gogogo.txt");
for ii from 1 to 2040 do
for jj from 1 to 2040 do
for kk from 1 to 2040 do
if ii < jj and jj < kk then
if mm[ii] <> mm[jj] and mm[jj] <> mm[kk] then
sol := dsolve([mm[ii],mm[jj],mm[kk]],[a(t),b(t),c(t)]);
print([mm[ii],mm[jj],mm[kk]]);
print(sol);
print("***");
end if:
end if:
od:
gc();
od:
od:

 

i can only think

use c# to generate code with file write, nested for loop every 1000 records 

is it the only way to do?

 

int counter = 0;
for (int i=1; i<=2040; ++i)
{
for (int j = 1;j <= 2040; ++j)
{
for (int k = 1; k <= 2040; ++k)
{
if(counter % 1000 == 0)
{
Console.Write(i.ToString() + ","+j.ToString() + ","+k.ToString());
Console.Write("\r\n");
}
counter = counter + 1;
}
}
}

 

for ii from 1 to 2039 do
for jj from 1 to 2039 do
for kk from 1 to 2039 do
print(hello(ii,jj,kk));
od:
od:
od:

Dear Community,

Would someone have a good and easy to understand/implement description of the Den Iseger algorithm for the numerical inversion of Laplace transform? Even better if someone would have a Maple script to do it, that would be superb.

Tx in advance,

best regards

Andras

Hello guys

 

I am trying to solve a non-linear 2nd order diffential equation in order to describe the turbulence phenomena for tip vortices .

The initial conditions that i have are

b(0)=1

b(8)=exp(.1842068074*ln(-(2.*(-1.+b(x)))/(2.*b(8)-7.))+9.670857391)

Moreover i know that b from 0 to 1 equals one and then it  increases propably parabololoid(from the experimntal measurements)with respect o Re(Reynolds)

I tried to solve the eqution using maple solve de interactivelly but doesnt seem to work

Do you have any suggestions?

Thank you

giorgio_3.mw

 

I updated the OSX from El Capitan 10.11 to Sierra 10.12.

After then, whenever I try to input Japanese characters, Internal Error notification appears, and it crashes.
As it is not realistic to go back to the old OS, as the time when I did back up was a little before, I wouldn't do so.
Then, I might end up torelating not using Japanese until the next version of maple appears next year.
Are the situations like this ?

I know there might not be any response as this is about inputting Japanese characters.

Best wishes.

taro

Hello

I am trying to solve interactively the following DE with initial condiotion b(0)=1 and maple freezes as well as my pc.Can you help me?

Hello, 

I have a PDE system. When I use pdsolve it gets me the messege " pdsolve->Warning: System is inconsistent". Is there a way I can see which equations breaks the system down? 
For this system, it's difficult to see from ayeball where the problem is. 
Thank you! 

test.mw

When I use the Determinant function on a matrix with (single variable) polynomial entries with real coefficients I often get an incorrect answer. I know the answers are incorrect because they have a higher degree or a lower lowest degree than is possible given the matrix elements.

However, when I replace the coefficients in the polynomials with rational numbers or I put in the option method=minor, I get the correct answer.

The problem seems to be roundoff error. However, the important error is not simply small changes in the resulting polynomial. The important error is the presence of entirely incorrect powers of the variable and not with very small coefficients.

How does this happen and why does the help page for Determinant( ) not warn of this behavior? In particuiar, why does the help page not say that using Gaussian elimination (i.e., the default) will often give incorrect answers in such cases, but using method=minor will work? Is this behavior known? I cannot find any reference to it on the internet.

 

Hello all. Is there any solution for the indefinite integralBadIntegral.mw

``

``

``

``

int((B*alpha[1]+(1/8)*`&Nu;u`)*HeunT((1/16)*alpha[1]*(8*B*alpha[1]+`&Nu;u`)*3^(2/3)*(2*alpha[2]+1)^2/(alpha[2]*(alpha[1]*alpha[2]*(8*B*alpha[1]+`&Nu;u`))^(1/3)), 0, (1/2)*alpha[1]*(8*B*alpha[1]+`&Nu;u`)*3^(1/3)/(alpha[1]*alpha[2]*(8*B*alpha[1]+`&Nu;u`))^(2/3), (1/3)*3^(2/3)*(alpha[1]*alpha[2]*(8*B*alpha[1]+`&Nu;u`))^(1/6)*y), y)

int((B*alpha[1]+(1/8)*`&Nu;u`)*HeunT((1/16)*alpha[1]*(8*B*alpha[1]+`&Nu;u`)*3^(2/3)*(2*alpha[2]+1)^2/(alpha[2]*(alpha[1]*alpha[2]*(8*B*alpha[1]+`&Nu;u`))^(1/3)), 0, (1/2)*alpha[1]*(8*B*alpha[1]+`&Nu;u`)*3^(1/3)/(alpha[1]*alpha[2]*(8*B*alpha[1]+`&Nu;u`))^(2/3), (1/3)*3^(2/3)*(alpha[1]*alpha[2]*(8*B*alpha[1]+`&Nu;u`))^(1/6)*y), y)

(1)

``

``

``

``

``

``

``

 

Download BadIntegral.mw

?

 

Thanks

i copy maple code from notepad to maple in maple window,

there is no error

my function in the code

explicit define parameters are Local type

for example

appendto("...");

func1(aaa)

Local aaa;

 

but when i run cmaple to read the code text file in window 8

it return error

missing operator, syntax error

at Local aaa;

originally 

i have defined

Local aaa, ii;

for ii from 1 to nops(aaa) do

etc.

but it has error too,

then i change to one by one

Local aaa;

Local ii;

still have error at Local aaa;

if i m working in maple 13 i have to solve a non linear integral equation then what will be the steps to use the do loop.

 

 

Hello all,

 

Is here anyone who could help me with this Integral.

I want to determine P(y=2) in which the Nu is a constant and also is unknown.

In deed I want to derive the P which only has one unknown parameter(Nu). But I cannot solve this indefinite integral although the \alpha_1 , \alpha_2 and B are known constants. The Nu parameter mustn't be valued preferely.Indefinite_INtegral.mw

 


Here is my code. I am deriving soil compression line for oedometric testing with porosity, below is a part of the derivation.

I would like to differentiate the (Sst/Vst)/(Spt/Vpt) by dsm first then integrate it with dsm ranges from 0 to dmax to get my final answer...

Hope someone can help me on writing the code... Thank you very much!!!

restart

``

Sst := `&beta;ss`*[dsx(dsm)^(2-Ds)-dsm^(2-Ds)]/(2-Ds);

`&beta;ss`*[dsx(dsm)^(2-Ds)-dsm^(2-Ds)]/(2-Ds)

(1)

Spt := `&beta;ps`*[dpx(dsm)^(2-Dp)-dpm(dsm)^(2-Dp)]/(2-Dp);

`&beta;ps`*[dpx(dsm)^(2-Dp)-dpm(dsm)^(2-Dp)]/(2-Dp)

(2)

Vst := `&beta;sv`*[dsx(dsm)^(3-Ds)-dsm^(3-Ds)]/(3-Ds);

`&beta;sv`*[dsx(dsm)^(3-Ds)-dsm^(3-Ds)]/(3-Ds)

(3)

Vpt := `&beta;pv`*[dpx(dsm)^(3-Dp)-dpm(dsm)^(3-Dp)]/(3-Dp);

`&beta;pv`*[dpx(dsm)^(3-Dp)-dpm(dsm)^(3-Dp)]/(3-Dp)

(4)

``

(dpm(dsm)/dpx(dsm))^(2-Dp) = (dsm/dsx(dsm))^(2-Ds);

(dpm(dsm)/dpx(dsm))^(2-Dp) = (dsm/dsx(dsm))^(2-Ds)

(5)

dpx := proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds) end proc;

proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds) end proc

(6)

dpm := proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds) end proc;

proc (dsm) options operator, arrow; C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds) end proc

(7)

eval(Sst*Vpt/(Vst*Spt));

`&beta;ss`*[dsx(dsm)^(2-Ds)-dsm^(2-Ds)]*(3-Ds)*(2-Dp)*`&beta;pv`*[(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds))^(3-Dp)-(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds))^(3-Dp)]/((2-Ds)*`&beta;sv`*[dsx(dsm)^(3-Ds)-dsm^(3-Ds)]*`&beta;ps`*[(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds))^(2-Dp)-(C1*dsx(dsm)*(dsm/dsx(dsm))^((3-Ds)*Ds+Dp-Ds))^(2-Dp)]*(3-Dp))

(8)

``

``

``

``


Download GEO_Assignment_3.mw

First 156 157 158 159 160 161 162 Last Page 158 of 2218