Laboratorio 2

install.packages(“MASS”) library(MASS)

eqn_2o<-function(A1, A2, y0, y1, x){ if(A1^2-4A2<0){ print(“Imaginarias”)
} else if(A1^2-4A2 == 0){ print(“Reales iguales”) #Resolviendo le eq cuadrática m1 <-(-A1+sqrt(A1^2-4A2))/2 m2 <-(-A1-sqrt(A1^2-4A2))/2 #Resolviendo el sistema de ecuaciones a1<-c(1,m1) a2<-c(1,m2x) A<-cbind(a1,a2) b1<-c(y0,y1) B=cbind(b1) C=ginv(A)%%B yx = C[1]m1^x+C[2]xm2x } else if (A1^2-4A2 > 0){ print(“Reales diferentes”) #Resolviendo le eq cuadrática m1 <-(-A1+sqrt(A1^2-4A2))/2 m2 <-(-A1-sqrt(A1^2-4A2))/2 #Resolviendo el sistema de ecuaciones a1<-c(1,m1) a2<-c(1,m2) A<-cbind(a1,a2) b1<-c(y0,y1) B=cbind(b1) C=ginv(A)%%B yx = C[1]m1^x+C[2]*m2^x } }

y_x<-eqn_2o(8,4,5,7,10) y_x