I have calculated coefficients in maple using a "for...end do" command. For example,

for n from 0 to 1 do A[n] = int(((2*n+1)*(1/2))*simplify(LegendreP(n, cos(theta)))*sin(theta), theta = 0 .. (1/2)*Pi) end do;


A[0] = 1/2

A[1] = 3/4

Here, I can calculate as many coefficients as I want. But, how do I use these coefficients in the following line? For example, I need to calculate a summation. But for each term in the summation, I want to input the coefficients above to the corresponding term in the sum. So. the following summation

 A[0]+A[1]*r*cos(theta)+...

will automaticlly attach the numerical values of the coefficients calculated above to the corresponding term in the sum.

I hope this makes sense.

I know the Wronskian command. I want to use this command for detecting linearly dependence or independence of some polynomials. I know that the polynomials  are  linearly independent if the Wronskian is not zero.  Conversely, if the Wronskian vanishes  then the polynomials are linearly dependent. Now I want to know that how can I find the coefficient vector if the polynomials are dependent by Wronskian command?

For example if [f1,f2,f3] be a list of polynomials s.t. a1f1+a2f2+a3f3=0. How can I find a1, a2, a3? 

I would like to overwrite a plot that's updated within a loop over time. 

I used animate already, but the problem is that it's  only updated at specific points in time (not regularly distributed) and with animate I can only use frames that are equally divided over the wholre range. To update each time step costs unneeded memory and time. 

I hope that there exists a better solution.

thanks for your help.

hi,

how can I convert the decimal representation of the symbolic expressions: I dont want the result in terms of symbols like e^5 or ln(3) 

thanks

Hi,

I have a system of nonlinear equations in conjuction with boundary conditions, as below:

(all coefficients are constant.)

dsys3 := {a2*(diff(f1(x), x, x))+a3*(diff(f2(x), x, x, x))+a4*(diff(f2(x), x))+a6*(diff(f3(x), x))+a7*f1(x)+a1N*(diff(f3(x), x, x))*(diff(f3(x), x))+a2N*(diff(f3(x), x))*f3(x)+a3N*(diff(f3(x), x))*f3(x) = 0, a8*(diff(f2(x), x, x, x, x))+a9*(diff(f2(x), x, x))+a10*f2(x)+a11*(diff(f1(x), x, x, x))+a12*(diff(f1(x), x))+a13*(diff(f3(x), x, x))+a14*f3(x)+a4N*(f3(x)*f3(x))+a5N*((diff(f3(x), x))*(diff(f3(x), x)))+a6N*(diff(f3(x), x, x))*f3(x) = 0, (diff(f3(x), x, x, x, x))*a16+(diff(f3(x), x, x))*a18+(diff(f3(x), x, x))*a19+(diff(f1(x), x))*a22+(diff(f1(x), x))*a23+(diff(f2(x), x, x))*a24+f2(x)*a25+f2(x)*a26+f3(x)*a27+f3(x)*a29 = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}:
dsol5 := dsolve(dsys3, 'maxmesh' = 500, abserr = .1, numeric, range = 0 .. 1, output = listprocedure)

but i encounter an error as:

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

what is your suggestion to solve this nonlinear equations?

Thank you.

Hi everyone,

I have a great problem with the evaluation of following definite integral

> restart;

> int((t-x)^2)/(1+2t+(1/2)t^2-ln(t^2+2t+2)t-ln(t^2+2t+2)+arctan(1+t)t^2+2arctan(1+t)t+ln(2)t+ln(2)-(pi/4)t^2-(pi/2)t)^2,t=0..x)

I have tried different classical commands but Maple doesn't give an answer. Probably, it's just a silly fault.

Does anyone knows how to solve it?

Thanks.

Is it possible to make scatter-histogram plots in Maple13 like these

http://stackoverflow.com/questions/6508769/matplotlib-scatter-hist-with-stepfilled-histtype-in-histogram ?

I achieved the scatter-plot, but what about the additional histograms?

Hi,

I'm using Java OpenMaple interface to call Maple functions from Java on Ubuntu Linux. I'm following the examples provided on Maple website. My Java code:

import com.maplesoft.openmaple.*;

import com.maplesoft.externalcall.MapleException;

class test {

    public static void main( String args[] ) throws Exception {

        String[] a = {"java"};

        Engine  t = new Engine( a, new EngineCallBacksDefault(), null, null );

        try {

            t.evaluate( "factor(-8*((16*I)*g**4*(-4*mc**2+u)**(-2)*(4*mb**2-u)**(-2)*u**(-1)*(12*mb**2-u)*(-4*s*mb**2+16*mb**4-4*u*mb**2+t*u-4*t*mb**2)*s*(4*s*mb**2+4*mc**2*s+16*mc**2*mb**2-s**2-4*t*mc**2+t**2-4*t*mb**2)*(12*mc**2-u)+(16*I)*(4*s*mb**2+4*mc**2*s+u**2-4*mc**2*u+16*mc**2*mb**2-4*u*mb**2-s**2)*(12*mc**2-t)*g**4*(t-4*mb**2)**(-2)*(4*mc**2-t)**(-2)*(-t+12*mb**2)*(-4*s*mb**2+16*mb**4-4*u*mb**2+t*u-4*t*mb**2)*s*t**(-1))**2);" );

        }

        catch ( MapleException e ) {

            System.out.println(e);

        }

    }

}

Maple throws exception:

com.maplesoft.externalcall.MapleException: Error, (in evala/Factors) input must be polynomials over the integers

The input polynoial is valid (when coping and pasting `factor(...)` directly into Maple interface, all works). When trying to factor simple polynomials like "x**2 - I*x" all goes fine too. Tried both Maple 17 and 18 (Linux).

PS. There is no problem when running this code on OS X (I tried Maple 16, 17, 18),  but I need to be able to run Maple with Java on my Linux cluster.

Hello everyone,

How can unknown values in boundary conditions of system ODEs be determine by using shooting technique. See the work sheet here Sht.mw

Thanks.

How  we can decide whether a polynomial f is a linear combination of some polynomials g1,...,gm?

For example if f=x^2+y^2 and g1=y+x^2 , g2=y^2-y then f=g1+g2.

I have a set of equations ics_i and a set of variables v_i . How do i assign each solution to the variable before solving the next equation?

I tried- for i from 1 to 3 do solved_i := solve(ics_i, v_i);assign(solved_i); end. This didn't work for me.

Which Maple versions will run on Windows 10? Windows 10 is offered free to users of Windows 7 and 8.1. Several (including old) Maple versions run on Windows 7.

Hi,

I tried document blocks, copied my codes in an area then created a block however in collapse mode all of my codes are seen in a line in series! How can i hide my codes and only show output?

Hi everyone. Can enuone help with small parameter method?

restart;

eq := diff(x(t), `$`(t, 2))-epsilon(alpha*x(t)^4-beta)*(diff(x(t), t))-x(t)+x(t)^3;

subs(x(t) = 1+epsilon*x[1](t)+epsilon^2*x[2](t), eq);

pls help with this...