Items tagged with error error Tagged Items Feed

Hi all,

A while ago, I had asked a question relating to calling to an MapleSim .msim simulation within maple .mw worksheet. (Thank you Acer! :) I got it working but then another issue just popped up...The basic structure is as follows:

********************************************************************

A := MapleSim:-LinkModel('filename' = "TriplePendulum_WORKS_v2.msim")

for z from 1 to 50 do

     initialize_variables;

Hey everyone.

I have been trying to solve an ODE with BC with dsolve,numeric but always get the error below:

> Error, (in dsolve/numeric/BVPSolve) matrix is singular

My problem is not solving it (ie fine to do with LinearAlgebra package or similar), but why does this error occur? What is the problem with the ODE/BC that produces this error?

And is there a way of adjusting the ODE/BC to be able to solve it simply by using dsolve,numeric?

On my Dell Latitude E5400 laptop with up-to-date Fedora 16:

$ maple15
    |\^/|     Maple 15 (X86 64 LINUX)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2011
 \  MAPLE  /  All rights reserved. Maple is a trademark of
 <____ ____>  Waterloo Maple Inc.
      |       Type ? for help.
>  log[2](8);

I tried to solve a non linear coupled boundary value problem in MAPLE using DSolve command. The code is :

alias(eta = e, theta = t)

Eq[1] := 5*(diff(F(e), `$`(e, 3)))+(m+3)*F(e)*(diff(F(e), `$`(e, 2)))-(2*m+1)*(diff(F(e), e))^2-(4*m+2)*H(e)-(m-2)*e*(diff(H(e), e)) = 0

Eq[2] := diff(H(e), e) = t(e)

Eq[3] := 5*(diff(t(e), `$`(e, 2)))/Pr-(m+3)*F(e)*(diff(t(e), e))-5*m*(diff(F(e), e))*t(e) = 0

BCs := [F(0), (D(F))(0), (D(F))(infinity), t(0)-1, t(infinity), H(infinity)]

Postage_Stamp_proble.mw

Hi

I have a program which I need to exit a loop when an error is thrown. it should then output the value of a. but it doesn't....

I want to get the zeroth order and first order perturbation equations of the following nonlinear coupled ODE.

restart; N := 1

alias(eta = e, theta = t)

eq[1] := 5*(diff(F(e), `$`(e, 3)))+(1+3/m)*F(e)*(diff(F(e), `$`(e, 2)))-(2+1/m)*(diff(F(e), e))^2-(4+2/m)*H(e)-(1-2/m)*e*(diff(H(e), e)) = 0

eq[2] := diff(H(e), e) = t(e)

eq[3] := 5*(diff(t(e), `$`(e, 2)))/Pr+(1+3/m)*F(e)*(diff(t(e), e))-5*(diff(F(e), e))*t(e)

eqs := eq[1], eq[2], eq[3]

I wrote the following MAPLE code

restart; N := 2

alias(epsilon = e, eta = t)

eq1 := 5*(diff(F(t), `$`(t, 3)))+(1+3*e)*F(t)*(diff(F(t), `$`(t, 2)))-(2+e)*(diff(F(t), t))^2 = 0

F := proc (t) options operator, arrow; sum('F[i](t)*e^i', 'i' = 0 .. N) end proc

deqn := simplify(collect(eq1, e), {e^(N+1) = 0})

MAPLE says

Error, (in sum) too many levels of recursion

How do I get this error corrected?

I asked a question few hours ago. I corrected the error in the process I was suggested. How in writing the program to collect the coefficients of m^0, m^1 and m^2, I am again stuck in an error and could not find any help.

restart; N := 2

Eq[1] := diff(f(eta), `$`(eta, 3))+(m+4)*f(eta)*(diff(f(eta), `$`(eta, 2)))-(2*(m+1))*(diff(f(eta), eta))^2-(2*(m+1))*h(eta)-(.5*(m-2))*eta*(diff(h(eta), eta))

Eq[2] := g(eta)-(diff(h(eta), eta))

Eq[3] := (diff(g(eta...

I am trying to evaluate integral (-1/rho)*del (p)/del (x) dy within the limits y=0 to y=delta (x). For that I first evaluated (-1/rho)*del (p)/del (x) and stored it in res2. Then I tried to integrate it w.r.t dy within limits 0 to delta (x). I wrote the following Maple lines for the purpose:

Te := proc (x, y) options operator, arrow; T[inf]+(T[w]-T[inf])*(1-y/delta[t](x))^2/(1+2*A*H/delta[t](x)) end proc

PDE := diff(P(x, y), y)-rho*g*beta*(Te(x, y)-T[inf]) = 0

restart

u := proc (x, y) options operator, arrow; u[x](x)*(y/delta(x)-2*y^2/delta(x)^2+y^3/delta(x)^3+B*H*(1-3*y^2/delta(x)^2+2*y^3/delta(x)^3)/delta(x))/(3*B*H/delta(x)+2) end proc

int(u(x, y)^2, y = 0 .. delta(x))

I tried the following partial integration in MAPLE. But Maple says:

Error, (in u[x]) too many levels of recursion.

Analytically it is possible to evaluate the integral in the above lomits. Where lies the problem?

I wrote the following MAPLE lines to determine the three constants a, b and c from three boundary conditions. Whatever a, b and c I get as result are correct. Now I want to substitute back the values of obtained a, b and c in the expression for Te to obtain the resulting temperature distribution. I used both algsubs(a, b, c, Te) and subs(a, b, c, Te) one by one to substitute a, b and c in the expression for Te and get the result. But both way it is showing...

I wrote the following lines in MAPLE:

This is the procedure for the question I posted on http://www.mapleprimes.com/questions/128600-Romberg-Approximation-Procedure-In-Maple?sq=128600 

Can anyone see where the problem is?

RombInt:=proc(f,a,b,m,n)
local j,k,R;
if type(m, numeric) then
   if not type(m,nonnegint) then ERROR("m must be a non-negative integer"); fi;fi;

Hi,

I wanted to know if it possible to get maple to skip all operations to the bottom if it receives an error anyway through my program? I have a number of seperate lines and don't want to put it into one big if statement.

Thanks

This is the Pde I am trying to solve, with the initial profile function at t=0 called Piinitial, and IC=initialconditions,boundary conditions,

Piinitial := piecewise(r <= 1, (1/9)*r^2+(1/12150)*(31*(8807640902+1170*sqrt(24475122685023))^(2/3)+109499254+22553*(8807640902+1170*sqrt(24475122685023))^(1/3))/(8807640902+1170*sqrt(24475122685023))^(1/3), r <= 2 and 1 < r, (4/405)*(8807640902+1170*sqrt(24475122685023))^(1/3)*r+(14128936/405)*r/(8807640902+1170*sqrt(24475122685023...