I hereby claim:
- I am jdelacruz26 on github.
- I am jdelacruz26 (https://keybase.io/jdelacruz26) on keybase.
- I have a public key ASD4PAQaTi3mUMvmmTkMppNlBK27aSDq_utxwV4rqKurdQo
To claim this, I am signing this object:
I hereby claim:
To claim this, I am signing this object:
sudo apt-get update
sudo apt-get install terminator
%/////////////////////////////////////////////////// | |
%// Función Secante /// | |
%// Desarrollada por Jorge De La Cruz /// | |
%/ /// | |
%////////////////////////////////////////////////// | |
function xr=secante(fun,x0,x1,tol,imax) | |
x=[x0 x1]; | |
f=eval(fun); | |
err=tol+1; | |
n=1; |
%Función Gauss-Seidel | |
%Desarrollada por: Jorge De La Cruz | |
%Para un sistema de la forma 'Ax=b' | |
%imax: Número máxima de iteraciones en caso de | |
%que el método no converja | |
%tol: tolerancia deseada | |
function x=gseidel(A,b,tol,imax) | |
n=max(size(A)); % Obtenemos de la matriz A |
function yp=derivadanolin(t,y) | |
global L m P k c | |
yp(1,1)=y(3); | |
yp(2,1)=y(4); | |
yp(3,1)=(P*L*sin(y(1)-y(2))-k*y(1)-c*y(3)-0.5*m*L^2*y(4)^2*sin(y(1)-y(2))... |
function yp=derivada(t,y) | |
global L m P k c | |
yp(1,1)=y(3); | |
yp(2,1)=y(4); | |
% yp(3,1)=(-(P*L-(3/2)*k)*y(2)-(k-P*L)*y(1))*(12/(25*m*L^2)); | |
% yp(4,1)=(-k*y(2)-0.5*m*L^2*yp(3))*3/(m*L^2); | |
yp(3,1)=(-(P*L-(3/2)*k)*y(2)-(k-P*L)*y(1)-c*y(3)-(3/2)*c*y(4))*(12/(25*m*L^2)); | |
yp(4,1)=(-k*y(2)-0.5*m*L^2*yp(3)-c*y(4))*3/(m*L^2); |
%%%///////////////////////////////////////////////// | |
% / Función Newton-Raphson /// | |
%/ Desarrollada por Jorge De La Cruz /// | |
%/ /// | |
%//////////////////////////////////////////////// | |
% La función 'fun' debe ser introducida como un string | |
% x: representa el valor inicial necesario para iniciar | |
% las iteraciones. | |
% tol: La tolerancia deseada. | |
% imax: máximo número de iteraciones, evita que el |
global L m P k c | |
k=100; | |
c=0.4; | |
m=0.3; | |
L=1; | |
P=99; | |
t=0:0.01:20; | |
y=[0.01 0 0 0]; |