Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

how to singularize existing function?

any group theory can describe the singularity function?

Good day, please how can one solve these BVPs using FINITE DIFFERENCE METHOD in maple. Here is the problem FDM.mw

Hi;

Could anyone help me to convert the following codes into Mpale?

while n>= 1
    for i=1:n
        S(i)=resultant(F(i), F(i+1), Plex(i));
    end
    F(i)=S(i);
    n=n-1;
end

Thanks in advance.

Hello! I have written a code to solve 1D heat equation with Neumann B.C. using explicit method. But there is a problem with plot. I don't understant, what's wrong. Help me, please. I will be very grateful. The code is attached here.
pr3.mw
11.pdf

Hello,

I am writing a program in C that uses the open maple library. It is not the first time that I use it but now I am facing a strange problem that involves the simplify command: suppose a,x,y are symbols that are not previously used in maple, the following lines

1)  EvalMapleStatement(kv, "simplify((a*x^2-y^2)/(x^2*y^2-1));");

2)  EvalMapleStatement(kv, "simplify((2*x^2-y^2)/(x^2*y^2-1));");

only differ by the fact that the parameter a is replaced by 2 in the second line. But they return the following output:

1) (a*x^2-y^2)/(x^2*y^2-1)

that is correct, nothing to simplify..

2) Error, (in gcd/LinZip) input must be polynomials over the integers

I must be doing something wrong but I am getting nowhere...

Thanks...

 

P.S. This is the complete listing 

 

#include <stdio.h>

#include <stdlib.h>

 

#include "maplec.h"

 

static void M_DECL textCallBack( void *data, int tag, char *output )

{

    printf("%s\n",output);

}

 

int main( int argc, char *argv[] )

{

    char err[2048];  /* command input and error string buffers */

    MKernelVector kv;  /* Maple kernel handle */

    MCallBackVectorDesc cb = {  textCallBack,

                                0,   /* errorCallBack not used */

                                0,   /* statusCallBack not used */

                                0,   /* readLineCallBack not used */

                                0,   /* redirectCallBack not used */

                                0,   /* streamCallBack not used */

                                0,   /* queryInterrupt not used */

                                0    /* callBackCallBack not used */

                            };

    ALGEB r, l;  /* Maple data-structures */

        char *myargv[]={"maple"};

        int myargc=1;

 

    if( (kv=StartMaple(myargc,myargv,&cb,NULL,NULL,err)) == NULL ) {

        printf("Fatal error, %s\n",err);

        return( 1 );

    }

 

EvalMapleStatement(kv, "simplify((a*x^2-y^2)/(x^2*y^2-1));");

 

EvalMapleStatement(kv, "simplify((2*x^2-y^2)/(x^2*y^2-1));");

 

    StopMaple(kv);

 

    return( 0 );

}

 

compiled with

gcc prova.c -I /Library/Frameworks/Maple.framework/Versions/Current/extern/include/ -L /Library/Frameworks/Maple.framework/Versions/Current/bin.APPLE_UNIVERSAL_OSX/ -l maplec

Hello i want to solve the differentiel equation but i have these problem i don't understund  why !?

``

``

restart:with(plots):

mb:=765 : mp:=587 :Ib:=76.3*10^3 :Ip:=7.3*10^3 :l:=0.92 :d:=10:F:=0.042:omega:=0.56 :

eq1:=(mb+mp)*diff(x(t),t$2)+mp*l*cos(alpha(t))*diff(alpha(t),t$2)+mp*[l*[diff(alpha(t),t)]^2*sin(alpha(t))]-F*sin(omega*t)=0;

1352*(diff(diff(x(t), t), t))+540.04*cos(alpha(t))*(diff(diff(alpha(t), t), t))+[540.04*[diff(alpha(t), t)]^2*sin(alpha(t))]-0.42e-1*sin(.56*t) = 0

(1)

eq2:=-mp*l*sin(alpha(t))*diff(alpha(t),t$2)+mp*[l*[diff(alpha(t),t)]^2*cos(alpha(t))]-9.81*(mp+mb)-F*sin(omega*t)=0;

-540.04*sin(alpha(t))*(diff(diff(alpha(t), t), t))+[540.04*[diff(alpha(t), t)]^2*cos(alpha(t))]-13263.12-0.42e-1*sin(.56*t) = 0

(2)

eq3:=mp*[d+l*cos(alpha(t))]*diff(x(t),t$2)+[Ip+mp*l^2+mp*d*l*cos(alpha(t))]*diff(alpha(t),t$2)-mp*sin(alpha(t))*[-l*d*alpha(t)^2]+mp*[l*9.81*sin(alpha(t))]=0;

587*[10+.92*cos(alpha(t))]*(diff(diff(x(t), t), t))+[7796.8368+5400.40*cos(alpha(t))]*(diff(diff(alpha(t), t), t))-587*sin(alpha(t))*[-9.20*alpha(t)^2]+[5297.7924*sin(alpha(t))] = 0

(3)

eq4:=mp*l*cos(alpha(t))*diff(x(t),t$2)+(Ip+mp*l^2)*diff(alpha(t),t$2)-mp*9.81*l*sin(alpha(t))=0;

540.04*cos(alpha(t))*(diff(diff(x(t), t), t))+7796.8368*(diff(diff(alpha(t), t), t))-5297.7924*sin(alpha(t)) = 0

(4)

CI:= x(0)=0,alpha(0)=0,D(x)(0)=0,D(alpha)(0)=0;

x(0) = 0, alpha(0) = 0, (D(x))(0) = 0, (D(alpha))(0) = 0

(5)

sys := eq1, eq2, eq3, eq4:

``

solution:=dsolve([sys,CI],numeric);

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

 

 

NULL


thanks for your help

Download tangage.mw

Is it possible to numerically calculate  the integral

int((-12*y^2+1)*ln(abs(Zeta(x+I*y)))/(4*y^2+1)^3, [y = 0 .. infinity, x = 1/2 .. infinity])

in Maple?

The code

int((-12*y^2+1)*ln(abs(Zeta(x+I*y)))/(4*y^2+1)^3, [y = 0 .. infinity, x = 1/2 .. infinity],numeric,epsilon=0.1)

has been executed on my comp  without any output since this morning.

 

 

 

Hello everybody, i need to graphic a couple of functions just like this one:

 http://temasmatematicos.uniandes.edu.co/Casquetes_cilindricos/Pags/Anim_1.htm

i have been watching this:

 http://www.maplesoft.com/teachingconcepts/detail.aspx?cid=12 VISUALIZATION --> Animation 2

i've tried with (plots) (plottools) animate, etc. but i can't figure out how to do it. 

It would be very helpful if someone explain me how to do this.

Thank you all!

I have setup two groups in MapleCloud, to be used for distribution of materials (a readonly group) and as a drop for students to deposit their exercise sheets (an opaque group). It is my understanding that e.g. students have to request joining a group and I get to approve them.

How does one join a group managed by someone else? If I try from another account I cannot see the groups I just created, so I cannot request joining them. I seem also to not be able to just approve someone (by their email e.g.) without that person having requested to join first.

???

Thanks

USPAS2014

Is there a way of ploting y=x^x for when x<0 and y is a noncomplex solution?

 

I have some questions about Mobius Project:

  1. What is the license of the worksheets posted on Mobius Project's website?
    I can not find any information about the license of the materials submitted to the Mobius Project.
    Is my worksheets are protected by intellectual property laws and copyrighted by me?
    Or it is freely available for any use for everyone?
  2. What are the terms of use of Mobius Project and its service? I can not find it anywhere.
  3. What is the copyright for http://mobius.maplesoft.com/?
    I do not see any copyright in the footer of the site or elsewhere.
    So, no copyrights mean the materials/resources on this website can be freely copied/used by anyone?

with pointplot3d and 14,000 points when I enter symbol=point I get an empty plot.

Only when I set symbolsize=1 (a point) do I get points appearing in the graph.  Bug?

Dear Users

I have a problem for solving a system of linear equations that arise from collocation method for getting approximate solution of a coupled PDE and ODE in Food engineering problems.

When it reach to the fsolve command it takes long time!!!

I used maple 13.

If kindly is possible, please help me in this special case.

With kind regards,

Emran Tohidi.

 

> Restart;
print(`output redirected...`); # input placeholder
> h := 50; hm := 0.1e-3; rhodp := 1500; Y := 0.5e-1; T0 := 20; rhoair := 1.2041; Dair := 0.2e-8; DD := 0.85e-9; C := 3240; L := 0.4e-1; X0 := 1.5; V := .2; delta := 0.2e-2; Yair := 0.5e-1; nu := .2; Tair := 60; Hnu := 2400; rho := 1359; tt := 3;
%;
> N := 5; Digits := 20;
> X := sum(sum(a[m, n]*z^m*t^n, m = 0 .. N), n = 0 .. N); X := unapply(X, z, t); Xt := diff(X(z, t), `$`(t, 1)); Xt := unapply(Xt, z, t); Xz := diff(X(z, t), `$`(z, 1)); Xz := unapply(Xz, z, t); Xzz := diff(X(z, t), `$`(z, 2)); Xzz := unapply(Xzz, z, t); T := sum(b[n]*t^n, n = 0 .. (N+1)^2-1); T := unapply(T, t); Tt := diff(T(t), `$`(t, 1)); Tt := unapply(Tt, t); aw := exp(.914)*X(z, t)^.5639-.5*exp(1.828)*X(z, t)^(2*.5639); aw := unapply(aw, z, t); TT := 8.3036+3816.44*(1+T(t)/(46.13)+T(t)^2/46.13^2)/(46.13); TT := unapply(TT, t); pwv := 133.3*(1+TT(t)+(1/2)*TT(t)^2); pwv := unapply(pwv, t); Yi := .622*pwv(t)*aw(z, t)*(1+pwv(t)*aw(z, t)/rho+(pwv(t)*aw(z, t)/rho)^2)/rho; Yi := unapply(Yi, z, t);
%;
> S1 := {seq(seq(Xt(delta*i/N, tt*j/N)-DD*Xzz(delta*i/N, tt*j/N) = 0, i = 1 .. N-1), j = 1 .. N)};
> S2 := {seq(DD*rhodp*Xz(delta, tt*j/N)+hm*rhoair*Yi(delta, tt*j/N) = 0, j = 0 .. N)};
> S3 := {seq(Xz(0, tt*j/N) = 0, j = 0 .. N)};
> S4 := {seq(X(delta*i/N, 0) = 0, i = 1 .. N-1)};
> S5 := {seq(seq(rho*delta*C*Tt(tt*j/N)-h*(Tair-T(tt*j/N))+hm*Hnu*rhoair*(Yair-Yi(delta*i/N, tt*j/N)) = 0, j = 1 .. N), i = 0 .. N)};
print(`output redirected...`); # input placeholder
> S6 := {seq(rho*delta*C*Tt(0)-h*(Tair-T0)+hm*Hnu*rhoair*(Yair-Yi(delta*i/N, 0)) = 0, i = 0 .. N)};
%;
> SS := `union`(`union`(`union`(`union`(`union`(S1, S2), S3), S4), S5), S6);
> sol := fsolve(SS);

Hello, I have two plots:

Is it possible to display differene between these plots?

I don't understand how to create better parameters. Why is the graph not flush with the image (it appears that the graph is a set of points and NOT lines)

 

First 1343 1344 1345 1346 1347 1348 1349 Last Page 1345 of 2224