MaplePrimes Questions

I got decimal places problem,

not known correct or not

the value are different but difference is constant 

i do not know how many places needed to get exact result

i do not believe the difference is constant

because the matrix are different

but even if using 36 decimal places still constant, 

i notice increasing decimal places , the constant difference is changed

is it possible  to output fraction when calculate eigenvector?

Partial_for_R0.mw
Anybody knows what happen to this coding?

restart

unprotect(Pi)

k := .5;

.5

 

.95

 

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-1

 

.2

 

1.2

(1)

R := k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))

65688.31239

(2)

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), k)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ε`)

Error, invalid input: diff received 0.2e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), tau)

Error, invalid input: diff received .95, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), phi)

Error, invalid input: diff received 0.25e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), chi)

Error, invalid input: diff received 0.598e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), p)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ϑ`)

Error, invalid input: diff received 0.8e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), mu)

Error, invalid input: diff received 0.1e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), alpha)

Error, invalid input: diff received 0.57e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), rho)

Error, invalid input: diff received 0.5e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), beta)

Error, invalid input: diff received 0.115e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), eta)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), q)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

NULL

NULL


 

Download Partial_for_R0.mw
 

restart

unprotect(Pi)

k := .5;

.5

 

.95

 

0.1e-1

 

116.1

 

0.8e-2

 

0.25e-2

 

0.2e-2

 

0.5e-1

 

0.115e-1

 

0.598e-2

 

.5

 

.2

 

.1

 

0.57e-1

 

.2

 

1.2

(1)

R := k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))

65688.31239

(2)

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), k)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ε`)

Error, invalid input: diff received 0.2e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), tau)

Error, invalid input: diff received .95, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), phi)

Error, invalid input: diff received 0.25e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), chi)

Error, invalid input: diff received 0.598e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), p)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), `ϑ`)

Error, invalid input: diff received 0.8e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), mu)

Error, invalid input: diff received 0.1e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), alpha)

Error, invalid input: diff received 0.57e-2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), rho)

Error, invalid input: diff received 0.5e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), beta)

Error, invalid input: diff received 0.115e-1, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), eta)

Error, invalid input: diff received .2, which is not valid for its 2nd argument

 

diff(k*tau*(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))*(Pi*(-mu*p+mu+phi)/(mu*(`ϑ`+mu+phi))+`ε`*Pi*(mu*p+`ϑ`)/(mu*(`ϑ`+mu+phi)))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q)), q)

Error, invalid input: diff received .5, which is not valid for its 2nd argument

 

NULL

NULL


 

Download Partial_for_R0.mw

 

 

Basic_Reproduction_Number.mw
Hello, anybody knows how to compute sensitivity analysis for basic reproduction number using Maple coding? By taking the partial derivatives of all the parameters involved in the basic reproduction number using normalised forward sensitivity index formula. Thanks.

restart

unprotect(Pi)

R[0] := k*tau*[(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))]*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))

k*tau*[(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))]*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))

(1)

``


 

Download Basic_Reproduction_Number.mw

 

I am having difficulties with a recursive function call with named parameters.

The proc is defind as follows:

Subs:=proc(eqn::seq(equation),elemt::Element,{num::boolean:=false},$) option remember;
...
qs[i]=Subs(eqn,elemt[qs[i]]); # This is where I need to add something like 'num'=num
...
end proc;

My difficulty is with the "num" flag. In the code this flag governs whether to try to evalf() certain results or not. That and the recursion per se all works. What does not work is when I add the num option to the internal call. I have been trying 'num'=num, ''num''=num, "num"=num and `num`=num.  I either get an error that something like true=true is not a valid argument here, or that "num" is not valid.

Any hint?

Thanks,

Mac Dude

Let G be a permutation group like S3. I have a loop in which x runs over the elements of this group.

When I ask the cycle type of x or the number of orbits of x, Maple gives an incorrect result in some cases. This is because when x=(12), for example, Maple considers x to be an element of the permutation group S2 instead of S3, i.e. it fails to understand that (12) should really be (12)(3) in this case.

Is there a way I can instruct Maple to consider the permutation (12) as a member of S3 instead of S2? 

I am posting the code as requested. Something like this:

with(GroupTheory);

for x in Elements(S3) do print(x, numelems(Orbits(PermutationGroup(x)))) end do;

(I don't know how to format this, sorry)

This code is simply asking how many orbits exist in the action of each element of S3. The result will be 1 for x=(12), however, because Maple will assume that (12) is acting on the set {1,2}, even though I am considering (12) as an element of S3. I want it to consider (12)=(12)(3) and let it act on the set {1,2,3}, thereby producing 2 orbits, not 1.

Ok embarassing I posted this one and now it outputs this java error file log and yep i dont know that much about java at all. 

 

Sometimes it works, sometimes it closes maple automatically and spits out the .txt error log i have attached if thats any help... MAPLE_EXAMPLE_16.mw  i honestly have never ever seen maple behave this way in more than a decade of playing around in it with stuff like this.

EDIT: Here is a version with smaller plot components that seems to be working.... MAPLE_EXAMPLE_17.mw

It says i am not allowed to upload a log file but this is what it looks like:

#
# A fatal error has been detected by the Java Runtime Environment:
#
# java.lang.OutOfMemoryError: requested 1024000 bytes for GrET in C:\BUILD_AREA\jdk6_18\hotspot\src\share\vm\utilities\growableArray.cpp. Out of swap space?
#
 
V  [jvm.dll+0x15df8a]
V  [jvm.dll+0x1e1e14]
V  [jvm.dll+0x1a1aad]
V  [jvm.dll+0xc834f]
V  [jvm.dll+0xca01c]
V  [jvm.dll+0xca370]
V  [jvm.dll+0xce42a]
V  [jvm.dll+0x1d8592]
V  [jvm.dll+0xc9398]
V  [jvm.dll
Java Threads: ( => current thread )
  0x522e7c00 JavaThread "Timer-28" [_thread_blocked, id=23236, stack(0x5f2f0000,0x5f6f0000)]
  0x522e7400 JavaThread "Timer-22" [_thread_blocked, id=24248, stack(0x6b440000,0x6b840000)]
  0x522e2c00 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=11140, stack(0x6a2d0000,0x6a6d0000)]
  0x522e7000 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=6792, stack(0x64d30000,0x65130000)]
  0x522e6400 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=20280, stack(0x6a6d0000,0x6aad0000)]
  0x522e4400 JavaThread "Timer-18" daemon [_thread_blocked, id=23888, stack(0x63f20000,0x64320000)]
  0x522e6800 JavaThread "Timer-17" [_thread_blocked, id=17320, stack(0x69ed0000,0x6a2d0000)]
  0x522e5000 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=13480, stack(0x69ad0000,0x69ed0000)]
  0x522e4800 JavaThread "Timer-16" daemon [_thread_blocked, id=18184, stack(0x68ed0000,0x692d0000)]
  0x522e3000 JavaThread "Timer-15" [_thread_blocked, id=23296, stack(0x692d0000,0x696d0000)]
  0x522e3800 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=16388, stack(0x66820000,0x66c20000)]
  0x522e1400 JavaThread "Timer-14" daemon [_thread_blocked, id=24500, stack(0x61c50000,0x62050000)]
  0x522e0c00 JavaThread "Timer-13" [_thread_blocked, id=15848, stack(0x659b0000,0x65db0000)]
  0x522e0000 JavaThread "WMI:MapleClientSocket:Kernel Connection " daemon [_thread_blocked, id=7260, stack(0x655b0000,0x659b0000)]
  0x466cb000 JavaThread "Timer-12" daemon [_thread_blocked, id=1112, stack(0x64320000,0x64720000)]
  0x4699b800 JavaThread "Timer-11" [_thread_blocked, id=24540, stack(0x64930000,0x64d30000)]
  0x4699b400 JavaThread "Timer-10" daemon [_thread_blocked, id=21144, stack(0x60bf0000,0x60ff0000)]
  0x4699d400 JavaTh
 

Hi. I have an excel data file.

IDEALLY, I want to import the results of the formula in cells C3:C7.

But it doesn't like it. Reports 0's

When I convert the range to result to values in cells D3:D7 it imports no problem. But I don't want to do this step.

any ideas?

Book1.xls

import.mw

 

 

 

Lett_L_rot_AL_MapleQuestion.mws

If the attached program is run, at the end are two animation outputs.  One is of two letters A; the other of the single second A being animated.   The second letter A animation works fine, but in the former it "blips" on Continuous animation, and on 'single cycle' the A totally disappears.  I'm using the same number of frames, 48, for each. 
   Prior to that is the letter L which is meant to give the illusion of it rotating, and the letters A getting out of the way.

   Any help on ridding the "blip" most appreciated. 

Hello, 
I have this code to solve partial equation, and I need to writdata in file .dat or.txt using table and loops or if there is another way 
Any help please,
Thank you 

 

``

restart

``

eq := diff(Tf(x, y), x) = diff(Tf(x, y), y, y):

Cl := {1 = Tf(0, y), (D[2](Tf))(x, 0) = 1, (D[2](Tf))(x, 1) = 1};

{1 = Tf(0, y), (D[2](Tf))(x, 0) = 1, (D[2](Tf))(x, 1) = 1}

(1)

``

pds := pdsolve(eq, Cl, numeric);

_m221385792

(2)

``


 

Download mapleprime.mw

is there a command line method for exporting plots as jpg files and gif files?

How do I create a wave in Maple such that the end which the wave hits of (fixed boundary condition) is formed of triangles rather than just a straight line. I have, 

with(plots):

animate(plot3d,[sin(Pi*(t-x))sin(Piy)Heaviside(t-x)-sin(Pi(t+x))sin(Piy)*Heaviside(t+x),x=0..4,y=0..1],t=-2*Pi..2*Pi);

Would I have to create a whole new formula or can I just add something to the code to give this different boundary? Or maybe I can use a different software other than maple, I just want to know which directions the waves will travel when they hit a ridged surface.

 

Say I wanted to create a five sides polygon, so the wave hits these two boundary conditions so i can see what happens next. (I added the two black lines in paint so it was easier to visualise)

How to translate the elliptic curve from the 4th to the 3rd curve with maple?

thanks.

 

Hi, I am using fsolve to find the roots x of a nonpolynomial expression f(x,a) = 0, using a loop for do to vary the values of the parameter a in small identical steps. At some stage, fsolve cannot find a solution. According to the Maple help, in this case fsolve returns a NULL sequence. I would like to know how to use the fact that fsolve cannot find a solution to break the loop with an if-then statement. 

How can plot  f := proc (x, y) options operator, arrow; sin(x)*cos(y) end proc

as surface and contour (contour is project on xy plan)

HOW I can find the function that satisfies these boundary conditions?

please help me.

Thanks

FUNCTION.mw

Download FUNCTION.mw

 

First 858 859 860 861 862 863 864 Last Page 860 of 2427