From: Erick Macote Yparraguirre on
>> syms x L P1 P4 A1 A2 A3 A4 R E
>> N1=-9/2/L^3*((x-L/3)*(x-2*L/3)*(x-L));
>> N2=27/2/L^3*(x*(x-2*L/3)*(x-L));
>> N3=-27/2/L^3*(x*(x-L/3)*(x-L));
>> N4=9/2/L^3*((x-L/3)*(x-2*L/3)*x);
>> N = [N1,N2,N3,N4];
>> A=A1*N1+A2*N2+A3*N3+A4*N4;
>> P=P1*(1-x/L)+P4*x/L;
>> B=simplify(diff(N)) %matriz B
>> Nt=N.';
>> Bt=B.';
>> M=simplify(int(Nt*N*R*A,0,L)); %MATRIZ M
>> K=simplify(int(Bt*B*E*A,0,L)); %MATRIZ K
>> f=simplify(int(Nt*P,0,L)); %VECTOR F
>> det(K)
PREGUNTAS :
1)Como hago en matlab para verificar que la suma de los elementos de la matriz M es un numero positivo???
2)COMO ASIGNO VALORES a L P1 P4 A1 A2 A3 A4 R E , (conjuntos distintos de datos)SIENDO L LONGITUD Y A LA SECCION P la CARGA ASI COMO R EL VALOR DE LA DENSIdaD, y obtener las nuevas Matrices la ASIGNACION DE VALORES SON A ELECCION.????
3)
Evaluar las matrices [ M ], [ K ] y { f } obtenidas en el punto anterior para 3 (tres)
conjuntos distintos de datos, y con las matrices calculadas formar las otras matrices
de acuerdo con la tabla siguiente:

Conjunto de datos Fila / Columna Matriz Local Fila / Columna Matriz Global


1 3
1 2 4
3 7
4 8


1 1
2 2
2 3 5
4 6


1 4
3 2 5
3 6
4 7

4) esto es para el conjunto de datos 1 .Aqui no se como hacer,por ejemplo el elemeto 11 de la matriz M la cual se obtiene para el conjunto de datos, es decir por medio de la asignacion de valores, TIENE q estar ubicado en la nueva matriz en la posicion 33, el elemento 12 en la posicion 34, el elemento 13 en la posicion 37,el elemento 14 en la posicion 38, el elemento 21 en la posicion 43 el elemento 22 en la posision 44, el elemento 23 en la posicion 47,el elemento 24 en la posicion 48, del mismo modo para el 3 Y 4 , lo mismo para el conjunto de datos 2 y 3

NO SE COMO HACER ESO Y MOSTRAR LA MATRIZ FINAL CON LAS NUEVAS POSICIONES Y LOS VALORES CORRESPONDIENTES