Hello,

I have two sets of data representing two function that depends on x with a parameter A.

I need to do a fit on both data series at the same time so to fit with the best parameter A.

Here is how I do a fit on one function

> f(x):=A*cos(x-B)^(2);               

> g(x):=A*(cos(x-C)^(2)+ sin(x-C)^(2))^2; 

> fit1 := Fit(f(x), r, x, parameternames = [A, B, C...

hey friends..
Am trying to solve 2d poission equation with finite difference method and using different scheme is got

u_i-1,j -8u_i,j + u_i+1,j + 3u_i,j-1 +3u_i,j+1 = -h *16

and i have made code for this which is

> restart;
> with(LinearAlgebra);
> printf(`The lefthand endpoint on the X-axis is 0.\n`);
> OK := FALSE;
> while OK = FALSE do printf(`Input...

r-x^4
u-(x^3)*y
v-x*y^3
w-y^4

Step 1
with(Groebner):
K := {r-x^4,u-(x^3)*y,v-x*y^3,w-y^4};
G := Basis(K, 'tord', degrevlex(r,u,v,w));

Step 2
G intersect k[r,u,v,w]

G2 := Basis({w*G[1],w*G[2],w*G[3],w*G[4],w*G[5],w*G[6],w*G[7],w*G[8],w*G[9],w*G[10],w*G[11],w*G[12],w*G[13],w*G[14], (1-w)*K[1], (1-w)*K[2], (1-w)*K[3], (1-w)*K[4]}, 'tord', degrevlex(r,u,v,w));

G2 has too much though include the correct solution

correct solution...

3*rho1 - 2*rho2 + rho3 - rho4 = -1

4*rho1 +   rho2 - rho3        = 5

original without cost function:

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,(y1^1000)*(y2^1)*(y3^1)*(y4^100)- x1*x2*w + 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));

after have cost function 1000*rho1 + rho2...

with(Groebner):
K := {y1-(x1^3)*(x2^4),y2-(x2^(1+2))*(w^2),y3-(x1^(1+1))*(w^1),y4-(x2^1)*w,x1*x2*w - 1};
G := Basis(K, plex(x1, x2, w, y1, y2, y3, y4));
Reduce((x2^(5+1))*(w^1), G, plex(x1, x2, w, y1, y2, y3, y4));
(1,1,0,2)

subs({rho1=1,rho2=1,rho3=0,rho4=2},3*rho1 - 2*rho2 + rho3 - rho4);
subs({rho1=1,rho2=1,rho3=0,rho4=2},4*rho1 +   rho2 - rho3       );

can also reduce with G8,G7

R1 := Reduce((x2^(5+1))*(w^1), {G[7...

Here is my sheet. non_linear_P2_last_q.mw

In this problem I have achieved every thing what I aimed, thnaks to the wonderful people here!!!

I am looking for a way to do it in a more structured way meaning the calculation of Jacobian in the attached sheet is somewhat manuel, though I wonder does maple has some other way?

Hi there i have function as:

f(B,F)=

where: F,B>0 and 

,

such that: 

h=B^2;

I need to plot the surface (contour plot) f(B,F) at 0 and -ve points i.e f(B,F)<0.

Thanks

Hi all I have the following ODE

ODE_1:=diff(w(r),r)+R*(diff(w(r), r))^3=K/r;

ODE_T_1:=collect(algsubs(w(r)=u(r)+(r-1)*(U-0)/(delta-1), ODE_1),diff(u(r), r)) ;

 and when I try to do the following integration,

eq1:=int(phi[i](r)*ODE_T_1,r=1..delta) assuming delta > 1;

it gives the following error

"

Error, (in assuming) when calling 'int'. Received: 'wrong number (or type) of arguments:

i am doing a Q Laplace and finally summation it

however, i find an AccurateQSummation which is not like Summation

Should i use it, if not, use SumTools' summation how to simplify it?

restart;
with(QDifferenceEquations):
qexp := QPochhammer(-(-p*t), q, infinity);
f:= qexp*t;
b := subs(t=q^n,f);
QSimplify((1-q)*AccurateQSummation(b*q^n,n=-infinity..infinity));

with(SumTools):
QSimplify((1-q)*Summation(b*q^n,n=-infinity..infinity));

Is it possible for someone to provide the roller coaster maplesim model?

Hi

i know to make an asymptotic expansion one uses asympt(1+x,x); how do you make it go to -infity rather than plus infinity?

Hi,

I'm trying to modify the size of the window where an animation is shown. Tipically, this window has axis going from xmin to xmax and ymin to ymax but this causes some scaling problems when xmax-xmin>>ymax-ymin. I would like to focus in a certain region of my animation and change this viewpoint depending on the position of the elements of the animation. Can anybody help me?

This is my code:

> tf := 5;
> nfr := 201;

I have a huge PDE system included reaction-diffusion system with 4 independent and 12 independent variables. I read on the maple help for pdsolve/numeric/method. I found it is only for single partial differential equations in two independent variables. Is there anyone knows the ability of Maple for such a system?

1.

qRacah := qhyperterm([q^(-x), alpha*q^(-x)/(gamma*delta)], [alpha*q], q, gamma*delta*q^(x+1)*t, k)*qhyperterm([beta*delta*q^(x+1), gamma*q^(x+1)], [beta*q], q, q^(-x)*t, k);

i use inverse laplace to transform ordinary generating function into exponential generating function

can not further evaluated

2.

i do not know whether x is main variable or q is main variable or both are.

ff := eval(subs(x=0, exponential_genfun));

Hi,

I am aplying a simple energie balance with some simplification, and i am implementing it on maplesim.

The component has the following equations:

Am i missing something? Do i need any aditional equation to maplesim know to where should the fluid be moving when Qjacto is positive or negative? Why do i have...