Rewrite the code that counts the number of primes less than using an if-then statement.  Implement    your code where j goes from 2 to 15.

j:=1;
                               1
for i from 2 to 10
while ithprime(j)<2^i do
j:=j+1
end:
print(2^i,primes=j-1):
                        2048, primes = 9

I need to edit this code to satisfy a IF then Statement. can any one help me out?

\

regards "Geordi"

In doing multiple linear regression with Statistics:-LinearFit, how do I get or compute the values commonly called R^2 and R^2[adjusted], also known as the coefficient of determination? I know that residualsumofsquares is part of it. I also need the "total sum of squares" to compute R^2. And how do I modify that to get R^2[adjusted]? These things do not seem to be among the numerous output options to LinearFit. These values (R^2 and R^2[adjusted]) are typically part of the output of other statistics software when doing multiple regression.

I'd also like the p - values for the significance of the individual parameters and the p - value for the global utility---also things that are standardly reported by statistics software.

I am trying to create a procedure that can solve integrals using the Composite Simpson's 3/8 rule. However when I test my procedure against maple's ApproximateInt I am getting the wrong results.

Here is my attempt:

restart;

f:= x -> exp(x)*sin(4*x); # function I am using

simp := proc(a, b, n)
  local h, sum, i, single:
  h := (b-a)/n:
  sum := 0:
  single := (3*h/8) * (f(a) + f(b)): # this is the end points
    for i from a+h by h to b-h do
       sum := sum + (3*h/8) * (3*f(i)):
    end do:
print(evalf(sum + single));
end proc:

simp(0,1,12);
                                                                                0.6224486445
evalf(Student:-Calculus1:-ApproximateInt(f(x), 0..1, method = simpson[3/8], partition=12));

                                                                                0.5323516717

As you can see my answer is not very close to the answer given by Maple. I am not sure why my procedure simp is wrong.

If i use expression with function add() it runs normaly, but if I use expression with Threads:-Add (parallel implementation) it causes error "Error, continuation task already created for the current task" or "Kernel connection has been lost"

Expression:

Array(1 .. N,1 .. 1/2*N-NR-1,(i, m) -> evalf(Add(cosArr2[modp(k*i,N)]/kl[k],k = 1 .. NR+m-1)+Add(cosArr2[modp(k*i,N)]*alpha[k],k = NR+m .. 1/2*N)))

What am I doing wrong? Can I use two Add in one evalf?

How can I create a new object from a stanrd, defined in a given metric? For example I want to define an object by convolution of two indices of the Christoffel symbols? I cannot even to see components of such object.

Many thans for an explanation

Leonid

Hi, there i want to basically plot the graph attached. It is basically the graph y= 2x^2 transformed to y = x^2. I can plot them, but how do they get the arrows showing the transformation and the arrows on the axis. Any help would be appreciated.

Thanks

Composite1:=proc(g,f)
Flist := indets(f): Glist := indets(g): gof:=g:
Subslist1 := [seq( Flist[i]=q, i=1.. nops(Flist))]:
return subs(q=f,subs(Subslist1, gof)): # g(f)
end proc:
#F(g o f) = F(g) o F(f) = F o g o f = (F o g) o (F o f)
Functor := (1/2)*(-y*t2-x*t1-y*t3+sqrt(y^2*t2^2+2*y*t2*x*t1+2*y^2*t2*t3+x^2*t1^2+2*x*t1*y*t3+y^2*t3^2-4*x*t4*y*t9-4*x*t4*y*t8-4*x^2*t4*t7-4*y^2*t5*t9-4*y^2*t5*t8-4*y*t5*x*t7-4*y^2*t6*t9-4*y^2*t6*t8-4*y*t6*x*t7))/(x*t4+y*t5+y*t6);
F1:=x+2;
G1:=3*x+5;
gof:=subs(x=F1, G1);
osys := Composite1(gof, Functor) = Composite1(Composite1(G1, Functor),Composite1(F1, Functor));
sys1 := subs(x=3,subs(y=2, osys));sys2 := subs(x=5,subs(y=1, osys));
sys3 := subs(x=1,subs(y=5, osys));sys4 := subs(x=1,subs(y=2, osys));
sys5 := subs(x=2,subs(y=5, osys));sys6 := subs(x=5,subs(y=2, osys));
sys7 := subs(x=2,subs(y=1, osys));sys8 := subs(x=3,subs(y=5, osys));
sys9 := subs(x=5,subs(y=3, osys));
res:=solve([sys1, sys2, sys3, sys4, sys5, sys6, sys7, sys8, sys9], {t1,t2,t3,t4,t5,t6,t7,t8,t9});
simplify(%);
`~`[lhs](select(evalb, res));

alpha:= (1/2)*(-y*t2-x*t1-y*t3+sqrt(y^2*t2^2+2*y*t2*x*t1+2*y^2*t2*t3+x^2*t1^2+2*x*t1*y*t3+y^2*t3^2-4*x*t4*y*t9-4*x*t4*y*t8-4*x^2*t4*t7-4*y^2*t5*t9-4*y^2*t5*t8-4*y*t5*x*t7-4*y^2*t6*t9-4*y^2*t6*t8-4*y*t6*x*t7))/(x*t4+y*t5+y*t6);
g := -y/x;
f := (-x+sqrt(x^2-x*y-2*y^2))/(2*y+x);
subs(p=f,subs(q=f,subs(x=p,subs(y=q,g))));
g := (1/2)*(-x+sqrt(x^2-4*y*x-4*y^2))/(x+y);
f := x*y;
gof := subs(p=f,subs(q=f,subs(x=p,subs(y=q,g))));
lhsgofoalpha := subs(q= alpha,subs(p=alpha, subs(x=p,subs(y=q,gof))));
foalpha := subs(p= alpha,subs(q=alpha,subs(x=p,subs(y=q,f))));
rhsgofoalpha := subs(x= foalpha,subs(y= foalpha, g));
osys := lhsgofoalpha = rhsgofoalpha;
sys1 := subs(x=3,subs(y=2, osys));
sys2 := subs(x=5,subs(y=1, osys));
sys3 := subs(x=1,subs(y=5, osys));
sys4 := subs(x=1,subs(y=2, osys));
sys5 := subs(x=2,subs(y=5, osys));
sys6 := subs(x=5,subs(y=2, osys));
sys7 := subs(x=2,subs(y=1, osys));
sys8 := subs(x=3,subs(y=5, osys));
sys9 := subs(x=5,subs(y=3, osys));
res:=solve([Re(sys1)=0, Re(sys2) =0, Re(sys3) =0, Re(sys4) =0, Re(sys5) =0, Re(sys6) =0, Re(sys7) =0, Re(sys8) =0, Re(sys9) =0], {t1,t2,t3,t4,t5,t6,t7,t8,t9});

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");

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? 

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

i assign 

seta := [x+1, x^2]

setb := [x^3, 2*x+5]

does morphism mean that

i use card_prod

to get

(x+1, x^3)

(x+1, 2*x+5)

(x^2, x^3)

(x^2, 2*x+5)

such that i composite each of 4 sets still satisfy F(f o g) = f o g

example

subs(x=x^3, x+1)

(A o C) o Colimit = (B o C) o Colimit

if known A, B, C and framework of Colimit

can it be said colimit?

Bonus, what are A,B,C? <- this can be not answered

after integration,
result is x^2

f := x^2

we know in the range 5..-5 = 5^2 - (-5)^2 = 0

this example is not good enough, as i can be any number x, -x

if in another example, i use a and b to represent that would like to find

3*a^2 - 3*b^2 = 0, how to find possible a and b?

my example is not good, is there example that a and b are different number, not only in sign

Notepad displays alt+1, alt+2 and alt+3 as ☺☻♥

Just wondering why Maple displays them as square boxes?  Shouldn't it have the font capability to display them as well?