I am trying to solve 4 nonlinear equations for four variables using fsolve  and the output that i am getting is basically the same equations repeated after some time.  I even tried reducing one of the equations using assumptions from my side but it results in same behaviour..  Quite new to maple, would like some advice as to this behaviour. Thanks

 Here's the file

fsolve_1.mw

PS- using do loop is part of the solving so i cannot remove that

Occasionally I have need to import 1D Maple expressions into Mathematica, because I believe mma has superior simplify functionality (sorry Maple). Mathematica has its own syntax which (of course) differs to maple.  eg sin(x) in maple is Sin[x] in mma ; ln(x)=Log[x] ;....... I'm aware of the package MmaTranslator, but I think it only converts mma to maple.

so i need some code pls.

(i have bolded the functions which need to be converted)

eg: i/p:

(1/365)*(365*B21*EP*ln(1+i)*(1+i)^B21*hr*kw*p-365*EP*hr*kw*p*(1+i)^(-1+B21)*i-365*EP*hr*kw*p*(1+i)^(-1+B21)+365*EP*hr*kw*p+SC*i)/(hr*kw*(EP*ln(1+i)*(1+i)^B21*p-FIT*i*p+cos(omega)))

o/p:

(1/365)*(365*B21*EP*Log[1+i]*(1+i)^B21*hr*kw*p-365*EP*hr*kw*p*(1+i)^(-1+B21)*i-365*EP*hr*kw*p*(1+i)^(-1+B21)+365*EP*hr*kw*p+SC*i)/(hr*kw*(EP*Log[1+i]*(1+i)^B21*p-FIT*i*p+Cos[omega]))

Hi everyone

I've come across a RAM-issue.

My program creates a lot of plots using millions of variables-values and it eats up my computers memory since I use the plots to make an animation.

The plots are simple 3d plots created using the surfdata function.

How do I save plots into a single file as the program creates them concurrently such that after the program has terminated I can use the plots to create an animation?

and how do I efficiently delete variables I dont not need anymore?

Dear Friends

            Hope everything going fine with you. I want the numerical solution of nonlinear system of ordinary differential equations using RK method. The system of ODEs and their required results are present in attached file. I am waiting your quick response.

Help.docx

With my best regards and sincerely.

Mob #: 0086-13001903838

I read research paper about elasticity, and I don't not how can get the results for displacment ant shaer stress,, could you help me please.

please help me

Kind Regards

hi.i encounter error in pdsole equations with unknown parameter(N)??

please help me for solve it....thanks alotmaple_prime.mw

Hi everybody,

Ever since i updated to maple 2015, the output of calculations display a disproportionate amount of zeros. For instance:

If i write:

2.3*10^(-4)

i get:

0.000230000000000

Is there anyway to make maple NOT display the last 10 zeros? I've tried Tools ->Options ->Precision -> Round screen display to X. But then it rounds all results to that X amount of digits, wich is also anoying. Like, if i set X to be 10, it'll write 3.000000000 instead of just 3.

Hope someone can help me :-)


Anton

In Maple 2015.1 we have

restart;

solve([sin(2*x)/cos(x+3*Pi/2)=1,  x>-4*Pi, x<-5*Pi/2], x, allsolutions, explicit);

solve([sin(2*x)/cos(x+3*Pi/2)=1, x>0, x<2*Pi], x, allsolutions, explicit);

In the first example, the error message is not clear (actually there exists a unique root  x=-11*Pi/3), in the second example, one root  (x=5*Pi/3) is lost.

I'm trying to solve a system of equations thats expressed as a summation (the original has the summation symbol as opposed to 'sum'):

d_actual := solve(W_actual = sum(W_guess(def-asp_rad_inverse[i], E, asp_rad[i]), i = 1 .. n), def);

When n<5, I get an answer after a few seconds, but when n is higher, the program sits and 'evaluates' forever...I've waited up to 30 min. 

Background:

I have a plate with a number of hemispheres on the surface (# of impacted hemispheres given a force = n). Each has a unique radius and they're listed from largest to smallest in 'asp_rad[]'. I have an equation for the deflection of a single hemisphere as a function of Force and material properties that I have rearranged with respect to Force (W_guess(deflection, E, r)). There's an opposing plate that stays parallel to the original plate while pushing down on the asperities with Force W_actual. The total deflection of the opposing plate is the sum of: 1) the difference between the tallest radius and the impacted radius in question (asp_rad_inverse[]), and 2) the deflection of the impacted radius in question. 

I'm attempting to solve for the total deflection of the opposing plate via solving for the 'def' in the summation above, but when I run it, the program is not able to compute a solution.

Thanks in advance.

I want to solve an ODE from Game Theory, the Cournot competition.

It says

p(q1+r2(q1))+p'(q1+r2(q1))*r2(q1)-c2'(r2(q1))=0

 where, I think,

' means diff(,q1),

c2(q2)=c*q2 for a fixed c in [0,1]

and

p(q)=max(0,1-q).

So c2,p and r2 are functions.r2 goes from [0,inf) to [0,inf).

I look for r2, which should be r2(q1)=(1-q1-c)/2 when correctly solved.

However, the command dsolve says Error in dsolve (divison by 0).

 What is wrong? How do I obtain the solution for r2 in Maple?

Hello,

I tried to solve below equation, but it gives me zeros result. Please help me to find their inverse laplace. 

 It will be clearer if was pasted on Maple:

restart

Ps := [P[0], P[1], P[2], P[3], P[4]]:

eqs := [P[0](s) =~ (P[1](s)*mu[1]+P[2](s)*mu[2]+P[3](s)*mu[3]+P[4](s)*mu[4])/(s+lambda[1]+lambda[2]+lambda[3]+lambda[4]), P[1](s) = lambda[1]*P[0](s)/(s+mu[1]), P[2](s) = lambda[2]*P[0](s)/(s+mu[2]), P[3](s) = lambda[3]*P[0](s)/(s+mu[3]), P[4](s) = lambda[4]*P[0](s)/(s+mu[4])];

Ls := solve(eqs, Ps(s))[];

P(t)=~inttrans[invlaplace]~(rhs~(Ls), s, t);

Thank you

Wonder if this can be accomplished in Maple.

so I have a list of 100 items labeled {1..100} of various value {$100, $160, $220, ......  , }

the task is to distribute these items among 3 people A,B,C so they get an approximately equal share.

Adding the values and dividing by 3 gives the dollar total to aim for. 

This post has C.Love procedure for evenly sized groups

 http://www.mapleprimes.com/questions/200480-Product-Grouping

but what i want is a method for different sized groups. ie 25 items for A, 35 for B and 40 for C (user defined).

additionally there is a fixed constraint: A has been bequeathed items 1,4,8; B items 2 and 20; C item 50.

restart:
S:= {3, 4, 5, 6, 8, 9, 28, 30, 35}:
SL:= [A,B,C,D,E,F,G,H,I]:
assign(Labels ~ (S) =~ SL); #Create remember table.
AllP:= [seq(P, P= Iterator:-SetPartitions(S, [[3,3]], compile= false))]:
lnp:= evalf(ln((`*`(S[]))^(1/3))):

Var:= proc(P::({list,set}(set)))
local r:= evalf(`+`(map(b-> abs(ln(`*`(b[]))-lnp), P)[]));
end proc:

Min:= proc(S::{list,set}, P::procedure)
local M:= infinity, X:= (), x, v;
     for x in S do
          v:= P(x);
          if v < M then  M:= v;  X:= x  end if
     end do;
     X
end proc:

ans:= Min(AllP, Var);
              [{3, 9, 35}, {4, 8, 28}, {5, 6, 30}]
subsindets(ans, posint, Labels);
               [{I, A, F}, {B, E, G}, {C, D, H}]

I am attempting to write a Gaussian elimination routine to solve a system Ax = b using loops, but I have been having trouble.  Any help would be mcuh appreciated.  Thanks!

Need help for manipulating tensor with the physics package.

I ask some questions about this.  But each time, I am refer to the help pages.  If I ask again some help, it is because I can't not start with the information on the help file.  It is written for people that already know General Relativity (GR).

So this time, I have created a document (attach to this post) where I ask specific queations on manipulations.  My goal is to ccrreate a document that I will put on the Applications Center.  I promess that those who will help me on this will be cited in the document.  This way, I hope to create an introduction on how to use tensors for beginner like me.

Then, with this help, I am sure I will be able to better understand the help page of the packages.  I am doing this as someone who is starting to learn GR and have to be able to better understand the manipulations of tensor and getting the grasp of teh meaning of all those tensor.  For exemple, the concept of parallel transport on a curve surface.

Thank you in advance for all the troubling I give you with this demand.

Regards,Parallel_Transport.mw

--------------------------------------
Mario Lemelin

Maple 2015 Ubuntu 14.04 - 64 bits
Maple 2015 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

Are there any examples of network solutions such as mine systems in Maple.
I wrote in the Maple system of partial differential equations describing the process of filtration combustion.
I'm a novice. I do not quite understand how to solve it.
Online a lot of "simple" examples. I need something very similar to my case.

Model03.mw