I want to taylor-expand cos(x) to the powers of two, four, and six, then to make equations which does not contain O(x^2) etc. and to obtain its value at x=0.1.

For that, I wrote as the following.
(1):
a:=seq(taylor(cos(x),x=0,i),i=[2,4,6]);
convert(a),polynom);

(2):
a:=seq(taylor(cos(x),x=0,i),i=[2,4,6]);
seq(convert(op(i,a),polynom),i=[2,4,6]);

(3):
a:=seq(taylor(cos(x),x=0,i),i=[2,4,6]):
convert(a,polynom);

All of the above were not good.

On earth, I couldn't pick up each of the result of
a:=seq(taylor(cos(x),x=0,i),i=[2,4,6])
with op(1,a).

How can I solve this problem?

I have created a diagram with "plots" and "plottools" that involves the line command creating a line in 3-space.  It works just fine on all of my machines and versions of Maple but it fails on Maple 17 in the Windows 8 environment.   By "fai;s" I mean that it does not draw the line; it does not generate an error message.   Any idea what's happening? 

I wrote the following code inside a code region in order to implement a scalar product for lists:

m:= module() option package; export `*`;
    `*`:= proc( a::list, b::list) option overload;
        return add(i, i  in zip((x, y)-> x*y, a, b));
    end proc;    
end module;
l1:=[1,2,3,4];
l2:=[4,1,0,5];
with(m);
l1*l2;

when i run the code I sometimes get the desired result ->26

but sometimes I get

I don't have any clue, where this weird behaviour comes from.

Any help would be greatly appreciated.

Hi! I can't answer the following two simple questions.

1) Why do use different font sizes every other  files created by the attached worksheet?

2) What happens when we uncomment some (all) of the #restart lines, and execute the commands (execution groups)  in the original order? How many files should we get?

Update: sorry, I can't see the uploaded worksheet, but here it is as plain text:

#### on Fedora 18 Linux, Maple 17.01:

#restart;

currentdir(kernelopts(homedir)):
currentdir();

#restart;

plotsetup(ps,plotoutput= `sinplot1`, plotoptions = `width=500,height=500, font=[TIMES, BOLD, 16]`);
plot(sin(x), x=-Pi..Pi, title="Output: sinplot1.eps");

#restart;

plotsetup(ps,plotoutput= `sinplot2`, plotoptions = `width=500,height=500`);
plot(sin(x), x=-Pi..Pi, font=[TIMES, BOLD, 16], title="Output: sinplot2.eps");

#### Restart in the same execution group as the plot statements - no output file?

#restart;
plotsetup(ps,plotoutput= `sinplot3`, plotoptions = `width=500,height=500, font=[TIMES, BOLD, 16]`);
plot(sin(x), x=-Pi..Pi, title="Output: sinplot3.eps");

#restart;
plotsetup(ps,plotoutput= `sinplot4`, plotoptions = `width=500,height=500`);
plot(sin(x), x=-Pi..Pi, font=[TIMES, BOLD, 16], title="Output: sinplot4.eps");

hello, how to write this system?

> restart;
> x10:=0.2; x20:=0.4; U0:=0.25;
> eq3 :=subs(U(t)=(piecewise(U(t)<=0,0,U(t)>=1,1,U(t))),(1/4)*(U(t)^2*x1(t)-2*U(t)*x1(t)*sqrt(x2(t))-U(t)^2*x2(t)+2*U(t)*x2(t)*sqrt(x1(t))-2*(diff(U(t), t))*x1(t)^(3/2)*U(t)+4*(diff(U(t), t))*x1(t)^(3/2)*sqrt(x2(t))+2*(diff(U(t), t))*x2(t)^(3/2)*U(t)-4*(diff(U(t), t))*x2(t)^(3/2)*sqrt(x1(t)))/(x1(t)^(3/2)*x2(t)^(3/2)));
> sys := eq3=0, diff(x1(t), t) = (piecewise(U(t)<=0,0,U(t)>=1,1,U(t)))-sqrt(x1(t)), diff(x2(t), t) = (piecewise(U(t)<=0,0,U(t)>=1,1,U(t)))-sqrt(x2(t));
> dsolve({sys,x1(0)=x10,x2(0)=x20,U(0)=U0},{x1(t),x2(t),U(t)},type=numeric);

Male produces an error: Error, (in DEtools/convertsys) invalid input: eval expects its 2nd argument, eqns, to be of type {integer, equation, set(equation)}, but received newdeqns

Hi,

restart:with(plots):

L1 := (1/2)*(S+sqrt(S^2+4*a*b))/b;

L2 := (1/2)*(S-sqrt(S^2+4*a*b))/b;

solve(L1=L2,S);

         2*sqrt(-a*b), -2*sqrt(-a*b)

solve(L1=L2,a);

             -(1/4)*S^2/b

solve(L1=L2,b);

                -(1/4)*S^2/a

Here is my question. 

Now plotting L1, L2 

a:=1:b:=1: # say

plot({L1,L2},S=0..1);

Is there a way to utilize the conditions which we can get through solve and use it while plotting (not 3d plotting) for any given choice of the values of a, b, S? 

Hallo guys,

I am newbie in Mapple, i could not solve the problem yet, here is my little function. I want to estimate beta with the Maximum Likelihood Estimation with Weibull: 

> MaximumLikelihoodEstimate(Weibull(beta, eta), F1); F1 := Vector([1500, 1750, 2250, 4000, 4300, 5000, 7000]);
print(`output redirected...`); # input placeholder
[beta = HFloat(HFloat(undefined)), eta = HFloat(HFloat(undefined))

]

May someone help me why i just get HFloat(undefined)  ?

Kind regards

Can anyone tell me how to solve the equation above using Maple.  I know that there is a solution around x=0.995, y=0.743, but I cannot induce Maple to find it.  Any help or suggestions would be appreciated.

Ratch

solve matrix equation Ax+xB=C ? where X_(n\times n), and  A,B,C \ in M_n(R) we are given. write maple procedure and check by a 3* 3 case.

Hint( USE MATRIX  CRAMER RULE), OR GAUSS-ELIMINATION.

Hi there,

I'm quite new to Maple so please forgive me! I have a system of partial differential equations I'm trying to solve in Maple as such below 

df/dt = f(1-f) - f * h

dg/dt = g(1-g) * Gradient(1-f * gradient(g))

dh/dt = (g - h) + Laplacian(h),

where f,g,h are functions of space and time (i.e. f(x,y,z,t)). I guess my first question is - is this possible in Maple to evaluate? (I'm currently unsure on ICs as I'm figuring it out from the model - it's a model for cancer growth I'm trying to evaluate but have a rough idea of what I'd use).

If it is possible, can you please share how I'd write this? Everytime I've tried I seem to be failing to define anything properly, so your expertise would be greatly appreciated!

Hi,

I have a problem in which i have to integrate Bessels function of first kind of order 1 from 0 to infinity. The answer that i m getting is -1, but when i am solving it numerically then the answer is close to 1. I dont know why it is happening. Please help me out for it.

Regards

Sunit

The value of Pi can be approximated by the infinite series 4-(4/3)+(4/5)-(4/7)+(4/9)-(4/11)+... . Write a procedure that takes as an input the integer n and uses a for-loop to return the decimal point approximation of using the n^th partial sum of the infinite series. Determine the minimum value of n that gives the true value of to 5 decimal places.  

the 99 bottles of beer on the wall for loop is alot easier than this one. can some one walk me through it? i only have a general idea how to implement it.

What is a bit surpiring,  although the manual doesn't mention it,  the following code produces

the same plots in the 2. and 3. case.  And what about the use of FAIL in similar situations? 

for is_continouos in [true, false, FAIL] do
dc := not is_continouos;
plot(1/x, x=-1..1, discont= dc, thickness=2, color = blue, legend = '1/x')
end do;

im having issues with a problem on my homework. we are supposed to find all primes between 1000 through 1015. using the for loop, for in loop and while loop. 

the while loop seemed the easiest.

i:=1000;
1000
while i<=1015 do
i, "Is it a prime?", isprime(i);
i:=i+1;
end;

the for loop and the for in loop have gotten me abit comfused ive tried several times with each and nothing.

and  x have range=(0 to 2),the y have range = (2,3),follow when the x=(3 to 4),so the y=(3 to -3) how to implement the title's function 

still the x and y have variability value and variability number