litun

85 Reputation

7 Badges

7 years, 12 days

MaplePrimes Activity


These are replies submitted by litun

@Carl Love I am getting that sum by manual addition. Can you please rewrite my program correctly to get that sum. Also can we give a restriction to the number of roots satisfying the condition that component sum is a constant.

restart; with(LinearAlgebra); with(inttrans); with(RootFinding); with(PolynomialTools); with(powseries);
Digits := 10;
10
Bstar := 3.60/(12.0+s)+14.00/(20.0+s);
3.60 14.00
-------- + --------
12.0 + s 20.0 + s
Nrot := [-23.8279785337007-4.37234259951331*I, -23.8279785337007+4.37234259951331*I, -13.8480664440518, -12.2635728941426-2.65264168833545*I, -12.2635728941426+2.65264168833545*I, -4.24998018723859];
[-23.8279785337007 - 4.37234259951331 I,

-23.8279785337007 + 4.37234259951331 I, -13.8480664440518,

-12.2635728941426 - 2.65264168833545 I,

-12.2635728941426 + 2.65264168833545 I, -4.24998018723859]
rfw := unapply(product(-Bstar*gm[i]+1, i = 1 .. 6), s);
for i to 6 do b[i] := rfw(Nrot[i]) end do;
solve({b[1] = 0, b[2] = 0, b[3] = 0, b[4] = 0, b[5] = 0, b[6] = 0}, {gm[1], gm[2], gm[3], gm[4], gm[5], gm[6]});
Warning, returning only the first 100 solutions, increase _MaxSols to see more solutions

questn.mw       

please see the attachment above for a specific example.

 

why my upload file is not visible

 

I want to get solutions in terms of a and c.

 

When I plot the function in Maple I didn't get the actual curve.

i.e at point t=0 I am not getting the value f(0) rather I got the value being evaluated from the function f(t) putting t=0

In my case f(0) is not equal to f(t) at t=0.

When I plot the function in Maple I didn't get the actual curve.

i.e at point t=0 I am not getting the value f(0) rather I got the value being evaluated from the function f(t) putting t=0

In my case f(0) is not equal to f(t) at t=0.

Dear sir

 I tried with that but I failed.

I got the message "Error, invalid parameters for inline function"

Dear sir

 I tried with that but I failed.

I got the message "Error, invalid parameters for inline function"


 restart; with(LinearAlgebra):with(RootFinding):  Digits:=100;  lambda:=5,mu:=6;  erlang:=m->(lambda^(m)*t^((m-1))*(e)^(-(lambda*t)))/((m-1)!); erlang(3);

f := proc (z) options operator, arrow; lambda^3/(lambda+mu-mu*z)^3 end proc

g := taylor(f(z), z = 0, 100); evalf(g); b := eval(g, z = 1); evalf(b)

h := proc (z) options operator, arrow; z-f(z) end proc; r[0] := fsolve(h(z), z, complex)

q := proc (n) options operator, arrow; (1-r[0])*r[0]^n end proc; q(0), q(1), q(2)

``

p := proc (n) options operator, arrow; lambda*q(n-1)/mu end proc; p(1), p(2), p(3), p(100), p(1000); p(0); p(-1)

NULL


Download Assignment.mw

 

Thank you

 

Thank you

 

> g := taylor(f(z), z = 0, 100);
 evalf(g);
2 3
0.02700000000 + 0.05670000000 z + 0.07938000000 z + 0.09261000000 z

4 5 6 7
+ 0.09724050000 z + 0.09529569000 z + 0.08894264400 z + 0.08004837960 z

8 9 10
+ 0.07004233215 z + 0.05992510640 z + 0.05033708937 z

11 12 13
+ 0.04164250121 z + 0.03400804265 z + 0.02746803445 z

14 15 16
+ 0.02197442756 z + 0.01743304586 z + 0.01372852362 z

17 18 19
+ 0.01074055083 z + 0.008353761757 z + 0.006463173570 z

20 21 22
+ 0.004976643649 z + 0.003815426798 z + 0.002913598645 z

23 24 25
+ 0.002216868535 z + 0.001681125305 z + 0.001270930731 z

26 27 28
+ 0.0009580862433 z + 0.0007203389162 z + 0.0005402541872 z

29 30 31
+ 0.0004042591676 z + 0.0003018468452 z + 0.0002249245846 z

32 33 34
+ 0.0001672876598 z + 0.0001241984141 z + 0.00009205294222 z

35 36 37
+ 0.00006811917724 z + 0.00005033250319 z + 0.00003713722532 z

38 39 40
+ 0.00002736427129 z + 0.00002013729708 z + 0.00001480091335 z

41 42 43
+ 0.00001086603639 z + 0.000007968426685 z + 0.000005837335827 z

44 45 46
+ 0.000004271868492 z + 0.000003123210519 z + 0.000002281301597 z

47 48 -7 49
+ 0.000001664864782 z + 0.000001213963904 z + 8.844594156 10 z

-7 50 -7 51 -7 52
+ 6.438864546 10 z + 4.683958326 10 z + 3.404877399 10 z

-7 53 -7 54 -7 55
+ 2.473354337 10 z + 1.795472037 10 z + 1.302533351 10 z

-8 56 -8 57 -8 58
+ 9.443366792 10 z + 6.842299096 10 z + 4.954768311 10 z

-8 59 -8 60 -8 61
+ 3.585908591 10 z + 2.593807214 10 z + 1.875195052 10 z

-8 62 -9 63 -9 64
+ 1.354979650 10 z + 9.785964141 10 z + 7.064242864 10 z

-9 65 -9 66 -9 67
+ 5.097122928 10 z + 3.676106839 10 z + 2.650088960 10 z

-9 68 -9 69 -10 70
+ 1.909622927 10 z + 1.375482021 10 z + 9.903470554 10 z

-10 71 -10 72 -10 73
+ 7.127709089 10 z + 5.127990706 10 z + 3.687938521 10 z

-10 74 -10 75 -10 76
+ 2.651328775 10 z + 1.905421613 10 z + 1.368895001 10 z

-11 77 -11 78 -11 79
+ 9.831155005 10 z + 7.058265132 10 z + 5.065868772 10 z

-11 80 -11 81 -11 82
+ 3.634760844 10 z + 2.607155618 10 z + 1.869521345 10 z

-11 83 -12 84 -12 85
+ 1.340199037 10 z + 9.604759763 10 z + 6.881527877 10 z

-12 86 -12 87 -12 88
+ 4.929094387 10 z + 3.529684831 10 z + 2.526933458 10 z

-12 89 -12 90 -13 91
+ 1.808602936 10 z + 1.294155879 10 z + 9.258192055 10 z

-13 92 -13 93 -13 94
+ 6.621619970 10 z + 4.734814279 10 z + 3.384888506 10 z

-13 95 -13 96 -13 97
+ 2.419304522 10 z + 1.728794690 10 z + 1.235107959 10 z

-14 98 -14 99 / 100\
+ 8.822199704 10 z + 6.300298173 10 z + O\z /
 
how to get the coefficients as a(k)s for different values of k

@Markiyan Hirnyk 

Thank you sir for your kind help.

the problem was with the imaginary unit "I" which by wrong i ve entered "i".

@Markiyan Hirnyk 

Dear sir i have worked with the command

>plot(log(-1+exp(t)),t=0..2*Pi);

i got plot.

When i replced the argument by complex variable t ->i*θ

i got the same message.

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

and got a blank plot.

 

1 2 3 Page 3 of 3