tomleslie

6809 Reputation

17 Badges

10 years, 222 days

MaplePrimes Activity


These are answers submitted by tomleslie

The attached solves the system for lambda=0.1. It then uses this solution (via approxsoln=..) to obtain a solution for lambsa=0.2


 

restart

b := 25.4

f := 2.2

h := 10

E := 0.70e5

G := 0.26e5

Ec := 420

Gc := 220

nuc := .25

A11 := E*f

D11 := (1/12)*E*f^3

A55 := G*f

dt := fNULL

db := f

c := h

phi_ct := .55; delta0_ct := .12; phi_cb := .55; delta0_cb := .12

eq1_1 := (diff(y1_1(x), x))/(27.1) = 0

eq1_2 := (diff(y1_2(x), x))/(27.1) = 0

eq1_3 := (diff(y1_3(x), x))/(27.1) = y1_2(x)

eq1_4 := (diff(y1_4(x), x))/(27.1) = 2*b*phi_cb*(y1_17(x)+(1/2)*db*y1_19(x)-y1_20(x)-y1_21(x)*c-y1_22(x)*c^2-y1_23(x)*c^3)/delta0_cb^2

eq1_5 := (diff(y1_5(x), x))/(27.1) = b*phi_cb*(y1_18(x)-y1_24(x)-y1_25(x)*c-y1_26(x)*c^2)/delta0_cb^2

eq1_6 := (diff(y1_6(x), x))/(27.1) = y1_5(x)-2*b*db*phi_cb*(y1_17(x)+(1/2)*db*y1_19(x)-y1_20(x)-y1_21(x)*c-y1_22(x)*c^2-y1_23(x)*c^3)/delta0_cb^2*(1/2)

eq1_7 := (diff(y1_7(x), x))/(27.1) = -2*b*phi_cb*(y1_17(x)+(1/2)*db*y1_19(x)-y1_20(x)-y1_21(x)*c-y1_22(x)*c^2-y1_23(x)*c^3)/delta0_cb^2

eq1_8 := (diff(y1_8(x), x))/(27.1) = -2*b*c*phi_cb*(y1_17(x)+(1/2)*db*y1_19(x)-y1_20(x)-y1_21(x)*c-y1_22(x)*c^2-y1_23(x)*c^3)/delta0_cb^2+y1_11(x)

eq1_9 := (diff(y1_9(x), x))/(27.1) = -2*b*c^2*phi_cb*(y1_17(x)+(1/2)*db*y1_19(x)-y1_20(x)-y1_21(x)*c-y1_22(x)*c^2-y1_23(x)*c^3)/delta0_cb^2+2*y1_12(x)

eq1_10 := (diff(y1_10(x), x))/(27.1) = -2*b*c^3*phi_cb*(y1_17(x)+(1/2)*db*y1_19(x)-y1_20(x)-y1_21(x)*c-y1_22(x)*c^2-y1_23(x)*c^3)/delta0_cb^2+3*y1_13(x)

eq1_11 := (diff(y1_11(x), x))/(27.1) = -b*phi_cb*(y1_18(x)-y1_24(x)-y1_25(x)*c-y1_26(x)*c^2)/delta0_cb^2

eq1_12 := (diff(y1_12(x), x))/(27.1) = -b*c*phi_cb*(y1_18(x)-y1_24(x)-y1_25(x)*c-y1_26(x)*c^2)/delta0_cb^2+2*Ec*b*c*y1_25(x)+nuc*y1_7(x)

eq1_13 := (diff(y1_13(x), x))/(27.1) = -b*c^2*phi_cb*(y1_18(x)-y1_24(x)-y1_25(x)*c-y1_26(x)*c^2)/delta0_cb^2+2*(4*Ec*b*c^3*y1_26(x)*(1/3)+nuc*y1_8(x))

eq1_14 := (diff(y1_14(x), x))/(27.1) = y1_1(x)/(b*A11)

eq1_15 := (diff(y1_15(x), x))/(27.1) = y1_16(x)+y1_2(x)/(b*A55)

eq1_16 := (diff(y1_16(x), x))/(27.1) = -y1_3(x)/(b*D11)

eq1_17 := (diff(y1_17(x), x))/(27.1) = y1_4(x)/(b*A11) 

eq1_18 := (diff(y1_18(x), x))/(27.1) = y1_19(x)+y1_5(x)/(b*A55)

eq1_19 := (diff(y1_19(x), x))/(27.1) = -y1_6(x)/(b*D11)

eq1_20 := (diff(y1_20(x), x))/(27.1) = (-nuc^2+1)*(9*c^2*y1_7(x)-15*y1_9(x))/(8*Ec*b*c^3)-nuc*y1_25(x)

eq1_21 := (diff(y1_21(x), x))/(27.1) = (-nuc^2+1)*(75*c^2*y1_8(x)-105*y1_10(x))/(8*Ec*b*c^5)-2*nuc*y1_26(x)

eq1_22 := (diff(y1_22(x), x))/(27.1) = (-nuc^2+1)*(45*y1_9(x)-15*c^2*y1_7(x))/(8*Ec*b*c^5)

eq1_23 := (diff(y1_23(x), x))/(27.1) = (-nuc^2+1)*(175*y1_10(x)-105*c^2*y1_8(x))/(8*Ec*b*c^7)

eq1_24 := (diff(y1_24(x), x))/(27.1) = (9*c^2*y1_11(x)-15*y1_13(x))/(8*Gc*b*c^3)-y1_21(x)

eq1_25 := (diff(y1_25(x), x))/(27.1) = 3*y1_12(x)/(2*Gc*b*c^3)-2*y1_22(x)

eq1_26 := (diff(y1_26(x), x))/(27.1) = (45*y1_13(x)-15*c^2*y1_11(x))/(8*Gc*b*c^5)-3*y1_23(x)

NULL

eq2_1 := (diff(y2_1(x), x))/(23.7) = 0; eq2_2 := (diff(y2_2(x), x))/(23.7) = 0; eq2_3 := (diff(y2_3(x), x))/(23.7) = y2_2(x); eq2_4 := (diff(y2_4(x), x))/(23.7) = 2*b*phi_cb*(y2_17(x)+(1/2)*db*y2_19(x)-y2_20(x)-y2_21(x)*c-y2_22(x)*c^2-y2_23(x)*c^3)/delta0_cb^2; eq2_5 := (diff(y2_5(x), x))/(23.7) = b*phi_cb*(y2_18(x)-y2_24(x)-y2_25(x)*c-y2_26(x)*c^2)/delta0_cb^2; eq2_6 := (diff(y2_6(x), x))/(23.7) = y2_5(x)-2*b*db*phi_cb*(y2_17(x)+(1/2)*db*y2_19(x)-y2_20(x)-y2_21(x)*c-y2_22(x)*c^2-y2_23(x)*c^3)/delta0_cb^2*(1/2); eq2_7 := (diff(y2_7(x), x))/(23.7) = -2*b*phi_cb*(y2_17(x)+(1/2)*db*y2_19(x)-y2_20(x)-y2_21(x)*c-y2_22(x)*c^2-y2_23(x)*c^3)/delta0_cb^2; eq2_8 := (diff(y2_8(x), x))/(23.7) = -2*b*c*phi_cb*(y2_17(x)+(1/2)*db*y2_19(x)-y2_20(x)-y2_21(x)*c-y2_22(x)*c^2-y2_23(x)*c^3)/delta0_cb^2+y2_11(x); eq2_9 := (diff(y2_9(x), x))/(23.7) = -2*b*c^2*phi_cb*(y2_17(x)+(1/2)*db*y2_19(x)-y2_20(x)-y2_21(x)*c-y2_22(x)*c^2-y2_23(x)*c^3)/delta0_cb^2+2*y2_12(x); eq2_10 := (diff(y2_10(x), x))/(23.7) = -2*b*c^3*phi_cb*(y2_17(x)+(1/2)*db*y2_19(x)-y2_20(x)-y2_21(x)*c-y2_22(x)*c^2-y2_23(x)*c^3)/delta0_cb^2+3*y2_13(x); eq2_11 := (diff(y2_11(x), x))/(23.7) = -b*phi_cb*(y2_18(x)-y2_24(x)-y2_25(x)*c-y2_26(x)*c^2)/delta0_cb^2; eq2_12 := (diff(y2_12(x), x))/(23.7) = -b*c*phi_cb*(y2_18(x)-y2_24(x)-y2_25(x)*c-y2_26(x)*c^2)/delta0_cb^2+2*Ec*b*c*y2_25(x)+nuc*y2_7(x); eq2_13 := (diff(y2_13(x), x))/(23.7) = -b*c^2*phi_cb*(y2_18(x)-y2_24(x)-y2_25(x)*c-y2_26(x)*c^2)/delta0_cb^2+2*(4*Ec*b*c^3*y2_26(x)*(1/3)+nuc*y2_8(x)); eq2_14 := (diff(y2_14(x), x))/(23.7) = y2_1(x)/(b*A11); eq2_15 := (diff(y2_15(x), x))/(23.7) = y2_16(x)+y2_2(x)/(b*A55); eq2_16 := (diff(y2_16(x), x))/(23.7) = -y2_3(x)/(b*D11); eq2_17 := (diff(y2_17(x), x))/(23.7) = y2_4(x)/(b*A11); eq2_18 := (diff(y2_18(x), x))/(23.7) = y2_19(x)+y2_5(x)/(b*A55); eq2_19 := (diff(y2_19(x), x))/(23.7) = -y2_6(x)/(b*D11); eq2_20 := (diff(y2_20(x), x))/(23.7) = (-nuc^2+1)*(9*c^2*y2_7(x)-15*y2_9(x))/(8*Ec*b*c^3)-nuc*y2_25(x); eq2_21 := (diff(y2_21(x), x))/(23.7) = (-nuc^2+1)*(75*c^2*y2_8(x)-105*y2_10(x))/(8*Ec*b*c^5)-2*nuc*y2_26(x); eq2_22 := (diff(y2_22(x), x))/(23.7) = (-nuc^2+1)*(45*y2_9(x)-15*c^2*y2_7(x))/(8*Ec*b*c^5); eq2_23 := (diff(y2_23(x), x))/(23.7) = (-nuc^2+1)*(175*y2_10(x)-105*c^2*y2_8(x))/(8*Ec*b*c^7); eq2_24 := (diff(y2_24(x), x))/(23.7) = (9*c^2*y2_11(x)-15*y2_13(x))/(8*Gc*b*c^3)-y2_21(x); eq2_25 := (diff(y2_25(x), x))/(23.7) = 3*y2_12(x)/(2*Gc*b*c^3)-2*y2_22(x); eq2_26 := (diff(y2_26(x), x))/(23.7) = (45*y2_13(x)-15*c^2*y2_11(x))/(8*Gc*b*c^5)-3*y2_23(x)

NULL

eq3_1 := (diff(y3_1(x), x))/(101.6) = -2*b*phi_ct*(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))*(1+(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct)*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2

eq3_2 := (diff(y3_2(x), x))/(101.6) = -b*phi_ct*(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2

eq3_3 := (diff(y3_3(x), x))/(101.6) = y3_2(x)-2*b*dt*phi_ct*(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))*(1+(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct)*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2*(1/2)

eq3_4 := (diff(y3_4(x), x))/(101.6) = 2*b*phi_cb*(y3_17(x)+(1/2)*db*y3_19(x)-y3_20(x)-y3_21(x)*c-y3_22(x)*c^2-y3_23(x)*c^3)/delta0_cb^2

eq3_5 := (diff(y3_5(x), x))/(101.6) = b*phi_cb*(y3_18(x)-y3_24(x)-y3_25(x)*c-y3_26(x)*c^2)/delta0_cb^2

eq3_6 := (diff(y3_6(x), x))/(101.6) = y3_5(x)-2*b*db*phi_cb*(y3_17(x)+(1/2)*db*y3_19(x)-y3_20(x)-y3_21(x)*c-y3_22(x)*c^2-y3_23(x)*c^3)/delta0_cb^2*(1/2)

eq3_7 := (diff(y3_7(x), x))/(101.6) = -2*b*phi_cb*(y3_17(x)+(1/2)*db*y3_19(x)-y3_20(x)-y3_21(x)*c-y3_22(x)*c^2-y3_23(x)*c^3)/delta0_cb^2+2*b*phi_ct*(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))*(1+(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct)*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2

eq3_8 := (diff(y3_8(x), x))/(101.6) = -2*b*c*phi_ct*(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))*(1+(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct)*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2-2*b*c*phi_cb*(y3_17(x)+(1/2)*db*y3_19(x)-y3_20(x)-y3_21(x)*c-y3_22(x)*c^2-y3_23(x)*c^3)/delta0_cb^2+y3_11(x)

eq3_9 := (diff(y3_9(x), x))/(101.6) = -2*b*c^2*phi_cb*(y3_17(x)+(1/2)*db*y3_19(x)-y3_20(x)-y3_21(x)*c-y3_22(x)*c^2-y3_23(x)*c^3)/delta0_cb^2+2*b*c^2*phi_ct*(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))*(1+(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct)*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2+2*y3_12(x)

eq3_10 := (diff(y3_10(x), x))/(101.6) = -2*b*c^3*phi_ct*(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))*(1+(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct)*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2-2*b*c^3*phi_cb*(y3_17(x)+(1/2)*db*y3_19(x)-y3_20(x)-y3_21(x)*c-y3_22(x)*c^2-y3_23(x)*c^3)/delta0_cb^2+3*y3_13(x)

eq3_11 := (diff(y3_11(x), x))/(101.6) = -b*phi_cb*(y3_18(x)-y3_24(x)-y3_25(x)*c-y3_26(x)*c^2)/delta0_cb^2+b*phi_ct*(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2

eq3_12 := (diff(y3_12(x), x))/(101.6) = -b*c*phi_ct*(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2-b*c*phi_cb*(y3_18(x)-y3_24(x)-y3_25(x)*c-y3_26(x)*c^2)/delta0_cb^2+2*Ec*b*c*y3_25(x)+nuc*y3_7(x)

eq3_13 := (diff(y3_13(x), x))/(101.6) = -b*c^2*phi_cb*(y3_18(x)-y3_24(x)-y3_25(x)*c-y3_26(x)*c^2)/delta0_cb^2+b*c^2*phi_ct*(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))*exp(-(y3_24(x)-y3_25(x)*c+y3_26(x)*c^2-y3_15(x))/delta0_ct-(y3_20(x)-y3_21(x)*c+y3_22(x)*c^2-y3_23(x)*c^3-y3_14(x)+(1/2)*dt*y3_16(x))^2/delta0_ct^2)/delta0_ct^2+2*(4*Ec*b*c^3*y3_26(x)*(1/3)+nuc*y3_8(x))

eq3_14 := (diff(y3_14(x), x))/(101.6) = y3_1(x)/(b*A11); eq3_15 := (diff(y3_15(x), x))/(101.6) = y3_16(x)+y3_2(x)/(b*A55); eq3_16 := (diff(y3_16(x), x))/(101.6) = -y3_3(x)/(b*D11); eq3_17 := (diff(y3_17(x), x))/(101.6) = y3_4(x)/(b*A11); eq3_18 := (diff(y3_18(x), x))/(101.6) = y3_19(x)+y3_5(x)/(b*A55); eq3_19 := (diff(y3_19(x), x))/(101.6) = -y3_6(x)/(b*D11); eq3_20 := (diff(y3_20(x), x))/(101.6) = (-nuc^2+1)*(9*c^2*y3_7(x)-15*y3_9(x))/(8*Ec*b*c^3)-nuc*y3_25(x); eq3_21 := (diff(y3_21(x), x))/(101.6) = (-nuc^2+1)*(75*c^2*y3_8(x)-105*y3_10(x))/(8*Ec*b*c^5)-2*nuc*y3_26(x); eq3_22 := (diff(y3_22(x), x))/(101.6) = (-nuc^2+1)*(45*y3_9(x)-15*c^2*y3_7(x))/(8*Ec*b*c^5); eq3_23 := (diff(y3_23(x), x))/(101.6) = (-nuc^2+1)*(175*y3_10(x)-105*c^2*y3_8(x))/(8*Ec*b*c^7); eq3_24 := (diff(y3_24(x), x))/(101.6) = (9*c^2*y3_11(x)-15*y3_13(x))/(8*Gc*b*c^3)-y3_21(x); eq3_25 := (diff(y3_25(x), x))/(101.6) = 3*y3_12(x)/(2*Gc*b*c^3)-2*y3_22(x); eq3_26 := (diff(y3_26(x), x))/(101.6) = (45*y3_13(x)-15*c^2*y3_11(x))/(8*Gc*b*c^5)-3*y3_23(x)

NULL

eq1_27 := diff(y1_27(x), x) = 0

``

NULL

odesys := [eq1_1, eq1_2, eq1_3, eq1_4, eq1_5, eq1_6, eq1_7, eq1_8, eq1_9, eq1_10, eq1_11, eq1_12, eq1_13, eq1_14, eq1_15, eq1_16, eq1_17, eq1_18, eq1_19, eq1_20, eq1_21, eq1_22, eq1_23, eq1_24, eq1_25, eq1_26, eq1_27, eq2_1, eq2_2, eq2_3, eq2_4, eq2_5, eq2_6, eq2_7, eq2_8, eq2_9, eq2_10, eq2_11, eq2_12, eq2_13, eq2_14, eq2_15, eq2_16, eq2_17, eq2_18, eq2_19, eq2_20, eq2_21, eq2_22, eq2_23, eq2_24, eq2_25, eq2_26, eq3_1, eq3_2, eq3_3, eq3_4, eq3_5, eq3_6, eq3_7, eq3_8, eq3_9, eq3_10, eq3_11, eq3_12, eq3_13, eq3_14, eq3_15, eq3_16, eq3_17, eq3_18, eq3_19, eq3_20, eq3_21, eq3_22, eq3_23, eq3_24, eq3_25, eq3_26, y1_1(0) = 0, y1_2(0) = 0, y1_3(0) = 0, y1_4(0) = 0, y1_5(0) = 0, y1_6(0) = 0, y1_7(0) = 0, y1_8(0) = 0, y1_9(0) = 0, y1_10(0) = 0, y1_11(0) = 0, y1_12(0) = 0, y1_13(0) = 0, y1_2(1)-y2_2(0)-y1_27 = 0, y1_3(1)-y2_3(0) = 0, y1_6(1)-y2_6(0) = 0, y1_7(1)-y2_7(0) = 0, y1_8(1)-y2_8(0) = 0, y1_9(1)-y2_9(0) = 0, y1_10(1)-y2_10(0) = 0, y1_11(1)-y2_11(0) = 0, y1_12(1)-y2_12(0) = 0, y1_13(1)-y2_13(0) = 0, y1_14(1) = y2_14(0), y1_15(1) = y2_15(0), y1_16(1) = y2_16(0), y1_17(1) = y2_17(0), y1_18(1) = y2_18(0), y1_19(1) = y2_19(0), y1_20(1) = y2_20(0), y1_21(1) = y2_21(0), y1_22(1) = y2_22(0), y1_23(1) = y2_23(0), y1_24(1) = y2_24(0), y1_25(1) = y2_25(0), y1_26(1) = y2_26(0), y1_14(1) = 0, y1_15(1) = -lambda, y1_17(1) = 0, y1_18(1) = 0, y2_1(1)-y3_1(0) = 0, y2_2(1)-y3_2(0) = 0, y2_3(1)-y3_3(0) = 0, y2_4(1)-y3_4(0) = 0, y2_5(1)-y3_5(0) = 0, y2_6(1)-y3_6(0) = 0, y2_7(1)-y3_7(0) = 0, y2_8(1)-y3_8(0) = 0, y2_9(1)-y3_9(0) = 0, y2_10(1)-y3_10(0) = 0, y2_11(1)-y3_11(0) = 0, y2_12(1)-y3_12(0) = 0, y2_13(1)-y3_13(0) = 0, y2_14(1) = y3_14(0), y2_15(1) = y3_15(0), y2_16(1) = y3_16(0), y2_17(1) = y3_17(0), y2_18(1) = y3_18(0), y2_19(1) = y3_19(0), y2_20(1) = y3_20(0), y2_21(1) = y3_21(0), y2_22(1) = y3_22(0), y2_23(1) = y3_23(0), y2_24(1) = y3_24(0), y2_25(1) = y3_25(0), y2_26(1) = y3_26(0), y3_1(1) = 0, y3_2(1) = 0, y3_3(1) = 0, y3_4(1) = 0, y3_5(1) = 0, y3_6(1) = 0, y3_7(1) = 0, y3_8(1) = 0, y3_9(1) = 0, y3_10(1) = 0, y3_11(1) = 0, y3_12(1) = 0, y3_13(1) = 0], [y1_1(x), y1_2(x), y1_3(x), y1_4(x), y1_5(x), y1_6(x), y1_7(x), y1_8(x), y1_9(x), y1_10(x), y1_11(x), y1_12(x), y1_13(x), y1_14(x), y1_15(x), y1_16(x), y1_17(x), y1_18(x), y1_19(x), y1_20(x), y1_21(x), y1_22(x), y1_23(x), y1_24(x), y1_25(x), y1_26(x), y1_27(x), y2_1(x), y2_2(x), y2_3(x), y2_4(x), y2_5(x), y2_6(x), y2_7(x), y2_8(x), y2_9(x), y2_10(x), y2_11(x), y2_12(x), y2_13(x), y2_14(x), y2_15(x), y2_16(x), y2_17(x), y2_18(x), y2_19(x), y2_20(x), y2_21(x), y2_22(x), y2_23(x), y2_24(x), y2_25(x), y2_26(x), y3_1(x), y3_2(x), y3_3(x), y3_4(x), y3_5(x), y3_6(x), y3_7(x), y3_8(x), y3_9(x), y3_10(x), y3_11(x), y3_12(x), y3_13(x), y3_14(x), y3_15(x), y3_16(x), y3_17(x), y3_18(x), y3_19(x), y3_20(x), y3_21(x), y3_22(x), y3_23(x), y3_24(x), y3_25(x), y3_26(x)]
``

solnumeric1 := dsolve(eval(odesys, lambda = .1), numeric)

solnumeric1(1)

[x = 1., y1_1(x) = HFloat(0.0), y1_2(x) = HFloat(-3.420984483671309e-31), y1_3(x) = HFloat(-9.266873428961157e-30), y1_4(x) = HFloat(-2.8265796778082293), y1_5(x) = HFloat(4.585656258526918), y1_6(x) = HFloat(29.73356557713309), y1_7(x) = HFloat(2.8265796778082293), y1_8(x) = HFloat(1.6414688476893966), y1_9(x) = HFloat(78.13532261272604), y1_10(x) = HFloat(134.29458529790443), y1_11(x) = HFloat(-4.585656258526918), y1_12(x) = HFloat(-12.84587679453849), y1_13(x) = HFloat(-141.040214445319), y1_14(x) = HFloat(0.0), y1_15(x) = HFloat(-0.09999999999999995), y1_16(x) = HFloat(0.004086375487434791), y1_17(x) = HFloat(0.0), y1_18(x) = HFloat(0.0), y1_19(x) = HFloat(-2.458791784171863e-4), y1_20(x) = HFloat(-0.001957542333257033), y1_21(x) = HFloat(1.9956662771272542e-4), y1_22(x) = HFloat(-1.5079719874448682e-6), y1_23(x) = HFloat(3.3273624675095036e-8), y1_24(x) = HFloat(-7.164450668824305e-4), y1_25(x) = HFloat(1.1943543119401487e-5), y1_26(x) = HFloat(1.0902056689470582e-6), y1_27(x) = HFloat(-11.236234112909697), y2_1(x) = HFloat(7.740295426479526e-11), y2_2(x) = HFloat(11.236234112909697), y2_3(x) = HFloat(266.2987484570526), y2_4(x) = HFloat(-13.486001518995927), y2_5(x) = HFloat(-0.02999714213517951), y2_6(x) = HFloat(-22.43800529304631), y2_7(x) = HFloat(13.486001518918538), y2_8(x) = HFloat(-94.16612630071822), y2_9(x) = HFloat(785.015722665203), y2_10(x) = HFloat(-6775.009356439386), y2_11(x) = HFloat(-11.206236970774473), y2_12(x) = HFloat(77.76745611873598), y2_13(x) = HFloat(-642.9939052968533), y2_14(x) = HFloat(4.689768932388945e-16), y2_15(x) = HFloat(-0.0187710269573956), y2_16(x) = HFloat(0.0020861962139203305), y2_17(x) = HFloat(-5.312274942589151e-5), y2_18(x) = HFloat(-0.005830316541460037), y2_19(x) = HFloat(-1.227143374954293e-4), y2_20(x) = HFloat(-0.0011198874012043685), y2_21(x) = HFloat(7.602496797322661e-5), y2_22(x) = HFloat(-2.932591600936832e-6), y2_23(x) = HFloat(5.109018589948548e-7), y2_24(x) = HFloat(-0.0066417834273388736), y2_25(x) = HFloat(1.1855664085790522e-4), y2_26(x) = HFloat(-6.471131431633229e-6), y3_1(x) = HFloat(0.0), y3_2(x) = HFloat(0.0), y3_3(x) = HFloat(0.0), y3_4(x) = HFloat(0.0), y3_5(x) = HFloat(0.0), y3_6(x) = HFloat(0.0), y3_7(x) = HFloat(0.0), y3_8(x) = HFloat(0.0), y3_9(x) = HFloat(0.0), y3_10(x) = HFloat(0.0), y3_11(x) = HFloat(0.0), y3_12(x) = HFloat(0.0), y3_13(x) = HFloat(0.0), y3_14(x) = HFloat(-7.007769576850207e-6), y3_15(x) = HFloat(-0.004786448735922665), y3_16(x) = HFloat(4.689280929440655e-6), y3_17(x) = HFloat(-1.115140128648413e-4), y3_18(x) = HFloat(-0.004788926306328681), y3_19(x) = HFloat(4.716347003469115e-6), y3_20(x) = HFloat(-5.966958988267552e-5), y3_21(x) = HFloat(-4.710193875348028e-6), y3_22(x) = HFloat(3.723242564645378e-9), y3_23(x) = HFloat(1.7383541695404693e-11), y3_24(x) = HFloat(-0.004787700776648968), y3_25(x) = HFloat(-8.409570657167222e-8), y3_26(x) = HFloat(3.920653061680723e-13)]

(1)

[x = 1., y1_1(x) = HFloat(0.0), y1_2(x) = HFloat(-3.420984483671309e-31), y1_3(x) = HFloat(-9.266873428961157e-30), y1_4(x) = HFloat(-2.8265796778082293), y1_5(x) = HFloat(4.585656258526918), y1_6(x) = HFloat(29.73356557713309), y1_7(x) = HFloat(2.8265796778082293), y1_8(x) = HFloat(1.6414688476893966), y1_9(x) = HFloat(78.13532261272604), y1_10(x) = HFloat(134.29458529790443), y1_11(x) = HFloat(-4.585656258526918), y1_12(x) = HFloat(-12.84587679453849), y1_13(x) = HFloat(-141.040214445319), y1_14(x) = HFloat(0.0), y1_15(x) = HFloat(-0.09999999999999995), y1_16(x) = HFloat(0.004086375487434791), y1_17(x) = HFloat(0.0), y1_18(x) = HFloat(0.0), y1_19(x) = HFloat(-2.458791784171863e-4), y1_20(x) = HFloat(-0.001957542333257033), y1_21(x) = HFloat(1.9956662771272542e-4), y1_22(x) = HFloat(-1.5079719874448682e-6), y1_23(x) = HFloat(3.3273624675095036e-8), y1_24(x) = HFloat(-7.164450668824305e-4), y1_25(x) = HFloat(1.1943543119401487e-5), y1_26(x) = HFloat(1.0902056689470582e-6), y1_27(x) = HFloat(-11.236234112909697), y2_1(x) = HFloat(7.740295426479526e-11), y2_2(x) = HFloat(11.236234112909697), y2_3(x) = HFloat(266.2987484570526), y2_4(x) = HFloat(-13.486001518995927), y2_5(x) = HFloat(-0.02999714213517951), y2_6(x) = HFloat(-22.43800529304631), y2_7(x) = HFloat(13.486001518918538), y2_8(x) = HFloat(-94.16612630071822), y2_9(x) = HFloat(785.015722665203), y2_10(x) = HFloat(-6775.009356439386), y2_11(x) = HFloat(-11.206236970774473), y2_12(x) = HFloat(77.76745611873598), y2_13(x) = HFloat(-642.9939052968533), y2_14(x) = HFloat(4.689768932388945e-16), y2_15(x) = HFloat(-0.0187710269573956), y2_16(x) = HFloat(0.0020861962139203305), y2_17(x) = HFloat(-5.312274942589151e-5), y2_18(x) = HFloat(-0.005830316541460037), y2_19(x) = HFloat(-1.227143374954293e-4), y2_20(x) = HFloat(-0.0011198874012043685), y2_21(x) = HFloat(7.602496797322661e-5), y2_22(x) = HFloat(-2.932591600936832e-6), y2_23(x) = HFloat(5.109018589948548e-7), y2_24(x) = HFloat(-0.0066417834273388736), y2_25(x) = HFloat(1.1855664085790522e-4), y2_26(x) = HFloat(-6.471131431633229e-6), y3_1(x) = HFloat(0.0), y3_2(x) = HFloat(0.0), y3_3(x) = HFloat(0.0), y3_4(x) = HFloat(0.0), y3_5(x) = HFloat(0.0), y3_6(x) = HFloat(0.0), y3_7(x) = HFloat(0.0), y3_8(x) = HFloat(0.0), y3_9(x) = HFloat(0.0), y3_10(x) = HFloat(0.0), y3_11(x) = HFloat(0.0), y3_12(x) = HFloat(0.0), y3_13(x) = HFloat(0.0), y3_14(x) = HFloat(-7.007769576850207e-6), y3_15(x) = HFloat(-0.004786448735922665), y3_16(x) = HFloat(4.689280929440655e-6), y3_17(x) = HFloat(-1.115140128648413e-4), y3_18(x) = HFloat(-0.004788926306328681), y3_19(x) = HFloat(4.716347003469115e-6), y3_20(x) = HFloat(-5.966958988267552e-5), y3_21(x) = HFloat(-4.710193875348028e-6), y3_22(x) = HFloat(3.723242564645378e-9), y3_23(x) = HFloat(1.7383541695404693e-11), y3_24(x) = HFloat(-0.004787700776648968), y3_25(x) = HFloat(-8.409570657167222e-8), y3_26(x) = HFloat(3.920653061680723e-13)]
``

(2)

solnumeric2 := dsolve(eval(odesys, lambda = .2), numeric, approxsoln = solnumeric1)

solnumeric2(1)

[x = 1., y1_1(x) = HFloat(0.0), y1_2(x) = HFloat(-3.030574676828747e-31), y1_3(x) = HFloat(-8.154092896420496e-30), y1_4(x) = HFloat(-2.8265796778082306), y1_5(x) = HFloat(4.585656258526919), y1_6(x) = HFloat(29.733565577133103), y1_7(x) = HFloat(2.8265796778082306), y1_8(x) = HFloat(1.6414688476894284), y1_9(x) = HFloat(78.13532261272591), y1_10(x) = HFloat(134.2945852979062), y1_11(x) = HFloat(-4.585656258526919), y1_12(x) = HFloat(-12.845876794538517), y1_13(x) = HFloat(-141.04021444531904), y1_14(x) = HFloat(0.0), y1_15(x) = HFloat(-0.09999999999999996), y1_16(x) = HFloat(0.004086375487434796), y1_17(x) = HFloat(0.0), y1_18(x) = HFloat(0.0), y1_19(x) = HFloat(-2.458791784171864e-4), y1_20(x) = HFloat(-0.0019575423332570297), y1_21(x) = HFloat(1.9956662771272528e-4), y1_22(x) = HFloat(-1.50797198744487e-6), y1_23(x) = HFloat(3.3273624675095235e-8), y1_24(x) = HFloat(-7.164450668824302e-4), y1_25(x) = HFloat(1.1943543119401463e-5), y1_26(x) = HFloat(1.0902056689470559e-6), y1_27(x) = HFloat(-11.2362341129097), y2_1(x) = HFloat(7.741244519455937e-11), y2_2(x) = HFloat(11.2362341129097), y2_3(x) = HFloat(266.2987484570527), y2_4(x) = HFloat(-13.486001518995932), y2_5(x) = HFloat(-0.029997142135179216), y2_6(x) = HFloat(-22.438005293046288), y2_7(x) = HFloat(13.486001518918526), y2_8(x) = HFloat(-94.16612630071823), y2_9(x) = HFloat(785.0157226652033), y2_10(x) = HFloat(-6775.0093564393865), y2_11(x) = HFloat(-11.206236970774471), y2_12(x) = HFloat(77.76745611873608), y2_13(x) = HFloat(-642.9939052968541), y2_14(x) = HFloat(4.690343978496297e-16), y2_15(x) = HFloat(-0.018771026957395617), y2_16(x) = HFloat(0.002086196213920331), y2_17(x) = HFloat(-5.312274942589153e-5), y2_18(x) = HFloat(-0.005830316541460037), y2_19(x) = HFloat(-1.2271433749542907e-4), y2_20(x) = HFloat(-0.001119887401204369), y2_21(x) = HFloat(7.602496797322656e-5), y2_22(x) = HFloat(-2.9325916009368313e-6), y2_23(x) = HFloat(5.109018589948565e-7), y2_24(x) = HFloat(-0.006641783427338871), y2_25(x) = HFloat(1.1855664085790523e-4), y2_26(x) = HFloat(-6.471131431633234e-6), y3_1(x) = HFloat(0.0), y3_2(x) = HFloat(0.0), y3_3(x) = HFloat(0.0), y3_4(x) = HFloat(0.0), y3_5(x) = HFloat(0.0), y3_6(x) = HFloat(0.0), y3_7(x) = HFloat(0.0), y3_8(x) = HFloat(0.0), y3_9(x) = HFloat(0.0), y3_10(x) = HFloat(0.0), y3_11(x) = HFloat(0.0), y3_12(x) = HFloat(0.0), y3_13(x) = HFloat(0.0), y3_14(x) = HFloat(-7.007769576850145e-6), y3_15(x) = HFloat(-0.00478644873592267), y3_16(x) = HFloat(4.6892809294406756e-6), y3_17(x) = HFloat(-1.115140128648413e-4), y3_18(x) = HFloat(-0.004788926306328684), y3_19(x) = HFloat(4.716347003469143e-6), y3_20(x) = HFloat(-5.966958988267562e-5), y3_21(x) = HFloat(-4.710193875348042e-6), y3_22(x) = HFloat(3.7232425646457026e-9), y3_23(x) = HFloat(1.7383541695557647e-11), y3_24(x) = HFloat(-0.004787700776648964), y3_25(x) = HFloat(-8.409570657166629e-8), y3_26(x) = HFloat(3.9206530681961753e-13)]

(3)

NULL

``


 

Download toughODE.mw

 

Four possible alternatives are shown in the attached

restart

with(Student:-VectorCalculus)

BasisFormat(true)

v := VectorField(`<,>`(x, y, z))

Vector(3, {(1) = x, (2) = y, (3) = z})

(1)

ScalarPotential(v)

(1/2)*x^2+(1/2)*y^2+(1/2)*z^2

(2)

 

Thats easy to use the ScalarPotential for a manual exercise

Example 2: finding a Potential function

Show that F = [exp(x)*cos(y)+y*z, x*z-exp(x)*sin(y), x*y+z]

F = [exp(x)*cos(y)+y*z, x*z-exp(x)*sin(y), x*y+z]

(3)

 is conversative and find a potential function for it.

F := VectorField(`<,>`(exp(x)*cos(y)+y*z, x*z-exp(x)*sin(y), x*y+z))

Vector(3, {(1) = exp(x)*cos(y)+y*z, (2) = x*z-exp(x)*sin(y), (3) = x*y+z})

(4)

"(->)"

y*z*x+exp(x)*cos(y)+(1/2)*z^2

(5)

ScalarPotential(F); pdsolve(F = Gradient(SP(x, y, z))); pdsolve(F = Nabla(SP(x, y, z))); eqs := [seq(F[i] = Nabla(SP(x, y, z))[i], i = 1 .. 3)]; pdsolve(eqs)

{SP(x, y, z) = (1/2)*z^2+y*z*x+exp(x)*cos(y)+_C1}

(6)

NULL

``


 

Download SP2.mw

If I open the help page DeepLearning:-Classify as a worksheet and then execcute it repeatedly - it always *seems* to work.

See the attached (which for some reason won't display inline here,  not sure why??)

So I think you will have to provide the worksheet which exhibits the failure (Use the big green  up-arrow in the Mapleprimes toolbar)

BTW I'm using 64-bit Win7

Download DLClass.mw

You just need to perform computations to much higher accuracy - in the attached I used Digits=50 (whihc may(?) be overkill

It is also probably a good idea to use RootFinding:-Nextzero() rather than RootFinding:-Analytic(), since you are only searching for roots along the real line.

It is also much simpler to use the LinearAlgebra:-GenerateMatrix() command to construct the  matrix whose determinant is needed to obtain the roots.

See the attached

restart; Digits := 50; interface(displayprecision = 10); with(LinearAlgebra); with(RootFinding)

Phi := unapply(C1*cos(beta*x)+C2*sin(beta*x)+C3*cosh(beta*x)+C4*sinh(beta*x), x); bc1 := EI*((D@@2)(Phi))(0)-Kr*(D(Phi))(0); bc2 := EI*((D@@3)(Phi))(0)+Kt*Phi(0); bc3 := ((D@@2)(Phi))(L)/beta^2; bc4 := ((D@@3)(Phi))(L)/beta^3; M := GenerateMatrix([bc1, bc2, bc3, bc4], [C1, C2, C3, C4])[1]; detM := eval(simplify(Determinant(M)), [EI = 20*10^9*7243.32, L = 121, Kt = 4.5*10^9, Kr = 5.7*10^9, rho = 400, A = 30]); rts := Array(0 .. 5); resids := Array(1 .. 5); for i to 5 do rts[i] := NextZero(unapply(detM, beta), rts[i-1]); resids[i] := eval(detM, beta = rts[i]) end do; rts[1 .. 5]; resids[1 .. 5]

bc1 := EI*(-C1*beta^2+C3*beta^2)-Kr*(C2*beta+C4*beta)

 

bc2 := EI*(-C2*beta^3+C4*beta^3)+Kt*(C1+C3)

 

bc3 := (-C1*beta^2*cos(beta*L)-C2*beta^2*sin(beta*L)+C3*beta^2*cosh(beta*L)+C4*beta^2*sinh(beta*L))/beta^2

 

bc4 := (C1*beta^3*sin(beta*L)-C2*beta^3*cos(beta*L)+C3*beta^3*sinh(beta*L)+C4*beta^3*cosh(beta*L))/beta^3

 

Matrix(4, 4, {(1, 1) = -EI*beta^2, (1, 2) = -Kr*beta, (1, 3) = EI*beta^2, (1, 4) = -Kr*beta, (2, 1) = Kt, (2, 2) = -EI*beta^3, (2, 3) = Kt, (2, 4) = EI*beta^3, (3, 1) = -cos(beta*L), (3, 2) = -sin(beta*L), (3, 3) = cosh(beta*L), (3, 4) = sinh(beta*L), (4, 1) = sin(beta*L), (4, 2) = -cos(beta*L), (4, 3) = sinh(beta*L), (4, 4) = cosh(beta*L)})

 

Vector[row](5, {(1) = 0.28561117575040674433472658728587202863832735961218e-2, (2) = 0.27953120389950134576687262910544157423896370383064e-1, (3) = 0.43550314889288492708146739421724924920212998363216e-1, (4) = 0.65908358098826847376941652358667073390929887695958e-1, (5) = 0.91228672905353004995231758761707573469828658718856e-1})

 

Array(%id = 18446744074366223774)

(1)

``

 


 

Download roots.mw

you might want to consider the attached very carefully

  restart:
#
# Define some parameters
#
   params:=[phi1 = .1, phi2 = .1, rhos1 = 2720, rhos2 = 2810,
            rhosf = 997.1, khnf = 1.083061737, kf = .613,
            cp1 = 893, cp2 = 960, cpf = 4179, Pr = 6.2,
            knf = .8154646474, S=0.5, R=0.5, gamma=0.5
          ]:
#
# Define the ODES
#
  ODES := (diff(f(eta), `$`(eta, 4)))/((1-phi1)^2.5*(1-phi2)^2.5*((1-phi2)*(1-phi1+phi1*rhos1/rhosf)+phi2*rhos2/rhosf))+S*(f(eta)*(diff(f(eta), `$`(eta, 3)))-3*(diff(f(eta), `$`(eta, 2)))-eta*(diff(f(eta), `$`(eta, 3)))-(diff(f(eta), eta))*(diff(f(eta), `$`(eta, 2)))) = 0,
          (khnf/kf+(4/3)*R)*(diff(theta(eta), `$`(eta, 2)))/((1-phi2)*(1-phi1+phi1*rhos1*cp1/(rhosf*cpf))+phi2*rhos2*cp2/(rhosf*cpf))+S*Pr*(f(eta)*(diff(theta(eta), eta))-eta*(diff(theta(eta), eta))-gamma*(eta^2*(diff(theta(eta), `$`(eta, 2)))-2*eta*f(eta)*(diff(theta(eta), `$`(eta, 2)))-eta*(diff(f(eta), eta))*(diff(theta(eta), eta))+f(eta)*(diff(f(eta), eta))*(diff(theta(eta), eta))+f(eta)^2*(diff(theta(eta), `$`(eta, 2))))) = 0:
#
# Define the boundary conditions. Notice that in the OPs
# original definition ((D^2)(f))(0) = 0 is a a bit dumb.
# Apart from anything else it could be more simply written
# written as D(f)(0) = 0
#
# Can only suggest that OP reads the output of these three
# commands very carefully and decide which is required

  BCs0:= f(0) = 0, ((D^2)(f))(0) = 0, (D(theta))(0) = 0, f(1) = 0, (D(f))(1) = 0, theta(1) = 1;
  BCs1:= f(0) = 0, D(f)(0) = 0, D(theta)(0) = 0, f(1) = 0, D(f)(1) = 0, theta(1) = 1;
  BCs2:= f(0) = 0, D[1,1](f)(0) = 0, D(theta)(0) = 0, f(1) = 0, D(f)(1) = 0, theta(1) = 1;

f(0) = 0, (D(f))(0)^2 = 0, (D(theta))(0) = 0, f(1) = 0, (D(f))(1) = 0, theta(1) = 1

 

f(0) = 0, (D(f))(0) = 0, (D(theta))(0) = 0, f(1) = 0, (D(f))(1) = 0, theta(1) = 1

 

f(0) = 0, ((D@@2)(f))(0) = 0, (D(theta))(0) = 0, f(1) = 0, (D(f))(1) = 0, theta(1) = 1

(1)

#
# Since I don't think the first set of BCs above is plausible
# because of the existence of ((D^2)(f))(0) = 0, just solve the
# system with the second and third set of boundary conditions.
# Neither produces very interesting solutions
#
  sol1:=dsolve( [eval([ODES], params)[], BCs1], numeric);
  sol2:=dsolve( [eval([ODES], params)[], BCs2], numeric);
  plots:-odeplot(sol1, [[eta, f(eta)], [eta,theta(eta)]], eta=0..1, axes=boxed, color=[red,blue]);
  plots:-odeplot(sol2, [[eta, f(eta)], [eta,theta(eta)]], eta=0..1, axes=boxed, color=[red,blue]);

proc (x_bvp) local res, data, solnproc, _ndsol, outpoint, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; _EnvDSNumericSaveDigits := Digits; Digits := 15; if _EnvInFsolve = true then outpoint := evalf[_EnvDSNumericSaveDigits](x_bvp) else outpoint := evalf(x_bvp) end if; data := Array(1..4, {(1) = proc (outpoint) local X, Y, YP, yout, errproc, L, V, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; X := Vector(8, {(1) = .0, (2) = .14285714285714277, (3) = .28571428571428553, (4) = .4285714285714284, (5) = .5714285714285713, (6) = .7142857142857142, (7) = .857142857142857, (8) = 1.0}, datatype = float[8], order = C_order); Y := Matrix(8, 6, {(1, 1) = .0, (1, 2) = .0, (1, 3) = .0, (1, 4) = .0, (1, 5) = 1.0, (1, 6) = .0, (2, 1) = .0, (2, 2) = .0, (2, 3) = .0, (2, 4) = .0, (2, 5) = 1.0, (2, 6) = .0, (3, 1) = .0, (3, 2) = .0, (3, 3) = .0, (3, 4) = .0, (3, 5) = 1.0, (3, 6) = .0, (4, 1) = .0, (4, 2) = .0, (4, 3) = .0, (4, 4) = .0, (4, 5) = 1.0, (4, 6) = .0, (5, 1) = .0, (5, 2) = .0, (5, 3) = .0, (5, 4) = .0, (5, 5) = 1.0, (5, 6) = .0, (6, 1) = .0, (6, 2) = .0, (6, 3) = .0, (6, 4) = .0, (6, 5) = 1.0, (6, 6) = .0, (7, 1) = .0, (7, 2) = .0, (7, 3) = .0, (7, 4) = .0, (7, 5) = 1.0, (7, 6) = .0, (8, 1) = .0, (8, 2) = .0, (8, 3) = .0, (8, 4) = .0, (8, 5) = 1.0, (8, 6) = .0}, datatype = float[8], order = C_order); YP := Matrix(8, 6, {(1, 1) = .0, (1, 2) = .0, (1, 3) = .0, (1, 4) = .0, (1, 5) = .0, (1, 6) = .0, (2, 1) = .0, (2, 2) = .0, (2, 3) = .0, (2, 4) = .0, (2, 5) = .0, (2, 6) = .0, (3, 1) = .0, (3, 2) = .0, (3, 3) = .0, (3, 4) = .0, (3, 5) = .0, (3, 6) = .0, (4, 1) = .0, (4, 2) = .0, (4, 3) = .0, (4, 4) = .0, (4, 5) = .0, (4, 6) = .0, (5, 1) = .0, (5, 2) = .0, (5, 3) = .0, (5, 4) = .0, (5, 5) = .0, (5, 6) = .0, (6, 1) = .0, (6, 2) = .0, (6, 3) = .0, (6, 4) = .0, (6, 5) = .0, (6, 6) = .0, (7, 1) = .0, (7, 2) = .0, (7, 3) = .0, (7, 4) = .0, (7, 5) = .0, (7, 6) = .0, (8, 1) = .0, (8, 2) = .0, (8, 3) = .0, (8, 4) = .0, (8, 5) = .0, (8, 6) = .0}, datatype = float[8], order = C_order); errproc := proc (x_bvp) local outpoint, X, Y, yout, L, V, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; Digits := 15; outpoint := evalf(x_bvp); X := Vector(8, {(1) = .0, (2) = .14285714285714277, (3) = .28571428571428553, (4) = .4285714285714284, (5) = .5714285714285713, (6) = .7142857142857142, (7) = .857142857142857, (8) = 1.0}, datatype = float[8], order = C_order); Y := Matrix(8, 6, {(1, 1) = .0, (1, 2) = .0, (1, 3) = .0, (1, 4) = .0, (1, 5) = .0, (1, 6) = .0, (2, 1) = .0, (2, 2) = .0, (2, 3) = .0, (2, 4) = .0, (2, 5) = .0, (2, 6) = .0, (3, 1) = .0, (3, 2) = .0, (3, 3) = .0, (3, 4) = .0, (3, 5) = .0, (3, 6) = .0, (4, 1) = .0, (4, 2) = .0, (4, 3) = .0, (4, 4) = .0, (4, 5) = .0, (4, 6) = .0, (5, 1) = .0, (5, 2) = .0, (5, 3) = .0, (5, 4) = .0, (5, 5) = .0, (5, 6) = .0, (6, 1) = .0, (6, 2) = .0, (6, 3) = .0, (6, 4) = .0, (6, 5) = .0, (6, 6) = .0, (7, 1) = .0, (7, 2) = .0, (7, 3) = .0, (7, 4) = .0, (7, 5) = .0, (7, 6) = .0, (8, 1) = .0, (8, 2) = .0, (8, 3) = .0, (8, 4) = .0, (8, 5) = .0, (8, 6) = .0}, datatype = float[8], order = C_order); if not type(outpoint, 'numeric') then if outpoint = "start" or outpoint = "left" then return X[1] elif outpoint = "right" then return X[8] elif outpoint = "order" then return 2 elif outpoint = "error" then return HFloat(-0.0) elif outpoint = "errorproc" then error "this is already the error procedure" elif outpoint = "rawdata" then return [6, 8, [f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)], X, Y] else return ('procname')(x_bvp) end if end if; if outpoint < X[1] or X[8] < outpoint then error "solution is only defined in the range %1..%2", X[1], X[8] end if; V := array([1 = 4, 2 = 0]); if Digits <= trunc(evalhf(Digits)) then L := Vector(4, 'datatype' = 'float'[8]); yout := Vector(6, 'datatype' = 'float'[8]); evalhf(`dsolve/numeric/lagrange`(8, 6, X, Y, outpoint, var(yout), var(L), var(V))) else L := Vector(4, 'datatype' = 'sfloat'); yout := Vector(6, 'datatype' = 'sfloat'); `dsolve/numeric/lagrange`(8, 6, X, Y, outpoint, yout, L, V) end if; [eta = outpoint, seq('[f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)]'[i] = yout[i], i = 1 .. 6)] end proc; if not type(outpoint, 'numeric') then if outpoint = "start" or outpoint = "left" then return X[1] elif outpoint = "method" then return "bvp" elif outpoint = "right" then return X[8] elif outpoint = "order" then return 2 elif outpoint = "error" then return HFloat(-0.0) elif outpoint = "errorproc" then return eval(errproc) elif outpoint = "rawdata" then return [6, 8, "depnames", X, Y, YP] else error "non-numeric value" end if end if; if outpoint < X[1] or X[8] < outpoint then error "solution is only defined in the range %1..%2", X[1], X[8] end if; if Digits <= trunc(evalhf(Digits)) and (_EnvInFsolve <> true or _EnvDSNumericSaveDigits <= trunc(evalhf(Digits))) then V := array( 1 .. 6, [( 1 ) = (7), ( 2 ) = (0), ( 3 ) = (false), ( 4 ) = (false), ( 5 ) = (false), ( 6 ) = (false)  ] ); L := Matrix(7, 2, {(1, 1) = .0, (1, 2) = .0, (2, 1) = .0, (2, 2) = .0, (3, 1) = .0, (3, 2) = .0, (4, 1) = .0, (4, 2) = .0, (5, 1) = .0, (5, 2) = .0, (6, 1) = .0, (6, 2) = .0, (7, 1) = .0, (7, 2) = .0}, datatype = float[8], order = C_order); yout := Vector(6, {(1) = .0, (2) = .0, (3) = .0, (4) = .0, (5) = .0, (6) = .0}, datatype = float[8]); evalhf(`dsolve/numeric/hermite`(8, 6, X, Y, YP, outpoint, var(yout), var(L), var(V))) else if _EnvInFsolve = true then Digits := _EnvDSNumericSaveDigits end if; V := array( 1 .. 6, [( 1 ) = (7), ( 2 ) = (0), ( 3 ) = (false), ( 4 ) = (false), ( 5 ) = (false), ( 6 ) = (false)  ] ); L := Matrix(7, 2, {(1, 1) = 0., (1, 2) = 0., (2, 1) = 0., (2, 2) = 0., (3, 1) = 0., (3, 2) = 0., (4, 1) = 0., (4, 2) = 0., (5, 1) = 0., (5, 2) = 0., (6, 1) = 0., (6, 2) = 0., (7, 1) = 0., (7, 2) = 0.}, order = C_order); yout := Vector(6, {(1) = 0., (2) = 0., (3) = 0., (4) = 0., (5) = 0., (6) = 0.}); `dsolve/numeric/hermite`(8, 6, X, Y, YP, outpoint, yout, L, V) end if; [outpoint, seq(yout[i], i = 1 .. 6)] end proc, (2) = Array(0..0, {}), (3) = [eta, f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)], (4) = 0}); solnproc := data[1]; if not type(outpoint, 'numeric') then if outpoint = "solnprocedure" then return eval(solnproc) elif member(outpoint, ["start", "left", "right", "errorproc", "rawdata", "order", "error"]) then return solnproc(x_bvp) elif outpoint = "sysvars" then return data[3] elif procname <> unknown then return ('procname')(x_bvp) else _ndsol := pointto(data[2][0]); return ('_ndsol')(x_bvp) end if end if; try res := solnproc(outpoint); [eta = res[1], seq('[f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)]'[i] = res[i+1], i = 1 .. 6)] catch: error  end try end proc

 

proc (x_bvp) local res, data, solnproc, _ndsol, outpoint, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; _EnvDSNumericSaveDigits := Digits; Digits := 15; if _EnvInFsolve = true then outpoint := evalf[_EnvDSNumericSaveDigits](x_bvp) else outpoint := evalf(x_bvp) end if; data := Array(1..4, {(1) = proc (outpoint) local X, Y, YP, yout, errproc, L, V, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; X := Vector(8, {(1) = .0, (2) = .14285714285714277, (3) = .28571428571428553, (4) = .4285714285714284, (5) = .5714285714285713, (6) = .7142857142857142, (7) = .857142857142857, (8) = 1.0}, datatype = float[8], order = C_order); Y := Matrix(8, 6, {(1, 1) = .0, (1, 2) = .0, (1, 3) = .0, (1, 4) = .0, (1, 5) = 1.0, (1, 6) = .0, (2, 1) = .0, (2, 2) = .0, (2, 3) = .0, (2, 4) = .0, (2, 5) = 1.0, (2, 6) = .0, (3, 1) = .0, (3, 2) = .0, (3, 3) = .0, (3, 4) = .0, (3, 5) = 1.0, (3, 6) = .0, (4, 1) = .0, (4, 2) = .0, (4, 3) = .0, (4, 4) = .0, (4, 5) = 1.0, (4, 6) = .0, (5, 1) = .0, (5, 2) = .0, (5, 3) = .0, (5, 4) = .0, (5, 5) = 1.0, (5, 6) = .0, (6, 1) = .0, (6, 2) = .0, (6, 3) = .0, (6, 4) = .0, (6, 5) = 1.0, (6, 6) = .0, (7, 1) = .0, (7, 2) = .0, (7, 3) = .0, (7, 4) = .0, (7, 5) = 1.0, (7, 6) = .0, (8, 1) = .0, (8, 2) = .0, (8, 3) = .0, (8, 4) = .0, (8, 5) = 1.0, (8, 6) = .0}, datatype = float[8], order = C_order); YP := Matrix(8, 6, {(1, 1) = .0, (1, 2) = .0, (1, 3) = .0, (1, 4) = .0, (1, 5) = .0, (1, 6) = .0, (2, 1) = .0, (2, 2) = .0, (2, 3) = .0, (2, 4) = .0, (2, 5) = .0, (2, 6) = .0, (3, 1) = .0, (3, 2) = .0, (3, 3) = .0, (3, 4) = .0, (3, 5) = .0, (3, 6) = .0, (4, 1) = .0, (4, 2) = .0, (4, 3) = .0, (4, 4) = .0, (4, 5) = .0, (4, 6) = .0, (5, 1) = .0, (5, 2) = .0, (5, 3) = .0, (5, 4) = .0, (5, 5) = .0, (5, 6) = .0, (6, 1) = .0, (6, 2) = .0, (6, 3) = .0, (6, 4) = .0, (6, 5) = .0, (6, 6) = .0, (7, 1) = .0, (7, 2) = .0, (7, 3) = .0, (7, 4) = .0, (7, 5) = .0, (7, 6) = .0, (8, 1) = .0, (8, 2) = .0, (8, 3) = .0, (8, 4) = .0, (8, 5) = .0, (8, 6) = .0}, datatype = float[8], order = C_order); errproc := proc (x_bvp) local outpoint, X, Y, yout, L, V, i; option `Copyright (c) 2000 by Waterloo Maple Inc. All rights reserved.`; Digits := 15; outpoint := evalf(x_bvp); X := Vector(8, {(1) = .0, (2) = .14285714285714277, (3) = .28571428571428553, (4) = .4285714285714284, (5) = .5714285714285713, (6) = .7142857142857142, (7) = .857142857142857, (8) = 1.0}, datatype = float[8], order = C_order); Y := Matrix(8, 6, {(1, 1) = .0, (1, 2) = .0, (1, 3) = .0, (1, 4) = .0, (1, 5) = .0, (1, 6) = .0, (2, 1) = .0, (2, 2) = .0, (2, 3) = .0, (2, 4) = .0, (2, 5) = .0, (2, 6) = .0, (3, 1) = .0, (3, 2) = .0, (3, 3) = .0, (3, 4) = .0, (3, 5) = .0, (3, 6) = .0, (4, 1) = .0, (4, 2) = .0, (4, 3) = .0, (4, 4) = .0, (4, 5) = .0, (4, 6) = .0, (5, 1) = .0, (5, 2) = .0, (5, 3) = .0, (5, 4) = .0, (5, 5) = .0, (5, 6) = .0, (6, 1) = .0, (6, 2) = .0, (6, 3) = .0, (6, 4) = .0, (6, 5) = .0, (6, 6) = .0, (7, 1) = .0, (7, 2) = .0, (7, 3) = .0, (7, 4) = .0, (7, 5) = .0, (7, 6) = .0, (8, 1) = .0, (8, 2) = .0, (8, 3) = .0, (8, 4) = .0, (8, 5) = .0, (8, 6) = .0}, datatype = float[8], order = C_order); if not type(outpoint, 'numeric') then if outpoint = "start" or outpoint = "left" then return X[1] elif outpoint = "right" then return X[8] elif outpoint = "order" then return 2 elif outpoint = "error" then return HFloat(-0.0) elif outpoint = "errorproc" then error "this is already the error procedure" elif outpoint = "rawdata" then return [6, 8, [f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)], X, Y] else return ('procname')(x_bvp) end if end if; if outpoint < X[1] or X[8] < outpoint then error "solution is only defined in the range %1..%2", X[1], X[8] end if; V := array([1 = 4, 2 = 0]); if Digits <= trunc(evalhf(Digits)) then L := Vector(4, 'datatype' = 'float'[8]); yout := Vector(6, 'datatype' = 'float'[8]); evalhf(`dsolve/numeric/lagrange`(8, 6, X, Y, outpoint, var(yout), var(L), var(V))) else L := Vector(4, 'datatype' = 'sfloat'); yout := Vector(6, 'datatype' = 'sfloat'); `dsolve/numeric/lagrange`(8, 6, X, Y, outpoint, yout, L, V) end if; [eta = outpoint, seq('[f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)]'[i] = yout[i], i = 1 .. 6)] end proc; if not type(outpoint, 'numeric') then if outpoint = "start" or outpoint = "left" then return X[1] elif outpoint = "method" then return "bvp" elif outpoint = "right" then return X[8] elif outpoint = "order" then return 2 elif outpoint = "error" then return HFloat(-0.0) elif outpoint = "errorproc" then return eval(errproc) elif outpoint = "rawdata" then return [6, 8, "depnames", X, Y, YP] else error "non-numeric value" end if end if; if outpoint < X[1] or X[8] < outpoint then error "solution is only defined in the range %1..%2", X[1], X[8] end if; if Digits <= trunc(evalhf(Digits)) and (_EnvInFsolve <> true or _EnvDSNumericSaveDigits <= trunc(evalhf(Digits))) then V := array( 1 .. 6, [( 1 ) = (7), ( 2 ) = (0), ( 3 ) = (false), ( 4 ) = (false), ( 5 ) = (false), ( 6 ) = (false)  ] ); L := Matrix(7, 2, {(1, 1) = .0, (1, 2) = .0, (2, 1) = .0, (2, 2) = .0, (3, 1) = .0, (3, 2) = .0, (4, 1) = .0, (4, 2) = .0, (5, 1) = .0, (5, 2) = .0, (6, 1) = .0, (6, 2) = .0, (7, 1) = .0, (7, 2) = .0}, datatype = float[8], order = C_order); yout := Vector(6, {(1) = .0, (2) = .0, (3) = .0, (4) = .0, (5) = .0, (6) = .0}, datatype = float[8]); evalhf(`dsolve/numeric/hermite`(8, 6, X, Y, YP, outpoint, var(yout), var(L), var(V))) else if _EnvInFsolve = true then Digits := _EnvDSNumericSaveDigits end if; V := array( 1 .. 6, [( 1 ) = (7), ( 2 ) = (0), ( 3 ) = (false), ( 4 ) = (false), ( 5 ) = (false), ( 6 ) = (false)  ] ); L := Matrix(7, 2, {(1, 1) = 0., (1, 2) = 0., (2, 1) = 0., (2, 2) = 0., (3, 1) = 0., (3, 2) = 0., (4, 1) = 0., (4, 2) = 0., (5, 1) = 0., (5, 2) = 0., (6, 1) = 0., (6, 2) = 0., (7, 1) = 0., (7, 2) = 0.}, order = C_order); yout := Vector(6, {(1) = 0., (2) = 0., (3) = 0., (4) = 0., (5) = 0., (6) = 0.}); `dsolve/numeric/hermite`(8, 6, X, Y, YP, outpoint, yout, L, V) end if; [outpoint, seq(yout[i], i = 1 .. 6)] end proc, (2) = Array(0..0, {}), (3) = [eta, f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)], (4) = 0}); solnproc := data[1]; if not type(outpoint, 'numeric') then if outpoint = "solnprocedure" then return eval(solnproc) elif member(outpoint, ["start", "left", "right", "errorproc", "rawdata", "order", "error"]) then return solnproc(x_bvp) elif outpoint = "sysvars" then return data[3] elif procname <> unknown then return ('procname')(x_bvp) else _ndsol := pointto(data[2][0]); return ('_ndsol')(x_bvp) end if end if; try res := solnproc(outpoint); [eta = res[1], seq('[f(eta), diff(f(eta), eta), diff(diff(f(eta), eta), eta), diff(diff(diff(f(eta), eta), eta), eta), theta(eta), diff(theta(eta), eta)]'[i] = res[i+1], i = 1 .. 6)] catch: error  end try end proc

 

 

 

 

 

Download odeProb.mw

over 2D-input and 2D-math formats

Is the attached what you expected (For some reason this wont display inline  maybe the animations make the file size too big?)

anim.mw

@AHSAN 

the attached, maybe?

restart

A := beta*(-60*lambda*(1+(1/2)*x^2)^4+(-63*k^3+99*k^2-99*k+63)*(1+(1/2)*x^2)^3-288*Q*(k^2-(11/8)*k+1)*(1+(1/2)*x^2)^2-486*Q^2*(k-1)*(1+(1/2)*x^2)-324*Q^3)/(20*(1+(1/2)*x^2)^7)-(3*(k*(1+(1/2)*x^2)+2*Q-1-(1/2)*x^2))/(2*(1+(1/2)*x^2)^3)

opts := [[beta = .1, Q = .5516, lambda = -0.396e-1, k = .1], [beta = .2, Q = .4903, lambda = -0.411e-1, k = .2], [beta = .3, Q = .4290, lambda = -0.419e-1, k = .3], [beta = .4, Q = .3677, lambda = -0.417e-1, k = .4]]

yVals := [-.1, -.12, -.14, -.16]; N := plot([seq(eval(A, opts[j]), j = 1 .. 4)], x = -2 .. 2, color = [black, red, blue, green], thickness = 2); B := plot(yVals, x = -1 .. -.7, color = [black, red, blue, green], thickness = 3); T := plots:-textplot([seq([-1.1, yVals[j], typeset(k = rhs(opts[j, 3]), ", ", Q = rhs(opts[j, 2]))], j = 1 .. 4)], align = left)

plots:-display(N, B, T, size = [900, 700], gridlines = false)

 

 


 

Download legend2.mw

@user-JG 

as in

ifactor(725);
ifactor(1125);
ifactor(2048);

``(5)^2*``(29)

 

``(3)^2*``(5)^3

 

``(2)^11

(1)

 


 

Download fac.mw

the randomize() command with no arguments. This resets the "seed" for Maple's random number generators, based on the system clock, so the seed will be different each time your worksheet is executed.

Without using "randomize()", you will get the same random numbers every time you execute a worksheet as you have observed. This behaviour is by design (and very useful for "debug" purposes!)

which will return the "next" line as a string in the input file.

What you want to do with this string is up to you - parse() maybe?

Since you now have a string up to a linebreak, how you "signify" the linebreaks is also up to you

Solution to your problem depends very much on how you want "linebreaks" to be signified in the output - eg a new row in a matrix?, a new entry in a table? whatever? Ypu are going to have to get a lot more specific on

  • what to do with each line
  • what to do with the linebreak

working out your problem!

Is this a 1D math input versus 2D math input issue?

See the attached for a couple of solutions in both input formats - which of these is difficult to understand?

I have taken into account the remark by CarlLove that the VectorField() command is a bit "redundant" - so, in the attached, there are execution groups which use it, and execution groups which don't - doesn't affect the answer.

Which of these don't you understand?

restart; with(Student:-VectorCalculus)

g := x^2+y^2+z^2
NULL

x^2+y^2+z^2

(1)

G2D := unapply('evalVF(VectorField(Gradient(z)), p)', [z, p]); G2D(g, `<,>`(1, 1, 1))

G2D := proc (z, p) options operator, arrow; evalVF(VectorField(Gradient(z)), p) end proc

 

Vector[column](%id = 18446744074370475182)

(2)

G2D := unapply('evalVF(Gradient(z), p)', [z, p]); G2D(g, `<,>`(1, 1, 1))

G2D := proc (z, p) options operator, arrow; evalVF(Gradient(z), p) end proc

 

Vector[column](%id = 18446744074370472406)

(3)

  G1D:= unapply
        ( 'evalVF
           ( VectorField
             (  Gradient(z)
             ),
             p
           )',
           [z,p]
         );
  G1D( g, <1,1,1> );

G1D := proc (z, p) options operator, arrow; evalVF(VectorField(Gradient(z)), p) end proc

 

Vector[column](%id = 18446744074370467230)

(4)

  G1D:= unapply
        ( 'evalVF
          (  Gradient(z),
             p
          )',
          [z,p]
        );
  G1D( g, <1,1,1> );

G1D := proc (z, p) options operator, arrow; evalVF(Gradient(z), p) end proc

 

Vector[column](%id = 18446744074370461814)

(5)

 

NULL

Download SVC2.mw

  1. You had a couple of typos where you had written, for example,  diff*(v(z),z), rather than diff(v(z),z)
  2. Why are you using pdsolve() to solve a system of ODEs??? Changed this to dsolve().

Maple now reports no errors, although it does not povide a solution. You may have to consider generating a 'numerical' solutiion.

See the attache
 

restart

NULL

D21 := 3; D31 := 14; D41 := 12; D51 := 23; D61 := 12; D71 := 14; D81 := 14

D32 := 5; D42 := 12; D52 := 23; D62 := 123; D72 := 16; D82 := 24

D43 := 8; D53 := 15; D63 := 34; D73 := 18; D83 := 79

D54 := 9; D64 := 81; D74 := 20; D84 := 17

D65 := 10; D75 := 2; D85 := 48; D86 := 6; D87 := 13; D88 := 19

D11 := 1; D22 := 2; D33 := 4; D44 := 5; D55 := 7; D66 := 10; D77 := 21

D76 := 14

D12 := D21; D13 := D31; D14 := D41; D15 := D51; D16 := D61; D17 := D71; D18 := D81

D23 := D32; D24 := D42; D25 := D52; D26 := D62; D27 := D72; D28 := D82

D34 := D43; D35 := D53; D36 := D63; D37 := D73; D38 := D83

D45 := D54; D46 := D64; D47 := D74; D48 := D84

D56 := D65; D57 := D75; D58 := D85

D67 := D76; D68 := D86

D78 := D87

xs := 0; ys := 0

NULL

e1 := diff(w(z), z)+1/2*((diff(u(z), z))^2+(diff(v(z), z))^2)+(diff(f(z), z))*(ys*((diff(u(z), z))*cos(f(z))+(diff(v(z), z))*sin(f(z)))+xs*((diff(u(z), z))*sin(f(z))-(diff(v(z), z))*cos(f(z))))

diff(w(z), z)+(1/2)*(diff(u(z), z))^2+(1/2)*(diff(v(z), z))^2

(1)

NULL

e2 := diff(tx(z), z)-(diff(f(z), z))*((diff(u(z), z))*cos(f(z))+(diff(v(z), z))*sin(f(z)))

diff(tx(z), z)-(diff(f(z), z))*((diff(u(z), z))*cos(f(z))+(diff(v(z), z))*sin(f(z)))

(2)

NULL

e3 := diff(ty(z), z)+(diff(f(z), z))*((diff(u(z), z))*sin(f(z))-(diff(v(z), z))*cos(f(z)))

diff(ty(z), z)+(diff(f(z), z))*((diff(u(z), z))*sin(f(z))-(diff(v(z), z))*cos(f(z)))

(3)

NULL

e4 := dr*(diff(f(z), z, z))

dr*(diff(diff(f(z), z), z))

(4)

NULL

e7 := -ty(z)+(diff(v(z), z))*sin(f(z))+(diff(u(z), z))*cos(f(z))

-ty(z)+(diff(v(z), z))*sin(f(z))+(diff(u(z), z))*cos(f(z))

(5)

NULL

e6 := tx(z)-(diff(u(z), z))*sin(f(z))+(diff(v(z), z))*cos(f(x))

tx(z)-(diff(u(z), z))*sin(f(z))+(diff(v(z), z))*cos(f(x))

(6)

NULL

e8 := diff(f(z), z)

diff(f(z), z)

(7)

NULL

e5 := (1/2)*(diff(f(z), z))^2

(1/2)*(diff(f(z), z))^2

(8)

``

DD := matrix([[D11, D12, D13, D14, D15, D16, D17, D18], [D21, D22, D23, D24, D25, D26, D27, D28], [D31, D32, D33, D34, D35, D36, D37, D38], [D41, D42, D43, D44, D45, D46, D47, D48], [D51, D52, D53, D54, D55, D56, D57, D58], [D61, D62, D63, D64, D65, D66, D67, D68], [D71, D72, D73, D74, D75, D76, D77, D78], [D81, D82, D83, D84, D85, D86, D87, D88]])

array( 1 .. 8, 1 .. 8, [( 3, 6 ) = (34), ( 2, 8 ) = (24), ( 7, 5 ) = (2), ( 4, 4 ) = (5), ( 6, 7 ) = (14), ( 6, 6 ) = (10), ( 5, 8 ) = (48), ( 8, 3 ) = (79), ( 7, 4 ) = (20), ( 6, 1 ) = (12), ( 8, 8 ) = (19), ( 1, 5 ) = (23), ( 3, 4 ) = (8), ( 1, 2 ) = (3), ( 6, 5 ) = (10), ( 1, 3 ) = (14), ( 5, 7 ) = (2), ( 7, 6 ) = (14), ( 6, 8 ) = (6), ( 1, 4 ) = (12), ( 3, 5 ) = (15), ( 2, 7 ) = (16), ( 5, 1 ) = (23), ( 8, 2 ) = (24), ( 1, 1 ) = (1), ( 1, 6 ) = (12), ( 3, 2 ) = (5), ( 5, 3 ) = (15), ( 2, 1 ) = (3), ( 6, 4 ) = (81), ( 8, 7 ) = (13), ( 4, 6 ) = (81), ( 7, 1 ) = (14), ( 2, 3 ) = (5), ( 2, 5 ) = (23), ( 1, 7 ) = (14), ( 4, 5 ) = (9), ( 3, 3 ) = (4), ( 4, 1 ) = (12), ( 7, 7 ) = (21), ( 5, 6 ) = (10), ( 7, 2 ) = (16), ( 2, 4 ) = (12), ( 4, 7 ) = (20), ( 8, 6 ) = (6), ( 2, 2 ) = (2), ( 3, 8 ) = (79), ( 6, 3 ) = (34), ( 6, 2 ) = (123), ( 5, 4 ) = (9), ( 8, 1 ) = (14), ( 7, 3 ) = (18), ( 8, 4 ) = (17), ( 2, 6 ) = (123), ( 1, 8 ) = (14), ( 4, 3 ) = (8), ( 4, 2 ) = (12), ( 8, 5 ) = (48), ( 3, 7 ) = (18), ( 4, 8 ) = (17), ( 3, 1 ) = (14), ( 5, 5 ) = (7), ( 5, 2 ) = (23), ( 7, 8 ) = (13)  ] )

(9)

F := DD.e

e := Vector([e1, e2, e3, e4, e5, e6, e7, e8])

Vector[column](%id = 18446744074565103238)

(10)

N := F[1]; Mx := F[2]; My := F[3]; Mw := F[4]; gT := F[5]; Vy := F[6]; Vx := F[7]; Ts := F[8]

NULL

diff(N*(diff(u(z), z)), z)

eq1 := diff(N*(`*`(diff(u(z), z)+(diff(f(z), z))*(ys*cos(f(z))+xs*sin(f(z))))), z)+diff((Vx-Mx*(diff(f(z), z)))*cos(f(z))+(My*(diff(f(z), z))-Vy)*sin(f(z)), z) = 0

eq2 := diff(N*(`*`(diff(v(z), z)+(diff(f(z), z))*(ys*sin(f(z))-xs*sin(f(z))))), z)+diff((Vx-Mx*(diff(f(z), z)))*sin(f(z))+(Vy-My*(diff(f(z), z)))*cos(f(z)), z) = 0

eq3 := diff(N, z) = 0

eq4 := diff(Mx, z)-Vy = 0

eq5 := diff(My, z)+Vx = 0

eq6 := diff(N*(cos(f(z))*((diff(u(z), z))*ys-(diff(v(z), z))*xs)+sin(f(z))*((diff(v(z), z))*ys+(diff(u(z), z))*xs)), z)-N*(diff(f(z), z))*(cos(f(z))*((diff(v(z), z))*ys+(diff(u(z), z))*xs)+sin(f(z))*((diff(v(z), z))*xs-(diff(u(z), z))*ys))-(diff(cos(f(z))*(Mx*(diff(u(z), z))+My*(diff(v(z), z)))+sin(f(z))*(Mx*(diff(v(z), z))-My*(diff(u(z), z))), z))+cos(f(z))*(Mx*(diff(v(z), z))*(diff(f(z), z))-My*(diff(u(z), z))*(diff(f(z), z))+Vy*(diff(u(z), z))-Vx*(diff(v(z), z)))+sin(f(z))*(-Mx*(diff(u(z), z))*(diff(f(z), z))-My*(diff(v(z), z))*(diff(f(z), z))+Vy*(diff(v(z), z))+Vx*(diff(u(z), z)))-dr*(diff(Mw, z, z))+diff(gT*(diff(f(z), z)), z)+diff(Ts, z) = 0

NULL

NULL

NULL

NULL

NULL

NULL

ode := eq1, eq2, eq3, eq4, eq5, eq6

 

sol := dsolve([ode], {f(z), tx(z), ty(z), u(z), v(z), w(z)})

(11)

sol := dsolve({eq1, eq2, eq3, eq4, eq5, eq6, f(0) = f1, f(L) = f2, tx(0) = tx1, tx(L) = tx2, ty(0) = ty1, ty(L) = ty2, u(0) = u1, u(L) = u2, v(0) = v1, v(L) = v2, w(0) = w1, w(L) = w2}, {f(z), tx(z), ty(z), u(z), v(z), w(z)})

(12)

``

NULL


 

Download diffProb.mw

d

Just cheat using the 'view' option, whihc I have (sort of) generalized in the attached

  restart;
  xlim:=10:
  ylim:=10:
  expr:=4*x^2+9*y^2:
  plot3d( expr,
          x=-xlim..xlim,
          y=-ylim..ylim,
          view=[ -xlim..xlim,
                 -ylim..ylim,
                  0..min( [ eval(expr,[x=0, y=ylim]),
                            eval(expr,[y=0, x=xlim])
                          ]
                        )
               ],
         style=surface,
         colorscheme=["zcoloring", z->z]
       );

 

 


 

Download paraPlot.mw

 

 

The a, b, and c values in the above coordinate transformations can be given using the coordinate specification as a function, e.g., conical(a,b) or ellcylindrical(2). The values a, b, and c if necessary, should be specified.  If not specified, the default values used are a = 1, b = 1/2, and c = 1/3.

 

The toolbar on a Maple help page has navigation icons at the top left. The second thrid and fourth buttons (from the left) will allow you to navigate backwards, forwards, or upwards to the parent topic.

Are you suggesting that these buttons don't work?

Also not clear what your problem is in creating a "function". See the attached for a "toy" example

  restart;
  alias(SVC=Student[VectorCalculus]):
#
# Define a function which takes an expression and a point.
# It creates a vector field from the gradient of the
# supplied expression and then evaluates the vector field
# at the supplied point
#
  f:= (z,p)-> SVC:-evalVF
              ( SVC:-VectorField
                (  SVC:-Gradient(z)
                ),
                p
              );

 f( exp(x*y^3*z^2), <1,1,1> );

f := proc (z, p) options operator, arrow; SVC:-evalVF(SVC:-VectorField(SVC:-Gradient(z)), p) end proc

 

Vector[column](%id = 18446744074397177606)

(1)

 

 


 

Download SVC.mw

 

1 2 3 4 5 6 7 Last Page 1 of 140