В Maple нахожу решение:

> restart:
> with(LinearAlgebra):
> L:=1:x1:=L/6: x2:=(L/2): x3:=5*(L/6):m:=1: E:=1:J:=1:
> p1 := 7.433266872: p2:= 47.43856088: p3 := 81.34482671:
> k:=((m*p^2)/(E*J))^(1/4):
Функции:
> z1 := (k,x)->((1-1/2*(k*x)^2+1/24*(k*x)^4)+(1+1/2*(k*x)^2+1/24*(k*x)^4))/2:
> z2 := (k,x)->(1*(k*x)-1/6*(k*x)^3+1/120*(k*x)^5+1*(k*x)+1/6*(k*x)^3+1/120*(k*x)^5)/(2*k):
> z3 :=  (k,x)->(-(1-1/2*(k*x)^2+1/24*(k*x)^4)+(1+1/2*(k*x)^2+1/24*(k*x)^4))/(2*k^2):
> z4 :=  (k,x)->(-((1*(k*x)-1/6*(k*x)^3+1/120*(k*x)^5+O((k*x)^6)))+(1*(k*x)+1/6*(k*x)^3+1/120*(k*x)^5+O((k*x)^6)))/(2*k^3):

> z11:=z1(k,0): z21:=z2(k,0): z31:=z3(k,0): z41:=z4(k,0):
> z12:=z1(k,x1): z22:=z2(k,x1): z32:=z3(k,x1): z42:=z4(k,x1):
> z13:=z1(k,x2): z23:=z2(k,x2): z33:=z3(k,x2): z43:=z4(k,x2):
> z14:=z1(k,x3): z24:=z2(k,x3): z34:=z3(k,x3): z44:=z4(k,x3):
> z15:=z1(k,L): z25:=z2(k,L): z35:=z3(k,L): z45:=z4(k,L):
> with(linalg):
Warning, the previous binding of the name GramSchmidt has been removed and it now has an assigned value

Warning, the protected names norm and trace have been redefined and unprotected

Матрица из коээфициентов системы:
> A := matrix( [ [z11,z21,z31,z41,0,0,0,0,0,0,0,0,0,0,0,0], [k^4*z31, k^4*z41, z11, z21,0,0,0,0,0,0,0,0,0,0,0,0], [z12,z22,z32,z42,-z12,-z22,-z32,-z42,0,0,0,0,0,0,0,0], [k^4*z42,z12,z22,z32,-k^4*z42,-z12,-z22,-z32,0,0,0,0,0,0,0,0], [E*J*k^4*z32,E*J*k^4*z42, E*J*z12,E*J*z22, -E*J*k^4*z32, -E*J*k^4*z42, -E*J*z12, -E*J*z22, 0,0,0,0,0,0,0,0], [E*J*k^4*z22+m1*p^2*z12, E*J*k^4*z32+m1*p^2*z22, E*J*k^4*z42+m1*p^2*z32, E*J*z12+m1*p^2*z42, -E*J*k^4*z22, -E*J*k^4*z32, -E*J*k^4*z42, -E*J*z12,0,0,0,0,0,0,0,0], [0,0,0,0, z13, z23, z33, z43,-z13, -z23, -z33, -z43,0,0,0,0 ], [0,0,0,0,k^4*z43,z13,z23,z33,-k^4*z43,-z13,-z23,-z33, 0,0,0,0 ], [0,0,0,0,E*J*k^4*z33,E*J*k^4*z43, E*J*z13,E*J*z23, -E*J*k^4*z33, -E*J*k^4*z43, -E*J*z13, -E*J*z23,0,0,0,0], [0,0,0,0,E*J*k^4*z23+m1*p^2*z13, E*J*k^4*z33+m1*p^2*z23, E*J*k^4*z43+m1*p^2*z33, E*J*z13+m1*p^2*z43, -E*J*k^4*z23, -E*J*k^4*z33, -E*J*k^4*z43, -E*J*z13, 0,0,0,0], [0,0,0,0,0,0,0,0,0,0,0,0,z15, z25, z35, z45], [0,0,0,0,0,0,0,0,0,0,0,0,k^4*z35, k^4*z45, z15, z25], [0,0,0,0,0,0,0,0,z14, z24,z34,z44,-z14, -z24,-z34,-z44], [0,0,0,0,0,0,0,0,k^4*z44,z14,z24,z34,-k^4*z44,-z14,-z24,-z34], [0,0,0,0,0,0,0,0,E*J*k^4*z34,E*J*k^4*z44, E*J*z14,E*J*z24, -E*J*k^4*z34, -E*J*k^4*z44, -E*J*z14, -E*J*z24], [0,0,0,0,0,0,0,0,E*J*k^4*z24+m1*p^2*z14, E*J*k^4*z34+m1*p^2*z24, E*J*k^4*z44+m1*p^2*z34, E*J*z14+m1*p^2*z44, -E*J*k^4*z24, -E*J*k^4*z34, -E*J*k^4*z44, -E*J*z14]]):
Нахождение масс :
> p:=p1: 
> Dp1:=det(A):
> M:=solve(Dp1,m1);
      M := .3333332271, 1.291639979, 33.00472012
>

Нужно оформить это решение в виде процедуры.
Делаю так:

Процедура нахождения масс 
> Mass:=proc(L,x1,x2,x3,m,E,J,p1,p2,p3) local k,z1,z2,z3,z4,z11,z21,z31,z41,z12,z22,z32,z42,z13,z23,z33,z43,z14,z24,z34,z44,z15,z25,z35,z45,A,p,Dp1,Dp2,Dp3; global M;
>  k:=((m*p^2)/(E*J))^(1/4); 
Функции :
> z1 := (k,x)->((1-1/2*(k*x)^2+1/24*(k*x)^4)+(1+1/2*(k*x)^2+1/24*(k*x)^4))/2;
> z2 := (k,x)->(1*(k*x)-1/6*(k*x)^3+1/120*(k*x)^5+1*(k*x)+1/6*(k*x)^3+1/120*(k*x)^5)/(2*k);
> z3 :=  (k,x)->(-(1-1/2*(k*x)^2+1/24*(k*x)^4)+(1+1/2*(k*x)^2+1/24*(k*x)^4))/(2*k^2);
> z4 :=  (k,x)->(-((1*(k*x)-1/6*(k*x)^3+1/120*(k*x)^5+O((k*x)^6)))+(1*(k*x)+1/6*(k*x)^3+1/120*(k*x)^5+O((k*x)^6)))/(2*k^3);

> z11:=z1(k,0); z21:=z2(k,0); z31:=z3(k,0); z41:=z4(k,0);
> z12:=z1(k,x1); z22:=z2(k,x1); z32:=z3(k,x1); z42:=z4(k,x1);
> z13:=z1(k,x2); z23:=z2(k,x2); z33:=z3(k,x2); z43:=z4(k,x2);
> z14:=z1(k,x3); z24:=z2(k,x3); z34:=z3(k,x3); z44:=z4(k,x3);
> z15:=z1(k,L); z25:=z2(k,L); z35:=z3(k,L); z45:=z4(k,L);
> with(linalg):
Матрица из коэффициентов системы:
> A := matrix( [ [z11,z21,z31,z41,0,0,0,0,0,0,0,0,0,0,0,0], [k^4*z31, k^4*z41, z11, z21,0,0,0,0,0,0,0,0,0,0,0,0], [z12,z22,z32,z42,-z12,-z22,-z32,-z42,0,0,0,0,0,0,0,0], [k^4*z42,z12,z22,z32,-k^4*z42,-z12,-z22,-z32,0,0,0,0,0,0,0,0], [E*J*k^4*z32,E*J*k^4*z42, E*J*z12,E*J*z22, -E*J*k^4*z32, -E*J*k^4*z42, -E*J*z12, -E*J*z22, 0,0,0,0,0,0,0,0], [E*J*k^4*z22+m1*p^2*z12, E*J*k^4*z32+m1*p^2*z22, E*J*k^4*z42+m1*p^2*z32, E*J*z12+m1*p^2*z42, -E*J*k^4*z22, -E*J*k^4*z32, -E*J*k^4*z42, -E*J*z12,0,0,0,0,0,0,0,0], [0,0,0,0, z13, z23, z33, z43,-z13, -z23, -z33, -z43,0,0,0,0 ], [0,0,0,0,k^4*z43,z13,z23,z33,-k^4*z43,-z13,-z23,-z33, 0,0,0,0 ], [0,0,0,0,E*J*k^4*z33,E*J*k^4*z43, E*J*z13,E*J*z23, -E*J*k^4*z33, -E*J*k^4*z43, -E*J*z13, -E*J*z23,0,0,0,0], [0,0,0,0,E*J*k^4*z23+m1*p^2*z13, E*J*k^4*z33+m1*p^2*z23, E*J*k^4*z43+m1*p^2*z33, E*J*z13+m1*p^2*z43, -E*J*k^4*z23, -E*J*k^4*z33, -E*J*k^4*z43, -E*J*z13, 0,0,0,0], [0,0,0,0,0,0,0,0,0,0,0,0,z15, z25, z35, z45], [0,0,0,0,0,0,0,0,0,0,0,0,k^4*z35, k^4*z45, z15, z25], [0,0,0,0,0,0,0,0,z14, z24,z34,z44,-z14, -z24,-z34,-z44], [0,0,0,0,0,0,0,0,k^4*z44,z14,z24,z34,-k^4*z44,-z14,-z24,-z34], [0,0,0,0,0,0,0,0,E*J*k^4*z34,E*J*k^4*z44, E*J*z14,E*J*z24, -E*J*k^4*z34, -E*J*k^4*z44, -E*J*z14, -E*J*z24], [0,0,0,0,0,0,0,0,E*J*k^4*z24+m1*p^2*z14, E*J*k^4*z34+m1*p^2*z24, E*J*k^4*z44+m1*p^2*z34, E*J*z14+m1*p^2*z44, -E*J*k^4*z24, -E*J*k^4*z34, -E*J*k^4*z44, -E*J*z14]]):
Нахождение масс :
> p:=p1; Dp1:=Det(A);
> M:=solve(Dp,m1);
> end proc:

Вызываю процедуру:

 Mass(1,1/6,1/2,5/6,1,1,1,7.433266872,47.43856088,81.34482671);

И ничего не происходит,т.е. решение не ищется, хотя когда делаешь тоже самое без процедуры решение находится. В чем ошибка? Помогите, пожалуйста!

Отличий в коде программы и процедуры нет, кроме соответствующих команд для задания процедуры.

Смотрите также