%Rec 0.1 %KMWard clear; clc; L=10; EA=100; Amp=3; type=4; g=@(x) 3*x/L n=21 %Number of divisions along L h=L/(n-1); R=re
Posted: Tue Jul 12, 2022 8:21 am
23 24 25 26 7 23 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 Ma I1=1/EA*I1 Z=[10;11] d=L/EA B=[1,0,-1,0; d,1,0,-1] BC=[1,0,1,0] C=B(:,BC==1) S=C\z f=[0;0;0;0] f(BC==1)=s u0=f(2) NO=f(1) nT=51 beta=0.5; dX= L/(nT-1); Nold=NO; uold=u0; K XT=[0:dX:L]; pold-Amp Load FunctionCP (0, L, type); SA(1,:)=[0 pold Nold uold] for k=2:nT xold=xT(k-1) xnew=xT (k); Gold-Am-Eline-incolu etionCP.m