% Taking inputs % Length of link-1 (L1) and link-2 (L2) in meter L1=1; L2=0.5;   % Angle traced by link-1 in degrees (theta_1) theta_1=linspace(0,90,11);   % Angle traced by link-2 in degrees (theta_2) theta_2=linspace(0,90,21);   % Initialising base coordinates of link-1 x_0=0; y_0=0;   ct=1;  …

in the above equation g= gravity in m/s^2 L= length of the pendulum in meter. m= mass of the ball in kg. b= damping coeffecient.   let Θ=Θ1 dΘ1/dt=Θ2 dΘ^2/dt^2=dΘ/dt   Input: L=1 metre m= 1kg b=0.05 g= 9.81m/s^2 time span=0-20 sec angular displacement=0 angular velocity=0 angular…

CODE : Function code : function [f] = stalagmite(input_vector)x = input_vector(1);y = input_vector(2); f1_x = (sin((5.1*pi*x)+0.5)).^6; f1_y = (sin((5.1*pi*y)+0.5)).^6; f2_x = exp(((-4*log(2))*(x-0.0667).^2)/0.64); f2_y = exp(((-4*log(2))*(y-0.0667).^2)/0.64); [f] = -(f1_x.*f1_y.*f2_x.*f2_y); end Main code : clear all…

Where, Cp=specific heat R =molar gas constant T=temperature S=Entropy H=Enthalpy a=coefficients   BODY OF THE CONTENT: Code for calling function Specific_heat(): function Cp = Specific_heat(a1,a2,a3,a4,a5,a8,a9,a10,a11,a12,T,R,global_midtemp) %calculating specific_heat for higher and lower temperature respectively if(T>=…