null hypothesis test : statistic ...

Dear all;

I need a help to get a simple code about the null hypothesis test.

A drug is administrad to a population X of size 50 while a  placebo is given to a population Y of size 25.

Observed results of good bad and no effects are given in the following vectors for both population.

X=[ 20,11,19];

Y=[4,4,17];

test the null hypothesis H0: population independent of treatment versus the one tailed alternative that they are dependent by computing the theoretical contingency table with entries T[i,j] where i=1,2

for the two rows and j=1,2,3 for the three columns. At what p-value can we reject H0.

Many thanks

Curvature of DNA

Maple 18

With this application, you visualize the DNA chain using the position vector as a fundamental tool to describe its curvature and radius of curvature. I also show the native maple syntax for the graphics. You can download the maple center app to show different DNA trajectories based on the position vector. Developed for students of health sciences.

Lenin Araujo Castillo

how to write f''(0) in BC...

in this problem i used slip coundary condition.

the plot want to starts in 0.2, but the plot starts in 0

bc := f(0) = 0, (D(f))(0) = ak*((D^2)(f))(0), (D(f))(N) = 1, g(0) = -de*((D^2)(f))(0), g(N) = 0, T(0) = 1+be*(D(T))(0), T(N) = 0:

How to type double drivative in BC.

how to find X value using psi function and streaml...

ψ =-0.09,-0.07,-0.04,-0.01,0,0.01,0.04,0.07,0.09

these are the ψ values.

then X=

here we can take eta =0..2 and X=-15..15
using this relation how to plot streamlines for eta against X.

Code:
restart; with(plots); fcns := {T(eta), f(eta)}; ep := .1; M := 1; kp := .5; n := 1; ec := .1; pr := 1; s := .1; N := 5; sys := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))*(diff(f(eta), eta))+ep*ep+(M+1/kp)*(ep-(diff(f(eta), eta))) = 0, diff(T(eta), eta, eta)+pr*(f(eta)*(diff(T(eta), eta))-n*(diff(f(eta), eta))*T(eta))+pr*(ec*(diff(f(eta), eta, eta))*(diff(f(eta), eta, eta))+ec*(M+1/kp)*(diff(f(eta), eta))^2+s*T(eta)) = 0; bc := f(0) = 0, (D(f))(0) = 1, (D(f))(N) = ep, T(0) = 1, T(N) = 0; R := dsolve(eval({bc, sys}), numeric, method = bvp[midrich], abserr = 0.1e-9, output = operator); psi = [-0.9e-1, -0.7e-1, -0.4e-1, -0.1e-1, 0, 0.1e-1, 0.4e-1, 0.7e-1, 0.9e-1]; for i to 9 do X[i] = psi[i]/f(eta); print(plots:-contourplot(X[i](X, eta), eta = 0 .. N, X = 0 .. 6)) end do

How do I resolve this program whit statement "for ...

make a program that generates 20 numbers between 1 and 100, calculate the sum and the average of even numbers

please help, I do not know how to do it and the teacher wants this with "for, do and an external accountant

Regression with Maple

Maple 18

Insertion and processing of data (independent and dependent) with Maple syntax and using traditional equations. In the same way we reach the same result. We can also calculate a, b, Sa, Sb and r (pearson correlation coefficient). It can be used for students and researchers.

Regression_with_Maple.mw

Lenin Araujo Castillo

Force and Moment for Health Sciences

Maple 18

Hello student friends in this video we trained in vector operations on the plane and in the space using the native Maple syntax as if you were working on your notebook. Then I explain a model biomechanics exercise applying the vectors learned in the training, developed exclusively for students of health sciences.

Force_and_Moment_for_Health_Sciences.mw

Lenin Araujo Castillo

How to write an equation without compiling?...

I computed A and B matrices, now I want to write it in state space reperentation.

where

So I calculated and got this:

where but the formation is not anymore. How can I enforce Maple to output the result in the form of ?

Force and projection 2D & 3D

Maple 18

Set of vector mechanics exercises in the plane and space using the result technique in line (combining the key ALT + ENTER) also the unit package using the law of the triangle. It is observed that the solution is totally optimal. I leave your constructive criticism to the community's criteria. I hope that someone will raise an alternative solution using the minimum number of lines but that the students will learn. In spanish.

Ejercicios_de_Vectores_Fuerza_y_Proyecciones_2D_y_3D.mw

Lenin Araujo Castillo

Equilibrium with "Shortcut"

Maple 18

Important use of the embedded components called Shortcut applied to vector mechanics exercises for engineering students. This makes each solution of each problem open independently and thus this way to explain in class. To use this worksheet, first unzip all the files in a single folder. In spanish

Equilibrium_with_Shortcut.zip

Lenin Araujo Castillo

How do I solve first ODE using Langrange and Newto...

First, I want to say thank you to all who contributed to previous questions. God bless you.

I need Maple code in solving first order differential equation using Langrange and Newton's Interpolation Method.

The aim is to compare these two numerical result with the exact in tabular form and also to plot the graph.

Some questions are attached.

ODE_Questions.pdf

Distance between two points and midpoint

Maple 18

In this app you can visualize the location of the points in the different quadrants, also calculate the distance between two points. Finally the calculation of the coordinates of the midpoint. With these applications can be combined to study different cases between distance between two points and midpoint. Generated in Maple for students of secondary education and pre-calculation. In Spanish

Distance_between_two_points_and_midpoint.mw

Lenin Araujo Castillo

Increase speed of execution of the dsolve function...

Is it possible to increase the execution speed of the dsolve command? For example, is it possible to increase speed with the use of Multithreaded tools?

How to plot Contour...

with(plots); R1 := .1; R0 := .1; m := .1; a := .1; Ha := .1; Nt := .1; Nb := .1; Pr := 6.2; Le := .6; Bi := 1; Ec := .1; k := 1; r := .1; A := 1; fcns := {C(y), T(y), U(y), W(y)}; sys := diff(U(y), $(y, 2))+(R1*(diff(U(y), y))-2*R0*W(y))*exp(a*T(y))-a*(diff(U(y), y))*(diff(T(y), y))-Ha*(U(y)+m*W(y))*exp(a*T(y))/(m^2+1)-U(y)/k+A*exp(a*T(y)) = 0, diff(W(y), $(y, 2))+(R1*(diff(W(y), y))+2*R0*U(y))*exp(a*T(y))-a*(diff(W(y), y))*(diff(T(y), y))-Ha*(W(y)-m*U(y))*exp(a*T(y))/(m^2+1)-W(y)/k = 0, diff(T(y), $(y, 2))+R1*Pr*(diff(T(y), y))+Pr*Ec*exp(-a*T(y))*((diff(U(y), y))*(diff(U(y), y))+(diff(W(y), y))*(diff(W(y), y)))+Nt*(diff(T(y), y))*(diff(T(y), y))+Pr*Ec*(U(y)*U(y)+W(y)*W(y))*exp(-a*T(y))/k = 0, diff(C(y), $(y, 2))+Pr*Le*R1*(diff(C(y), y))+Nt*(diff(C(y), \$(y, 2)))/Nb-r*C(y) = 0; bc := U(0) = 0, W(0) = 0, C(0) = 0, (D(T))(0) = Bi*(T(0)-1), U(1) = 0, W(1) = 0, C(1) = 1, T(1) = 0; L := [.5, 1.0, 1.5, 2.0]; AP := NULL; for k to 4 do R := dsolve(eval({bc, sys}, Ha = L[k]), fcns, type = numeric, method = bvp[midrich], AP); AP := approxsoln = R; p1u[k] := odeplot(R, [y, U(y)], 0 .. 1, numpoints = 100, labels = ["y", "U"], linestyle = dash, color = black) end do; display({p1u[1], p1u[2], p1u[3], p1u[4]})

numeric integration fails?...

int((1-x^floor(u))/((1-x)*u^2), u = 1 .. infinity, numeric)