Hi:

i will seprate coefficients of a,a^2,a^3 and this coefficients equivalent zero(0) in EQ below,how can i do this?

EQ:= 3.522088319*10^(-10)*a^2*q21(T)+3.522088319*10^(-10)*a^3*q22(T)-6.094634327*10^(-16)*a*q41(T)-6.094634327*10^(-16)*a^2*q42(T)+6.094634327*10^(-16)*a*q51(T)+6.094634327*10^(-16)*a^2*q52(T)-1.106493792*10^(-9)*a*q11(T)-1.106493792*10^(-9)*a^2*q12(T)+2.429209785*10^(-24)*a^2*(diff(q21(T), T, T))+2.429209785*10^(-24)*a^3*(diff(q22(T), T, T)) = 0

I have a table of data arranged in the following columns: Year, Jan, Feb, Mar,..,Dec in Excel or csv

I would like to import this data into a time series in Maple so I can plot the result.

What is the simplest way of doing this?

Hi,

I have a mistake in these lines but I d'ont know the exact error. There is no numeric value of g(1,1)=4???????????????  Thanks for your help.

restart;
vvv:=n^2+3*m;  
g:=(n,m)->vvv;
g(1,1);  No display of the value 4
                               

Hi,

I have these lines in my code, the function alpha:=(m,n,1)->...... is defined in the last line of the code below, but when I want to calculate alpah(1,1,1), there is no numeric value. Why????  Many thinks

 restart:
with(LinearAlgebra):
with(plots):
with(PDEtools):
with(IntegrationTools):
interface(rtablesize=20):
d:=1: N:=2: a:=1: h:=a/N:
Kernel(x,y):=ln(abs(x-y)):

 eq2:=int(Kernel(x,y)*phi(y),y=-a..a)=sum(int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h),n=-N..N-d):

Approximate the integral
eq3:=phi->int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h)=add(beta[i,n]*phi((n+i-1)*h),i=1..d+1):
eq4:=int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h)=add(alpha[m,n,i]*phi((n+i-1)*h),i=1..d+1):

Compute the coefficients  beta[i]

Vct_basis:=[seq(x^i,i=0..d+1)]:
fct:=[seq(unapply(Vct_basis[i],x),i=1..d+2)]:


sys:=[seq(eq3(fct[i]),i=1..d+1)]:
x:=m*h:
w := [seq(beta[i,n],i=1..d+1)]:
M,b := GenerateMatrix(sys,w):

M1:=-M: V:=-b:

Vect_beta:=(M1)^(-1).V:# Vect_beta is a vector.

beta[1]:=Vect_beta[1]:  
alpha1:=(n,m)->beta[1];
(n, m) -> beta[1]
alpaha(1,1);
                                alpaha(1, 1)  ??????????????????? No numeric result

Hi all,

Command collect can view a expression a as a general polynomial in specific indeterminate x with the calling sequence:

collect(a, x, form, func)

where form and func are optional arguments.

Then, I encountered a form like this:

collect(a, x, normal)

I'm not sure what does normal means in this case and cannot find the related information from Maple Help.

Does anyone know this? I'd appreciate any help on this topic.

Thanks a lot :)

(Sorry to not offer the worksheet file due to a private issue)

I'm trying to interpolate values using spline (data,x,degrees=2).

So I'm setting curvefit:=spline(data,x,degrees=2);

fsolve(curvefit=4,x);

However, I'm trying to solve this over a range of values in matrix M:=1..2000;

I'm sifting through the online documentation, but I can't find a proper way to do this. Ideally, i'd want to store the results as a seperate matrix.

It's been a while since I touched Maple. Just picked it up again to do some stat work, and I'm feeling a bit out of my element having been to accustomed to things like Matl.... (I mean the software that shall not be named).

Is there a way to convert a single array/matrix into multiple arrays?

 Scenario:

It's easy just to copy and past this as a single element within a worksheet. However, I want to get it so that it's formated like data = [[1,1],[2,7],[3,8],etc...]

Is there a way to reformat into individual elements once everything is in one large array? I'm basically copy and pasting a huge mess of cells from Excel and am looking to see if there's a simplier solution than manual data entry.

Hi,

I have an error how can I fix the error. Thank you.

Matrix_Fred:=proc(N,d,lambda)
local Kernel,K,Fredholm,phi,MatA,eq2,eq3,eq4,Vct_basis,fct,sys,M,b,M1,Vect_beta,W,i,beta,alpha,eq5,eq6,x,Subs1,Fredholm_stencil,Stencil,w,V,sys1,sys2,Subs2,sys3,Sol_phi;
Fredholm:=phi(x)=f(x)+lambda*int(Kernel(x,y)*phi(y),y=-a..a);
Kernel:=(x,y)->ln(abs(x-y));
eq2:=int(Kernel(x,y)*phi(y),y=-a..a)=sum(int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h),n=-N..N-d);
eq3:=phi->int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h)=add(beta[i,n]*phi((n+i-1)*h),i=1..d+1);
eq4:=int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h)=add(alpha[m,n,i]*phi((n+i-1)*h),i=1..d+1);
Vct_basis:=[seq(x^i,i=0..d+1)];
fct:=[seq(unapply('Vct_basis[i]','x'),i=1..d+2)];
sys:=[seq(eq3(fct[i]),i=1..d+1)];
w := [seq(beta[i,n],i=1..d+1)];
M,b := GenerateMatrix(sys,w);
M1:=-M; V:=-b;
Vect_beta:=(M1)^(-1).V;
for i from 1 to d+1 do
beta[i,n]:=Vect_beta[i];  
end do;
x:=m*h;
for i from 1 to d+1 do
alpha[m,n,i]:=unapply(beta[i,n],n,m);
end do;
eq5:=lhs(eq2)=sum(rhs(eq4), n=-N..N-d);
eq6:=subs(x=m*h,subs(lhs(eq5)=rhs(eq5) ,Fredholm)):
Subs1 :=[seq(phi(m*h)=phi[m], m=-N..N)];
Fredholm_stencil:=subs(Subs1,eq6):
Stencil[1] := unapply(Fredholm_stencil,m,lambda,phi,f);
sys1 := [seq(Stencil[1](m,lambda,phi,f),m=-N..N)]:
sys2:=subs(Subs1,sys1):
Subs2:=[seq(f(m*h)=f[m], m=-N..N)]:
sys3:=subs(Subs2,sys2):
Sol_phi := [seq(phi[i],i=-N..N)]; # The unknown vector must be computed.
MatA := GenerateMatrix(convert(sys3,list),Sol_phi)[1];
return MatA;
end proc:




d:=1; N:=2; lambda:=3/Pi;
Matrix_Fred(N,d,lambda);
                                      1
                                      2
                                     3
                                     --
                                     Pi
Error, (in rtable/Power) singular matrix

Hi all I am truing to write a maple program for Newton's Method for finding a zero of a function.  I got pretty far but I am trying to get it to list all the values not just one.  Any help would be great:

Here is  what I got:

0.2000000000

0.2016393443

0.2016396757

What I want is if I type in 2 (like the last line I showed) I want the output to show for n=0, n=1, =2.  Instead it is only showing the output just for 2.  Any help would be great thanks, Matt.  

Hi,

I have a code compute some function : 

         alpha1:(n, m,1) -> (n + 1) (int(K(|m h - y|), y = n h .. (n + 1) h))

              int(K(|m h - y|) y, y = n h .. (n + 1) h)
            - -----------------------------------------
                                  h                    
         alpha2:  (n, m,2) -> -n (int(K(|m h - y|), y = n h .. (n + 1) h))

                int(K(|m h - y|) y, y = n h .. (n + 1) h)
              + -----------------------------------------
                                    h                    
and  I have a Matrix "MatA" .

My aim, when I give the value of the Kernel K used in alpha1, and alpha2, like K(x,y)=ln|x-y| , I want a numerci Matrix.

How can I do it.

Many thinks.

Fred.mw

Hi,

m := Matrix(8, 2, [1, 2, 1, 4, 4,1,4,7,1, 2,4,1,1,4,4,7]);
plot(m)

I would like to add a mame of each point in the graph as: alpha:=(1,2), beta:=(1,4) , gamma:=(4,1),eta:=(4,7)

and the symbol of each point is a bleu sphere.

and how to remove axis borders, tick marks.

Many thinks

Several years ago, I used to plot, in Maple 7 I think, 3D scalar functions by using procedures to create a 3D mesh and populate the data points before I could use plot3d.

I wonder if there is a more convenient (least coding) way to do it today? Consider for example f(x,y,z) = x^2 + y^2+z^2

A way to visualize a number of concentric isosurfaces of f is to loop with (is there a tag to write this in code block?): 

for i to 10 do iso[i] := implicitplot3d( f = i, x = -10 .. 10, y = -10 .. 10, z = -10 .. 10) end do

display(seq(iso[j], j = 1 .. 10))

Of course, visualizing the output is another issue.  I wish for an app to explore 3D data like Paraview. It does not have to be as sophisticated, but to display standard elements like isosurfaces and arbitrary cutting plane views would suffice.

If you know a package/app/procedure in Maple of this nature, please share it here. Thank you

I am trying to solve an equation by continuation method (one vairable and one parameter) containing more than 1000 terms and also hight exponents like 800. I want to get omega' values for different 'p' values.The maple file  is appended below. the link is maple_equation.mw

Many thanks in advance.

How do I find the constants in a solution containing whittaker functions?

the boundary conditions are:

c(x,0)=1

c(0,y)=1

[c'(1,y),x]=0