I solved IVP of ODE with solutions x(t),y(t),z(t) numerically.x(0)=y(0)=0,z(0)=0.25

How can I get the coefficients of the initial values in x(t)

i.e  assume the solution

x(t)=e1(t)x(0)+e2(t)y(0)+e3(t)z(0)+e4(t)

how can we get :

e1(t),e2(t),e3(t) and e4(t)

 

I recieved this comment on an older post of mine, and I was going to post a reply, but I noticed my reply was pretty long, and I was also missing a post for last friday, so I figured I'd cheat and make the reply a blog post of its own.

Very few people would ever need this, I think, even while programming. But sometimes the details of the call stack are just what one wants.

So, for example,...

> h:=proc(x)
>   debugopts('callstack');
> end proc:

> m:=module() export f; local g;
>   f:=proc(x)
>     g(x^2);
>   end proc:
>   g:=proc(x)
>     global h;
>     h(x^3);
>   end proc:
> end module:

> m:-f(a);
[DEBUGSTACK, h, `debugopts('callstack')`, [a^6...

I've been working with fsolve and I'm having the following issue: fsolve({dif_Fliq=((nFsol-nFliq)/(psi_sol-psi_liq)),dif_Fliq=dif_Fsol},{psi_sol=0.1..0.9,psi_liq=0.1..0.9}); {psi_liq = 0.4447971078, psi_sol = 0.3459845910} how do I get the value of psi_liq and psi_sol and assign them to a different variable without having to copy the whole thing? i.e. if I use the following command assign(result_liq,psi_liq); I get result_liq=psi_liq = 0.4447971078 not what I want result_liq = 0.4447971078 I also tried: result_liq:=psi_liq and I ended up with the same result

An attempt at question 2 of the xkcd Velociraptor Math problem (mentioned on this blog post). The parameters and events facilities of

There is no released Classic interface to accompany the 64 bit version of Maple (12, 13) for the 64 bit Windows XP64 operating system. Personally, I prefer running the Standard over the Classic interface, although sometimes I miss common subexpression display for lengthy symbolic output.
 
The Maple Classic interface appears to talk to the Maple kernel only over a socket (or similar), and the...

In this post I'll take a closer look at the ways in which Maple code can be thread unsafe. If you have not already seen my post on Thread Safety, consider reading that post first. As a brief review, a procedure is thread safe if it works correctly when run in parallel.

The most obvious way in which procedures can be thread unsafe is if they share data without synchronizing access (using a Mutex, for example). So how can two threads share data?

I am going to wander away from parallel programming in Maple, to talk about GPU programming. However I will show an example of connecting Maple to a CUDA accelerated external library, so that's close enough.  This post is not intended to be a tutorial on CUDA or OpenCL programming, but an introduction to how the technology works.

Download 10597_asymp.mws
View file details

Hi
I am facing problem with maple. I try to get a asymptotic series of a transcendental equation but I get failed again and again. If You people could help me.

I am waiting for ur kind responce.

Thanks

 

I'd like to start by thanking all those readers who left feedback on my last post. It was good hear that most of you enjoy reading my posts and that they are generally helpful. I would like to encourage you to continue posting feedback, especially questions or comments about anything that I fail to explain sufficiently.

The following is a discussion of the limitations of parallel programming in Maple. These are the issues that we are aware of and are hoping to fix in future releases.

I have two sets A={a,b,c} and B={w,x,y,z}. I have to make sets like following

{(w,{a,c}), (y,{a,b,c}), (z,{b})}

 How I can know the total number of such sets from A and B?

 

I noticed that maple's command Transpose can mean two different things:
 

 

ListTools

The Transpose function transposes a list of lists
 

 

LinearAlgebra

The Transpose function computes the transpose of a Matrix, Vector, or scalar.
 

 

 

To highligt this I have selected two examples:

 

 

 

I don't really have anything prepared for today, so I'd like to ask you a few questions about the posts I've made so far.  My goal for this blog was to give typical Maple programmers the information they needed to start trying parallel programming. 

1. My posts progressed fairly quickly, building up to the the Task Programming Model. Did I move too quickly?  Were there topics that I did not explain well enough or that you felt needed more explaination?
2. As my goal was to present the Task Programming Model, I skipped a deeper explaination of Threads:-Create style of programming.    Would you like to know more about that type of low level threaded programming?
3. Most of the examples I used were artifical ones that illustrated the points I was trying to make.  Would you have prefered real world examples instead?
4. Did reading my posts get you to actually try writing a parallel algorithm?  If yes, did you succeed?  If no, why not?
5. Was the formatting ok, especially the code?  Each post included an worksheet containing the examples from the post so I did not worry too much about ease of copy and paste.
6. What else would you like to know about?  I am definitely planning a post on GPU computing, but since it is not really a Maple topic I delayed it till after I am finished with the Maple topics.

Any other feedback you would like to provide would also be appreciated, although I'd like to keep focused on the topics discussed in my blog, and less about Maple in general.

In my previous posts I have discussed various difficulties encountered when writing parallel algorithms. At the end of the last post I concluded that the best way to solve some of these problems is to introduce a higher level programming model. This blog post will discuss the Task Programming Model, the high level parallel programming model introduced in Maple 13.

I have been trying to calculate the Lyapunov exponent of the Rossler oscillator with an intention of finding it at higher precision. When i calculate the same at 16 digit precision on my 64 bit workstation or on my macbook using maple it takes much time (15 hrs have already been passed) and slows down the system badly. Can somebody plz help me out how to make my code more efficient and it runs in such a way that minimum resources of my computer are exploited.