## How to plot heat transfer rate for this problem?...

HTR.mw

In above problem, Additionally How to  plot  heat transfer rate  Q versus L^2  for distinct porosity parmeters(Sh) , using  heat transfer rate formula, Q = (q*L)/(k*A*T[b])=theta'(1).

using  [Sh = 0.1, L^2 = 0.1, Nr =0 .1, Ha =0 .1, Pe = 0.1],  [Sh = 0.3, L^2 = 0.3, Nr = 0.1, Ha = 0.1, Pe =0 .1],   [Sh = 0.5, L^2 =0 .5, Nr =0 .1, Ha = 0.1, Pe =0 .1].

## How can I find a values of the coefficients in pol...

I tried to solve a equation f(x). Finally I got a polynomial with coefficients.  (see Solve_for_coefficients.mw)

I want to find the values of unknown coeffients in the last polynomial equation as procedure to equating each cofficients to zero.

For example, the coefficients of x^2 is 1/48*(-24*p[2]+48*q[1]+24*q[2]+24) = 0.

Likewise, collect all the coefficients from x^0 to x^4 and solve them. Finally we get the coefficients value.

My question is How to collect the coefficients in this polynomial and solve this coefficients?

## How to add vector fields ...

Hallo every body

How to add vector fields to the figure of this example of a three-dimensional differential system.

in maple 18

Porgram_of_corollary_1_in_Maple.mw

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)

Considérons le système suivant:

 >
 >
 >
 >
 >
 >
 (5)
 >
 (6)
 >
 (7)

Faisons le changement (x,y,z)=(εX,εY,εZ)

 >
 (8)
 >
 (9)
 >
 (10)
 >
 (11)

Faisons le changement

 >
 (12)
 >
 (13)
 >
 (14)
 >
 (15)

Utilisons le développpement de taylor

 >
 (16)
 >
 (17)
 >

Averaging d'ordre 1

Les fonctions F11 et F21 sont données comme suit:

 >
 >
 (18)
 >
 (19)
 >

Calculons les fonctions moyennées et

 >
 (20)
 >
 (21)
 >
 (22)
 >

Averaging d'ordre 2

 >
 >
 (23)
 >
 (24)
 >

Calculons les fonctions moyennées et

 >
 (25)
 >
 (26)
 >
 (27)
 >
 (28)
 >
 (29)
 >
 >
 >
 (30)
 >
 (31)
 >
 (32)
 >
 (33)
 >
 (34)
 >
 (35)
 >
 >
 >
 (36)
 >
 (37)
 >
 (38)
 >
 (39)
 >

## Error, (in dsolve/numeric/process_input) input sys...

Dears,

Can you look the code bellow and send me my error please? I used Maple 18.

restart;
with(plots);
theta(t) = 19.592+1.2697*cos(.5240*t+4.3391)-.6343*cos((2*.5240)*t-.6963);
omicron(t) = 99.4876+89.8581*cos(.5232*t+15.4500)+19.1069*sin((2*.5232)*t)-8.5891*cos((3*.5232)*t+3.7723)+6.4660*sin((5*.5232)*t);
`&varphi;`(theta):=0.000203*theta*(theta - 11.7)*sqrt(42.3-theta);
mu[v](theta,omicron):=0.0886*exp(((0.01*omicron +1.01*theta  -21.211)/(14.852))^(4));
p[0](theta):=(-0.153* theta*theta + 8.61*theta - 97)/(mu[v](theta,omicron)):
p[2](omicron):=(4*0.25)/(2500)*omicron*( 50 -omicron);
p[3](omicron):=(4*0.75)/(2500)*omicron*( 50 -omicron);
p[2](theta):=exp (0.06737 - 0.00554*theta);
theta[EA](theta):=1/(-0.00094*theta*theta + 0.049*theta - 0.552);
L[v](theta,omicron):=(3.375*(4*omicron*(50-omicron))^(3)*exp(0.0054*theta+0.6737))/(50^(6)*(2+(0.00554*theta-0.06737)^(-1)));
eta(theta,omicron):=(p[0](theta)*p[1](omicron)*p[2](omicron)*p[3](omicron)*p[2](theta))/(theta[EA](theta));
lambda[v] := beta[v]*`&varphi;`(theta)*i[v](t)/n[h](t);
lambda[h] := (beta[h]*`&varphi;`(theta)*i[h](t)+beta[h]*`&varphi;`(theta)*omega*r[h](t))/n[h](t);
n[v](t):=s[v](t)+i[v](t);
beta[h] := 0.9e-1; beta[v] := 0.2e-1; Lambda[h] := .50; sigma[1] := 0.15e-1; sigma[2] := 0.71e-1; Omega[h] := .50; mu[h] := 0.128e-1; delta[h] := .45; k[v] := .66; omega := .3; mu[d] := 0.14e-2;
sys := {diff(i[h](t), t) = Omega[h]+sigma[2]*r[h](t)+lambda[v]*s[h](t)-(delta[h]+mu[d]+mu[h])*i[h](t), diff(i[v](t), t) = lambda[h]*s[v](t)-mu[v](theta, omicron)*i[v](t), diff(j[v](t), t) = L[v](theta, omicron)*(1-j[v](t)/k[v])*n[v](t)-(eta(theta, omicron)+mu[j](theta, omicron))*j[v](t), diff(r[h](t), t) = delta[h]*i[h](t)-(sigma[1]+sigma[2]+mu[h])*r[h](t), diff(s[h](t), t) = Lambda[h]+sigma[1]*r[h](t)-(lambda[v]+mu[h])*s[h](t), diff(s[v](t), t) = eta(theta, omicron)*j[v](t)-(lambda[h]+mu[v](theta, omicron))*s[v](t), i[h](0) = 100, i[v](0) = 100, j[v](0) = 200, r[h](0) = 0, s[h](0) = 10000, s[v](0) = 5000};
p1 := dsolve(sys, numeric, method = rkf45, output = procedurelist);
Error, (in dsolve/numeric/process_input) input system must be an ODE system, found {mu[j](theta, omicron), mu[v](theta, omicron)}
p1o := odeplot(p1, [theta, omicron, i[h](t)], 0 .. 10, numpoints = 100, labels = ["Time (Days)", " infectious population"], labeldirections = [horizontal, vertical], style = line, color = red, axes = boxed, legend = [front, rear, ideal]);

## How can I solve a higher order ODE to first ODE?...

I tried to solve a hogher order ODE system to reduce first order ODE sytem. In this case I have assume a initial conditions as a variable. (sse the attachment, Conver_to_first_order.mw).

I got an error in specification of initial value...

How can I solve this problem in RKF45 method with shooting technique.

Here, I have assume f=x[1], f'=x[2]. f''=x[3], f'''=x[3]' and theta=x[4], theta'=x[5], theta''=x[5]'

## How do i find the solution or integral curve given...

Determine the solution or integral curve given by the differential equation x+yy'=0 by the method of isodine. hence display its 2Dplot

## How can I solve a system of ODE in rkf45 method wi...

I tried to solve a ODE system using rkf45 with shooting technique. But I have a lot of errors. How to resolve this...

See the attachment: Shoot_Blasius.mw

## How do I solve a ODE system using Newton's finite ...

How to solve a ODE system using Newton's finite difference method?
How to plot f', theta and phi functions for various alpha values.

Difference.mw

## How can I get a ploting values in data?...

How to convert a ploting values in a graph to excel?

restart;
with(PDETools): with(DETools): with(plots):with(plottools):
eq1 := ((D@@2)(f))(eta)*f(eta)*sin(alpha)+((D@@2)(f))(eta)*eta*cos(alpha)+2*((D@@3)(f))(eta) = 0;
ics := f(0) = 0, (D(f))(0) = 0, (D(f))(10) = 1; bcs := (D(f))(10) = 0, theta(10) = 0, phi(10) = 0;
Parameters1 := alpha = (1/3)*Pi;
sol1 := dsolve(eval({eq1, ics}, {Parameters1}), numeric);
p1 := odeplot(sol1, [[eta, ((D@@2)(f))(eta)]], eta = 0 .. 10, color = [red], axes = boxed);
display({p1});

## How to perform good quality graphics export......

Hello everyone,

I'm using Maple 18. I have a problem that I can't solve. It concerns the programmatic - from the Maple code level - export of animations (even single images) with good quality graphics, i.e. either high resolution or "large" size - which translates into the same. This is because Maple 18 does not have a size option in plot3d.

Even manually enlarged graphics - a bit too much, makes manual export impossible; Maple reports "..the file could not be created...". Most often, this ends with the creation of a gif file with a capacity of 0 bytes.

I am asking for help if anyone knows how to deal with this [normal export of "normal" graphics manually or with software can of course be done in Maple 18].

Regards.

## Classic DiffGeo - without using Maple's built-in...

by: Maple 18

Hello,

Attached I am sending several procedures for curves in 3D space. They were written without using Maple's built-in DiffGeo procedures and functions. As an example of their use, I made several animations - Maple worksheets are attached. I hope that maybe they will be useful to someone.

Regards.

ClsDGproc-Curves.zip

## How do I solve a Blasius problem in Shooting metho...

I tried to solve a Blasius problem (available in maple), but I have an error. How to solve this issue.

## how can i solve elliptic integral...

p1>0, p2<0 and x[0] are arbitrary constants. How can i solve this integral?

I want to see

## How do I solve ODE system in RK method?...

How to solve and plot a ODE system in RK method.
eq1 := diff(f(x), x, x, x)-(1/2)*Sc*sin(alpha)*g(x)*(diff(g(x), x, x))+(1/2)*x*cos(alpha)*(diff(f(x), x, x))+(1/2)*sin(alpha)*f(x)*(diff(f(x), x, x)) = 0; eq2 := (diff(g(x), x, x, x))/Pm+(1/2)*x*cos(alpha)*(diff(g(x), x, x))+sin(alpha)*f(x)*(diff(g(x), x, x))-sin(alpha)*(diff(f(x), x, x))*g(x) = 0; eq3 := (diff(theta(x), x, x))/Pr+(1/2)*x*cos(alpha)*(diff(theta(x), x))+(1/2)*x*(diff(f(x), x))*(diff(theta(x), x))+sin(alpha)*(x*(diff(f(x), x))-f(x))*(diff(theta(x), x))-Nb*(diff(s(x), x))*(diff(theta(x), x))-Nt*(diff(theta(x), x))^2+(1/4)*Sc*Br*sin(alpha)^2*(diff(f(x), x))^2*(x*(diff(g(x), x))-g(x))+(diff(g(x), x))^2*(x*(diff(f(x), x))-f(x)) = 0; eq4 := diff(s(x), x, x)+S*((1/2)*cos(alpha)*x*(diff(s(x), x))+(1/2)*sin(alpha)*f(x)*(diff(s(x), x)))+Nt*(diff(theta(x), x, x))/Nb = 0

ics := f(0) = 0, (D(f))(0) = 1, g(0) = 0, (D(g))(0) = 1, theta(0) = 1, s(0) = 1; bcs := (D(f))(100) = 0, (D(g))(100) = 0, theta(100) = 0, s(100) = 0

alpha = - 30 degree, Sc = 1.0, Pm = .1, Pr = 6.2, Nb = .1, Nt = .1, Br = .5, S = 1

## Nested for do loops...

Hello,

How to  outer loop index setup as the upper value of the inner loop index?

 1 2 3 4 5 6 7 Last Page 1 of 87
﻿