find Databases for the NASA Seven-Coefficient Poly...

Hi.

in the ThermophysicalData[Chemicals] package that compute the coefficients for different species how I can find that coefficients for seven coefficients not nine of them

in other words, I am seeking to find Databases for the NASA Seven-Coefficient Polynomial Fits for Calculating Thermodynamic Properties of Individual Species.

Best

Q1: how to solve nonlinear simultaneous equations...

hello,

My equations are not getting solved.

Please tell me, is there any other command to solve these type of equations.

Q1.mw

Why is plotting my function giving a different val...

I have an issue where I believe the values given by evalf (I have independent verification), but the plot of the same function is incorrect. Specifically, g(a) plotted for a <=0 is correct, and incorrect for a > 0.

What am I doing wrong?

g:=(a)->int(exp(-(x^2-a)^2),x=0..infinity) plot(g(a), a=-2..5) #I don't believe anything for a > 0

evalf(g(-1)) #I believe all of these values                           0.2059311656 evalf(g(0))                           0.9064024772 evalf(g(1))                           0.9868660750

Mattthew

ONE DIMENSIONAL HEAT EQUATION-PROBLEM...

DEAR SIR,

I REQUEST YOU TO PROVIDE THE METHOD TO SOLVE THE FOLLOWING PROBLEM.

THANK YOU

The Fitting Procedure...

The Fitting Procedure

The data are T,CP,H and S. I need to determine a polynomial with seven Coefficients.

The essential input to the fitting procedure is a table of specific heat, enthalpy, and entropy as a function of temperature.

The constrained linear least-squares fitting procedure is a three-stage process. The first step is to determine simultaneously a1, through a5 for temperature ranges by fitting the specific heats to Eq. (1).

The second step is to determine a6 for temperature ranges by fitting the H data. In this state, a1-a5 are held fixed, and as in the cp fitting, equality constraints on the fit and its first derivative are imposed at the common temperature. Finally, the a7 coefficients are determined by fitting the entropy (S) data.

R=1.987

DATA.txt

FITTING.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/FITTING.mw .

no plot obtained!!!...

After running the program without any error, no plot obtained.

What is the reason?

Thanks3.mw

DeepLearning[Classify] systematically returns conn...

Hi,

I try using the DeepLearning package.
I use the function Classify and, even in the simplest test case presented in the its help page (please look at it), I regularly get connection errors to the mpython server as soon as I execute classifier := Classify(...) or classifier(...) more than once.
Errors are one of these twos

Error, (in Train) unable to communicate with mpython server
or
Error (in Python:-EvalFunction) unable to communicate with mpython server

I work with Windows 7 Enterprise, on an 8 proc PC and 64 GB of memory. The worst situation happened when Maple didn't even return these errors and that I saw inflating the consumed memory in 2 minutes, forcing me to manually shut down my PC because the task manager wasn't no longer  operational.

Is it a known problem?
Could it be an installation problem?

Even if it's not the point here, I would like to say that trying to use the DeepLearning package is really challenging considering the poverty of the help pages.

Solving an optimal control problem...

Hello all

I wanna solve an optimal control problem and I have searched the Internet but I could not find any tutorial or video course on how to solve it with the Pontryagin maximum principle method. It is my first time that I want to use MAPLE for solving an optimal control problem and I would be thankful if someone can help me.

$$\max \int_{0}^{1} x_{2} [u(t)-u(t)^2] dt$$

$$\dot{x}_{0} = -(1-u(t)) x_{0}(t)+2 x_{1}(t)$$

$$\dot{x}_{1}(t) = (1-u_{t}) x_{0}(t) +2 x_{2}(t) -[3-u(t)]x_{1}(t)$$

$$\dot{x}_{2}(t) = (1-u(t))x_{1}(t) -2 x_{2}(t)$$

$$0 \le u_{t} \le \frac{1+t^2}{1+t}$$

Thanks

I understand this is asking for homework help but ...

Moving to online learning has proved difficult to get any real 1 on 1 time or with my teachers around this question, and the help I have got hasn't helped me as yet.

I understand that asking for homework help is probably frowned upon but I would appreciate any help, guidance or direction on how to answer the below in Maple2019.

The general formula for a plane in the 3D space is z = ax + by + c, where a, b, c are the parameters. Alice encodes three English words to three numbers by using the Maple command "text2num". She then set up a (3,4) secret sharing scheme with the idea of the Blakley method. The four shares are four planes given by:

z =   4x + 19y + 2515211725275120 (mod 2515211819051461)

z = 52x + 27y + 2515210613496048 (mod 2515211819051461)

z = 36x + 65y + 2515210981587340 (mod 2515211819051461)

z =   6x + 60y + 2515211676449260 (mod 2515211819051461).

Find the secret English words of Alice with Maple command

I'm absolutely hitting a wall. The most I have found in the word "you".

Creating a gradient with one level curve...

I have been stuggling with this for a few hours. I have my function 4x^2+9y^2. I have my point (2,1,25). I found my gradient to be 16i+18j. I want to plot the level curve that contains the point (2,1,25) and the gradient together.

I know how to graph multiple level curves but not one.

3D Graph of Parabloid...

I am so sorry for all of the questions. I am getting extremely frustrated with Maple.

I want to plot a 3D rendering of 4x^2+9y^2 but I want the full paraboloid to show. However, all I keep getting is cut off versions.

factoring with exponentials...

I would like to factor an algebraic expression that contains exponentials.
Maple spits it back unfactored. Did I enter it wrong?

factor(exp(4*y) + 2*exp(2*y) + 1);

Why is this happening?...

I want the legend to be produced to look like the third to last example. https://fr.maplesoft.com/support/help/Maple/view.aspx?path=plots/contourplot

However, mine keeps coming out like this.

can any one solve this?...

 >
 > restart:
 >
 (1)
 > eq11 := (10-9*lambda1)*n*diff(f(eta),eta,eta,eta)*(-diff(f(eta),eta,eta))^(n-1)-m*diff(f(eta),eta)^2+M^2*diff(f(eta),eta)+(m*(2*n-1)+1)/(n+1)*f(eta)*diff(f(eta),eta,eta)=0;
 (2)
 > eq21 := (10-9*lambda1)*n*diff(theta(eta),eta,eta)*(-diff(theta(eta),eta))^(n-1)+(m*(2*n-1)+1)/(n+1)*f(eta)*diff(theta(eta),eta)-m*diff(f(eta), eta)*theta(eta)=0;
 (3)
 > bc21:= f(0)=0, D(f)(0)=1+alpha*(D@@2)(f)(0), D(f)(N)=0, theta(0)=1+b*D(theta)(0), theta(N)=0;
 (4)
 >
 > n:=1.2: b:=0: M:=0: alpha:=0:
 > N:=6:
 >
 > m:=0.8:
 > A11:=dsolve({eq11,eq21,bc21},numeric,continuation=lambda1,maxmesh=15000, output=array([seq( i, i=0..N,0.01 )])):
 > m:=1.2:
 > A12:=dsolve({eq11,eq21,bc21},numeric,continuation=lambda1,maxmesh=15000, output=array([seq( i, i=0..N,0.01 )])):
 > m:=2:
 > A13:=dsolve({eq11,eq21,bc21},numeric,continuation=lambda1,maxmesh=15000, output=array([seq( i, i=0..N,0.01 )])):
 >
 >
 > with(plots):
 > p1:=odeplot(A11, [[eta,f(eta),color=red,linestyle=2]],0..N): p2:=odeplot(A12, [[eta,f(eta),color=blue,linestyle=2]],0..N): p3:=odeplot(A13, [[eta,f(eta),color=green,linestyle=2]],0..N):
 >
 > display({p1,p2,p3}, axes=boxed,  title=velerature Profiles:);
 >
 >
 > for j from 1 to 100*(N-1) do eta[j]:=A11[2,1][j,1] end do:
 > for j from 1 to 100*(N-1) do vel1[j]:=A11[2,1][j,2] end do:
 > for j from 1 to 100*(N-1) do vel2[j]:=A12[2,1][j,2] end do:
 > for j from 1 to 100*(N-1) do vel3[j]:=A13[2,1][j,2] end do:
 >
 > with(LinearAlgebra): E1:= <|||>;
 (5)
 >
 > #ExportMatrix("F:/D/Drive D/Local/Rashad/Paper 9/vel1.txt", E1, format = rectangular):
 >

 > with(plots):
 > p1:=odeplot(A11, [[eta,theta(eta), color=red,linestyle=2]],0..N): p2:=odeplot(A12, [[eta,theta(eta), color=blue,linestyle=2]],0..N): p3:=odeplot(A13, [[eta,theta(eta), color=green,linestyle=2]],0..N):
 >
 > display({p1,p2,p3}, axes=boxed,  title=temperature Profiles:);
 >
 >

 > for j from 1 to 100*(N-1) do eta[j]:=A11[2,1][j,1] end do:
 > for j from 1 to 100*(N-1) do temp1[j]:=A11[2,1][j,5] end do:
 > for j from 1 to 100*(N-1) do temp2[j]:=A12[2,1][j,5] end do:
 > for j from 1 to 100*(N-1) do temp3[j]:=A13[2,1][j,5] end do:
 >
 > with(LinearAlgebra): E2:= <|||>;
 (6)
 > #ExportMatrix("F:/D/Drive D/Local/Rashad/Paper 9/temp1.txt", E2, format = rectangular):
 >
 >

How to numerically solve a system of linear ODEs w...

Dear Community,

I have a four compartment flow model described with the following system of linear ODEs:

diff(L[1](t), t) = L[2](t)/T[21] - L[1](t)/T[12]

diff(L[2](t), t) = L[1](t)/T[12] - L[2](t)/T[21] + L[3](t)/T[32] - L[2](t)/T[23] + Q(t)

diff(L[3](t), t) = L[2](t)/T[23] - L[3](t)/T[32] + L[4](t)/T[43] - L[3](t)/T[34]

diff(L[4](t), t) = L[3](t)/T[34] - L[4](t)/T[43]

(For easier readability I’ve also described the problem in the attached FourCompartmentModelDescription.docx Word file.)

The time constants are as follows:

T12 = 23.1481 d

T21 = 5.4537 d

T23 =  9.752 d

T32 = 14.9007 d

T34 = 8.8235 d

T43 = 363.7255 d

Initial conditions are as follows:

L1(0) = 2500.

L2(0) = 589.

L3(0) = 900.

L4(0) = 37100.

Simulation should run from 0 to 400 d.

Could you pls. help me, to solve it numerically in Maple? As a solution I would need L1(t), L2(t), L3(t) and L4(t) both numerically and graphically. I used Maple to solve ODEs previously, but this time I don’t know, how to correctly specify the problem with an external, tabulated source. I’ve also prepared and attached a Maple worksheet, but did not try to run it yet. (FourCompartments.mw). The external Q(t) source for compartment 2 is attached as Source.xlsx. Should I've forgotten something, pls. let me know.

best regards

Andras

MaplePrimesFourCompartmentModelDescription.docx

FourCompartments.mw

Source.xlsx

 First 11 12 13 14 15 16 17 Last Page 13 of 45
﻿