clear all close all clc %inputs l1=1; l2=0.5; th1=linspace(0,90,10); th2=linspace(0,90,10); ct=1; for i=1:length(th1) TH1=th1(i); for j=1:length(th2) TH2=th2(j); x0=0; y0=0; x1=l1*cosd(TH1); y1=l1*sind(TH1); x2=x1+l2*cosd(TH2); y2=y1+l2*sind(TH2); plot([x0 x1],[y0 y1],[x1 x2],[y1 y2],\'linewidth\',3) axis([-0.1…

clear all close all clc b=0.05; g=9.81; l=1; m=1; %initial_con th_0=[0;3]; %time plot tspan=linspace(0,20,500); %solve ode [t,results]= ode45(@(time,theta) ode_func(t,th,b,g,l,m),tspan,th_0); subplot(5,4,[13,14,17,18]); plot(t,results(:,1)) hold on; plot(t,results(:,2)) hold off subplot(5,4,[15,16,19,20]); plot(t,results(:,2));…

The objective for this program is to write a function a function that extracts the 14 co-efficients and calculates the enthalpy, entropy and specific heats for all the species in the data file. NASA  came up with polynomials that can be used to evaluate thermodynamic properties such as Cp, H and S. They have…

%column vector a=[1 2 3 4 5] %row vector a=[1 2 3 4 5]        

clear all close all clc %inputs g=1.4; t3=2300; %state variables p1=101325; t1=500; %egnine geometry bore=0.1; strk=0.1; cnrod=0.15; cr=12; %swept volume and clearance volume vswept=(pi/4)*bore^2*strk; vclear=vswept/(cr-1); v1=vswept+vclear; v2=vclear; %state 2 p2=p1*cr^g; t2=p2*v2*t1/(p1*v1); cons=p1*v1^g; vcomp=engine_kin(bore,strk,cnrod,cr,180,0);…

clear all close all clc %loading data file cp_data= load(\'data1\'); temp= cp_data(:,1); cp=cp_data(:,2); %cp= a*T^2 + b*T +c and a*T^3 + b*T^2 +c*T+d coeff = polyfit(temp,cp,2); predcp= polyval(coeff, temp); coeff1 = polyfit(temp,cp,3); predcp1=polyval(coeff1,temp); %plotting subplot(2,1,[1]); plot(temp, cp,\'linewidth\',3);…

In fluid mechanics the drag force exerted by the surrounding medium( fluid) plays an important role in motion study. The drag force is calculated by the formula, dragforce=0.5*(density of fluid)*(frontal area)*(velocity)^2*(drag coefficient) 1. Drag force vs velocity function [dragforce]=myf(a,v,rho,cd) dragforce=0.5*a*rho*v.*v.*cd;…

clear all close all clc x=linspace(0,0.6,150); y=linspace(0,0.6,150); [xx yy]= meshgrid(x,y); for i=1:length(xx) for j=1:length(yy) input(1)=xx(i,j); input(2)=yy(i,j); f(i,j)=stalagmit(input); end end surfc(xx,yy,f) shading interp [inputs, fval]=ga(@stalagmit,2) This program is about solving a global optimization problem.…