below code is calculate basis of kernel and kernel

i guess basis of image is 

remove(has, Ga, [r,u,v,w]); if this correct, i eliminate this, i can get the image

however it include variable 'a'

is it correct? if not, how to calculate? 

my final goal is to make unexact sequence into exact sequence

after union ideal I and J, how to remove the intersection of I and J in order to do primedecomposition?

I have a Matrix 'A' in maple with symbolic entries and I need to use that Matrix in matlab M file. The symbols with in 'A' are already initialized in M file.

How is this possible?

For Matrices 'A' of small dimension, I can use the codegeneration package to convert 'A' to Matlab Matrix and then copy paste.

But for matrices of 1000x1000, lets say, does not seem as a option.

Hello. 
I have a workcheet, that is incompleete, when using

  > with(XMLTools);
  > ParseFile("FILE");

i get the wollowing error, but i cant seem to fix it by my selv in an xml editor. 
Here you have the file:

9.mw

Hope you can help

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?