MaplePrimes Questions

Hey, 

I am trying to write a single procedure to find the root of any function using the Newton-Raphson method, given the initial approximation and the tolerance. If this fails to converge, the program must then use the Bisection method to find the root. Need some help please. The current procedure i have done is only coming out with the first Iteration 

Thanks for the help!

hi dears. i cant solve this question by Analytical method such variational itteration method (VIM) or homotopy perturbation method (HPM), homotopy analysis method(HAM)

Hello;

Suppose we have a vector as follows:

Plex:=[u,v,w,x,y,z]

Could anyone help me to reproduce the following vectors via a "for" loop?

[ v, u, w, x, y, z]

 [ w, v, u, x, y, z]

 [ x, w, v, u, y, z]

 [ y, x, w, v, u, z]

I have a simple Matlab script which can do what I intend. So, my question is how can I implement this script in Maple?

Plex=[u,v,w,x,y,z]

for i=2:5

    Plex=[Plex(i) Plex];

    Plex(i+1)=[];

    disp(Plex)

end

 

Thanks in advance.

Hello all,

I'm an engineering student working on bicycle shifting. After an dynamic study i encounter the equation:

E:=5.37*theta12(t)=C*(diff(theta12(t),[t$2])+diff(theta10(t),[t$2]))+247.2*[diff(theta10(t),[t$2])*(-6.53*cos(theta12(t))+8.51*sin(theta12(t)))+diff(theta10(t),t)*diff(theta12(t),t)*(6.53*sin(theta12(t))+8.51*cos(theta12(t)))]-247.2*diff(theta10(t),t)*(-diff(theta12(t),t)+diff(theta10(t),t))*(-6,53*sin(theta12(t))-8,51*cos(theta12(t)));


where C is a real constant

More exactly, i would need to find a formule for theta12 using the succesive derivatives of theta10.

I heard that i need to search a numerical solution but,up to now, don't succeed to solve it.

I'm open to all idea.

Thanks

I have low vision, and the default fonts in Maple are too small for me to read.  I can enlarge the normal input and output using the View/Zoom Factor, and this works very well for me.  However, it seems to have no effect on the fonts in spreadsheets.

Can anyone tell my how to make the fonts in Maple spreadsheets larger?  I'm using Maple 18 on Windows 8.1. Thank you.

how to singularize existing function?

any group theory can describe the singularity function?

Hi, My goal is to compute the coefficient beta_i, so i will solve a system and get the coefficient beta_i. But my code return an error. Any help please. Many thinks

coef_approx:=proc(a,N,i,d)
local Fredholm,eq2,eq3,Vct_basis,fct,sys,eq4,M,w,b,M1,V,Vect_beta,h,x,phi,Kernel,lambda;
# Fredholm Integral equation
Fredholm:=phi(x)=f(x)+lambda*int(Kernel(x,y)*phi(y),y=-a..a);
# stepsize
h:=a/N;
# First Approximation of integral
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 (Method used)
eq3:=phi->int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h)=add(beta[i]*phi((n+i-1)*h),i=1..d+1):
eq4:=int(Kernel(x,y)*phi(y),y=n*h..(n+d)*h)=add(alpha[i](n,m)*phi((n+i-1)*h),i=1..d+1);
# Fct used to compute the coeffcient beta[i]
Vct_basis:=[seq(x^i,i=0..d+1)]:
fct:=[seq(unapply(Vct_basis[i],x),i=1..d+2)];
# system of equation must be solved
sys:=[seq(eq3(fct[i]),i=1..d+1)]:
x:='x';
x:=m*h:
w := [seq(beta[i],i=1..d+1)];
M,b := GenerateMatrix(sys,w);
M1:=-M: V:=-b:
Vect_beta:=(M1)^(-1).V:
return Vect_beta;
end proc;

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

First 1420 1421 1422 1423 1424 1425 1426 Last Page 1422 of 2427