Hi, the fsolve works perfectly fine when otaining SOL2 below. However, as soon as I change the values of paramters (from 0 values in SOL2 to t_prof = .16, t_r = .16, t_w = .16 in SOL3) fsolve simplmy reproduces the eval results and would not solve! Please your help is highly appreciated. Thanks in advance!

restart;
eq1 := ((1-s)*w)^(1-gamma)*r^gamma-P_v = 0;
eq2 := (1-s)*w*H/M-(1-gamma)*P_v*q/phi = 0;
eq3 := r*(K+K_f)/M-gamma*P_v*q/phi = 0;
eq4 := prof-M*p_d*q+(1-s)*w*H+r*(K+K_f) = 0;
eq5 := -tau*y_x+q-y_d-y_g = 0;
eq6 := p_d-sigma*P_v/((sigma-1)*phi) = 0;
eq7 := y_d-M^(lambda-1)*Y*(p_d/P)^(-sigma) = 0;
eq8 := y_x-F_f*(tau*p_d)^(-sigma_x) = 0;
eq9 := y_f-M_f^(lambda-1)*Y*((1+z)*tau*p_f/P)^(-sigma) = 0;
eq10 := (G*alpha_g+Y)^(-alpha_c)/P-A*H^alpha_h/((1-t_w)*w) = 0;
eq11 := P*Y-(1-t_w)*w*H-(1-t_r)*r*K-(1-t_prof)*prof+T = 0;
eq12 := P*Y-M*p_d*y_d-M_f*(1+z)*tau*p_f*y_f = 0;
eq13 := (1-t_r)*r-r_f = 0;
eq14 := G-y_g*M^((lambda_g-1+sigma_g)/(sigma_g-1)) = 0;
eq15 := s*w*H+M*p_d*y_g-T-t_w*w*H-t_r*r*(K+K_f)-t_prof*prof-z*M_f*tau*p_f*y_f = 0;

SOL2 := fsolve(eval({eq1, eq10, eq11, eq12, eq13, eq14, eq15, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, [A = 7.716049383, K = 3.2, K_f = 0, alpha_c = .8, alpha_g = .2, alpha_h = 4, lambda = .5, p_f = 1, phi = .625, s = 0, sigma = 5, sigma_g = 5, sigma_x = 5, t_prof = 0, t_r = 0, t_w = 0, tau = 1.1, y_d = 1, y_f = 1, y_g = .4, y_x = 1, z = 0, lambda_g = .5, gamma = .25, gamma = .25]), {F_f, G, H, M, M_f, P, P_v, T, Y, p_d, prof, q, r, r_f, w}, {F_f = 2 .. 4, G = .1 .. .4, H = .4 .. .8, M = .4 .. .6, M_f = .3 .. .7, P = .8 .. 2, P_v = .4 .. .8, T = -.2 .. .3, Y = .7 .. 3, p_d = .9 .. 1.6, prof = .1 .. .4, q = 1.5 .. 2.7, r = 0.4e-1 .. 0.9e-1, r_f = 0.4e-1 .. 0.9e-1, w = .8 .. 3});

SOL3 := fsolve(eval({eq1, eq10, eq11, eq12, eq13, eq14, eq15, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, [A = 7.716049383, K = 3.2, K_f = 0, alpha_c = .8, alpha_g = .2, alpha_h = 4, lambda = .5, p_f = 1, phi = .625, s = 0, sigma = 5, sigma_g = 5, sigma_x = 5, t_prof = .16, t_r = .16, t_w = .16, tau = 1.1, y_d = 1, y_f = 1, y_g = .4, y_x = 1, z = 0, lambda_g = .5, gamma = .25, gamma = .25]), {F_f, G, H, M, M_f, P, P_v, T, Y, p_d, prof, q, r, r_f, w}, {F_f = 2 .. 4, G = .1 .. .4, H = .4 .. .8, M = .4 .. .6, M_f = .3 .. .7, P = .8 .. 2, P_v = .4 .. .8, T = -.2 .. .3, Y = .7 .. 3, p_d = .9 .. 1.6, prof = .1 .. .4, q = 1.5 .. 2.7, r = 0.4e-1 .. 0.9e-1, r_f = 0.4e-1 .. 0.9e-1, w = .8 .. 3});

I have started to use Maple 18 and I have found that my scale starts at 2.*10^16 for all functions I graph no matter if its 3D or 2D. Is there a way to change this scale?

Dear, my attached file take alot time to exicute the result, Please solve this problem. I am waiting your kind response. 

Hirota_Bilinear_Method-3.mw

With my best regards and sincerely.

Muhammad Usman

PhD (Scholar) Mathematics

HITEC University, Taxila

Lecturer of Mathematics

Govt. Degree College Taxila Pakistan.

Visiting Lecturer of Mathematics

University of Wah, Pakistan.

Email: mughal1103@gmail.com

Mob #: +923152329276

1.

similar gap system, which function can output coset table ?

if find example below,

f := FreeGroup( "a", "b" );

g := f / [f.1^2, f.2^3, (f.1*f.2)^5];
tab := CosetTable(g, Subgroup(g, [g.1, g.2*g.1*g.2*g.1*g.2^-1]));

2.

if quotient group is for equivalence relations, what do relations [g.1, g.2*g.1*g.2*g.1*g.2^-1] for? which kind of relations is for input to subgroup parameter which is like relations?

3. equalivance relations in wiki mentioned is to partition cells, how do they satisfy symmetric, reflexive and transitive after partition? how to test these partition?

how to implement a free group with symmetric properties?

i find that equivalence relations = relexive, symmetric and transitive , is only these three ?

if implement symmetric properties, does it mean that adding equivalence relations to group?

how to implement free group with this?

i guess something like

f := FreeGroup( "a", "b" );
g := f / [equivalence relations....];

is there a library storing all kinds of equivalence relations?
https://en.wikipedia.org/wiki/Equivalence_relation
For quotient groups, is there IsEquivalenceRelation and IsEquivalenceClass like gap system?
www.gap-system.org/Manuals/doc/ref/chap33.html
like
f := FreeGroup( "a", "b" );
g := f / [f.1^2, f.2^3, (f.1*f.2)^5];

when i paste words here, most words missing, is there being hacked?

Hi,

  Suppose I have a probability distribution function f(x). I would like to generate points in a 

(1) given interval, can be up to infinity

(2) number of points

(3) satisfy the known probability distribution function f(x)

(4) possible other specifications

  Is there any build-in option in maple I could use? Or I have to write something manually?

Thank you very much!

I want to calculate

Lim _x-->1⁺ ([X²]-1)/([X]-1)

I let  x=1+c  and  _c-->0⁺  then we have [X²]=[c²+2c+1]=[2c]+1

Lim _x-->1⁺ ([X²]-1)/([X]-1) = Lim _c-->0⁺ [2c]/[c]

Maple shows the both limitations equal to zero

I want know how maple claculates it

Download Define_Question.mwDefine_Question.mw

If there is  an equation or are several equations, I need to obtain all the roots, how can I do???

fsolve ? rootfindings? or what?

If an examples of actual is given,  That will be perfect  !!!

Thanks 

with this simple circuit i found a wrong simulation result : look at the probe 3 value.  someone could help me please ?

The same error i found in this circuit (binary adder)

Hi, there

How can we plot the volume of revolution of the region between two curves y=2x-x^2 and y=sqrt(2x-x^2) around the x-axis?

Any help will be appreciated. I used a command in maple 13 but the result was a sphere!

Regards

Yegan

Hi, I am new to Maple. I want to plot a graph in a separate window from command-line maple. I tried the following code and getting this error

> plotsetup(window);
> plot(x,x=-10..10);
Plotting error, no plot device driver for plotdevice=window
>

I tried maplet but I am not getting the plot

> plotsetup(maplet);
> plot(x,x=-10..10);
Initializing Java runtime environment.
>

What should I do to get the plot??

Given an almost contact metric manifold M(\phi,\xi,\eta, g), we say
that M is a generalized Sasakian-space-form if there exist three functions f1, f2, f3
on M such that the curvature tensor R is given by

R(X,Y)Z=f_1{g(Y,Z)X-g(X,Z)Y}+f_2{g(X,\phiZ)\phiY-g(Y,\phiZ)\phiX+2g(X,\phiY)\phiZ}+f_3{\eta(X)\eta(Z)Y-\eta(Y)\eta(Z)X+gg(X,Z)\eta(Y)\xi-g(Y,Z)\eta(X)\xi}

In (2n+1) dimensional generalized Sasakian space form M²ⁿ⁺¹(f_1,f_2,f_3), we have the following relations.

S(X,Y)=(2nf_1+3f_2-f_3)g(X,Y)-(3f_2+(2n-1)f_3)\eta(X)\eta(Y)

S(X,\xi)=2n(f_1-f_3)\eta(X)

C\bar(\xi,X)Y=[f_1-f_3-(r/2n(2n-1))][g(X,Y)\xi-\eta(Y)X]

P(X,Y)Z=R(X,Y)Z-(1/(n-1))[S(Y,Z)X-S(X,Z)Y]

R(X,Y)\xi=(f_1-f_3){\eta(Y)X-\eta(X)Y}

R(\xi,X)Y=(f_1-f_3){g(X,Y)\xi-\eta(Y)X}

for any vector fields X, Y on M, where R, S, C\bar, and r denote the Riemannian curvature tensor, Ricci tensor, concircular curvature tensor and scalar curvature of M²ⁿ⁺¹(f1, f2, f3), respectively 

Using above equations I have to evaluate P(C\bar(\xi,X)Y,Z)U.

Manually It is tedious job. Can we find the value by maple? Is there any option to solve these type of problems?

If yes, I can solve many more, which helps a lot in my work.. Thanks in advance.

Hello,

I have a procedure wich is something like 

F:=proc(a,b,l) limit(f(a,b),b=l); end proc;

For my specific problem, the limit always exists and is a well defined function

After that, I would like to generate other functions, depending on F(a,b,l), e.g.:

G:=proc(a,b,N) add(F(a,b,l),l=1..N); end proc;

However, when I tried to do that, Maple inserts in G an unevaluated F, depending on the results of the limits. This makes G not work properly (unevaluated), since the limit is not computed before G. I tried to use some intermediate steps, like using unapply command, but it does not works.

If someone can help, I would appreciate.