I know the databases I am using are correct….however, seems like when I try to attach a dbc file to my simulink model, I’ll get a warning that the dbc file MAY NOT be a valid J1939 database……
again, sometimes I get the message and simulink WILL NOT select the database during config……but other times, I CAN SELECT the same database and get no warnings…???????I know the databases I am using are correct….however, seems like when I try to attach a dbc file to my simulink model, I’ll get a warning that the dbc file MAY NOT be a valid J1939 database……
again, sometimes I get the message and simulink WILL NOT select the database during config……but other times, I CAN SELECT the same database and get no warnings…??????? I know the databases I am using are correct….however, seems like when I try to attach a dbc file to my simulink model, I’ll get a warning that the dbc file MAY NOT be a valid J1939 database……
again, sometimes I get the message and simulink WILL NOT select the database during config……but other times, I CAN SELECT the same database and get no warnings…??????? j1939 network cnfig bug? MATLAB Answers — New Questions

​

Hello everyone!
I am currently working with Matlab R2021b on Ubuntu 22.04. I am writing a program that will solve linear sparse system of equations, and I am trying to solve it with Pardiso function from Intel MKL library.
I have successfully compiled Mex file in C (that is calling Pardiso function); however when I call Mex file, Matlab crashes without any output or message.

I am using GCC to compile my Mex file. This is the code for compilation:
mex CFLAGS="$CFLAGS -m64" LDFLAGS="$LDFLAGS -m64" pardiso_sparse.c -I/opt/intel/oneapi/mkl/2024.2/include -L/opt/intel/oneapi/mkl/2024.2/lib -lmkl_intel_lp64 -lmkl_gnu_thread -lmkl_core -lgomp -lpthread -lm -ldl

Here is the Mex file itself. It is quite short:
#include "mex.h"
#include "math.h"
#include "matrix.h"
#include "stdint.h"
#include "mkl.h"

/*#define PARDISO_MAX_SIZE 64*/

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{

double *vals, *RG;
int32_t *rows, *cols;
int32_t n;

// input
vals = mxGetPr(prhs[0]);
rows = (int32_t*) mxGetData(prhs[1]);
cols = (int32_t*) mxGetData(prhs[2]);
RG = mxGetPr(prhs[3]);
n = (int32_t) mxGetScalar(prhs[4]);

// Create the output solution vector x (n x 1)
plhs[0] = mxCreateDoubleMatrix(n, 1, mxREAL);
double *sol = mxGetPr(plhs[0]);

// PARDISO control parameters
void *pt[64] = {0}; // Internal solver memory pointer
MKL_INT iparm[64] ; // Control parameters
MKL_INT maxfct = 1, mnum = 1, phase, error = 0, msglvl = 1;
MKL_INT mtype = 2; // Real symmetric matrix
MKL_INT n_mkl = (MKL_INT)n;
MKL_INT nrhs_mkl = 1; // Number of right-hand sides (for solving Ax=b)
MKL_INT idum; // Dummy integer used for PARDISO
double ddum;

// Initialize PARDISO control parameters to default
iparm[0] = 0; // Use default PARDISO parameters

// Step 1: Reordering and symbolic factorization
printf("Analysis is running!");
phase = 11;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, NULL, NULL, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during symbolic factorization: %d", error);
}

printf("First phase done!");

// Step 2: Numerical factorization
phase = 22;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, &ddum, &ddum, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during numerical factorization: %d", error);
}

printf("Second phase done!");

// Step 3: Solve and iterative refinement
phase = 33;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, RG, sol, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during solution: %d", error);
}

printf("Thrid phase done!");

// Step 4: Release internal memory
phase = -1;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, &ddum, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, &ddum, &ddum, &error);
}

Can anyone give me a solution to this? I’ve been trying to run this code for more than a week!
AnteHello everyone!
I am currently working with Matlab R2021b on Ubuntu 22.04. I am writing a program that will solve linear sparse system of equations, and I am trying to solve it with Pardiso function from Intel MKL library.
I have successfully compiled Mex file in C (that is calling Pardiso function); however when I call Mex file, Matlab crashes without any output or message.

I am using GCC to compile my Mex file. This is the code for compilation:
mex CFLAGS="$CFLAGS -m64" LDFLAGS="$LDFLAGS -m64" pardiso_sparse.c -I/opt/intel/oneapi/mkl/2024.2/include -L/opt/intel/oneapi/mkl/2024.2/lib -lmkl_intel_lp64 -lmkl_gnu_thread -lmkl_core -lgomp -lpthread -lm -ldl

Here is the Mex file itself. It is quite short:
#include "mex.h"
#include "math.h"
#include "matrix.h"
#include "stdint.h"
#include "mkl.h"

/*#define PARDISO_MAX_SIZE 64*/

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{

double *vals, *RG;
int32_t *rows, *cols;
int32_t n;

// input
vals = mxGetPr(prhs[0]);
rows = (int32_t*) mxGetData(prhs[1]);
cols = (int32_t*) mxGetData(prhs[2]);
RG = mxGetPr(prhs[3]);
n = (int32_t) mxGetScalar(prhs[4]);

// Create the output solution vector x (n x 1)
plhs[0] = mxCreateDoubleMatrix(n, 1, mxREAL);
double *sol = mxGetPr(plhs[0]);

// PARDISO control parameters
void *pt[64] = {0}; // Internal solver memory pointer
MKL_INT iparm[64] ; // Control parameters
MKL_INT maxfct = 1, mnum = 1, phase, error = 0, msglvl = 1;
MKL_INT mtype = 2; // Real symmetric matrix
MKL_INT n_mkl = (MKL_INT)n;
MKL_INT nrhs_mkl = 1; // Number of right-hand sides (for solving Ax=b)
MKL_INT idum; // Dummy integer used for PARDISO
double ddum;

// Initialize PARDISO control parameters to default
iparm[0] = 0; // Use default PARDISO parameters

// Step 1: Reordering and symbolic factorization
printf("Analysis is running!");
phase = 11;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, NULL, NULL, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during symbolic factorization: %d", error);
}

printf("First phase done!");

// Step 2: Numerical factorization
phase = 22;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, &ddum, &ddum, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during numerical factorization: %d", error);
}

printf("Second phase done!");

// Step 3: Solve and iterative refinement
phase = 33;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, RG, sol, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during solution: %d", error);
}

printf("Thrid phase done!");

// Step 4: Release internal memory
phase = -1;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, &ddum, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, &ddum, &ddum, &error);
}

Can anyone give me a solution to this? I’ve been trying to run this code for more than a week!
Ante Hello everyone!
I am currently working with Matlab R2021b on Ubuntu 22.04. I am writing a program that will solve linear sparse system of equations, and I am trying to solve it with Pardiso function from Intel MKL library.
I have successfully compiled Mex file in C (that is calling Pardiso function); however when I call Mex file, Matlab crashes without any output or message.

I am using GCC to compile my Mex file. This is the code for compilation:
mex CFLAGS="$CFLAGS -m64" LDFLAGS="$LDFLAGS -m64" pardiso_sparse.c -I/opt/intel/oneapi/mkl/2024.2/include -L/opt/intel/oneapi/mkl/2024.2/lib -lmkl_intel_lp64 -lmkl_gnu_thread -lmkl_core -lgomp -lpthread -lm -ldl

Here is the Mex file itself. It is quite short:
#include "mex.h"
#include "math.h"
#include "matrix.h"
#include "stdint.h"
#include "mkl.h"

/*#define PARDISO_MAX_SIZE 64*/

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{

double *vals, *RG;
int32_t *rows, *cols;
int32_t n;

// input
vals = mxGetPr(prhs[0]);
rows = (int32_t*) mxGetData(prhs[1]);
cols = (int32_t*) mxGetData(prhs[2]);
RG = mxGetPr(prhs[3]);
n = (int32_t) mxGetScalar(prhs[4]);

// Create the output solution vector x (n x 1)
plhs[0] = mxCreateDoubleMatrix(n, 1, mxREAL);
double *sol = mxGetPr(plhs[0]);

// PARDISO control parameters
void *pt[64] = {0}; // Internal solver memory pointer
MKL_INT iparm[64] ; // Control parameters
MKL_INT maxfct = 1, mnum = 1, phase, error = 0, msglvl = 1;
MKL_INT mtype = 2; // Real symmetric matrix
MKL_INT n_mkl = (MKL_INT)n;
MKL_INT nrhs_mkl = 1; // Number of right-hand sides (for solving Ax=b)
MKL_INT idum; // Dummy integer used for PARDISO
double ddum;

// Initialize PARDISO control parameters to default
iparm[0] = 0; // Use default PARDISO parameters

// Step 1: Reordering and symbolic factorization
printf("Analysis is running!");
phase = 11;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, NULL, NULL, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during symbolic factorization: %d", error);
}

printf("First phase done!");

// Step 2: Numerical factorization
phase = 22;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, &ddum, &ddum, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during numerical factorization: %d", error);
}

printf("Second phase done!");

// Step 3: Solve and iterative refinement
phase = 33;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, vals, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, RG, sol, &error);
if (error != 0) {
mexErrMsgIdAndTxt("MATLAB:mexfile:pardisoError", "Error during solution: %d", error);
}

printf("Thrid phase done!");

// Step 4: Release internal memory
phase = -1;
pardiso(pt, &maxfct, &mnum, &mtype, &phase, &n_mkl, &ddum, (MKL_INT *)rows, (MKL_INT *)cols, &idum, &nrhs_mkl,
iparm, &msglvl, &ddum, &ddum, &error);
}

Can anyone give me a solution to this? I’ve been trying to run this code for more than a week!
Ante matlab, intel mkl, sparse, mex compiler, ubuntu MATLAB Answers — New Questions

​

Good morning, I am asking for a review of the code that I am leaving as an attachment. The objective of the code should be to take as input the position of 4 points inside the square (integral with respect to the square, fixed in the face). Through an analysis of the variation of the 4 x and y coordinates you want to construct a code that simulates the motion of the square built around these points. What is to appear on the screen will be a black square with 4 fixed white points inside moving in a gray plane. The data with the positions are extracted from a text file that comes from a camera capture. Thank you in advance for your help.Good morning, I am asking for a review of the code that I am leaving as an attachment. The objective of the code should be to take as input the position of 4 points inside the square (integral with respect to the square, fixed in the face). Through an analysis of the variation of the 4 x and y coordinates you want to construct a code that simulates the motion of the square built around these points. What is to appear on the screen will be a black square with 4 fixed white points inside moving in a gray plane. The data with the positions are extracted from a text file that comes from a camera capture. Thank you in advance for your help. Good morning, I am asking for a review of the code that I am leaving as an attachment. The objective of the code should be to take as input the position of 4 points inside the square (integral with respect to the square, fixed in the face). Through an analysis of the variation of the 4 x and y coordinates you want to construct a code that simulates the motion of the square built around these points. What is to appear on the screen will be a black square with 4 fixed white points inside moving in a gray plane. The data with the positions are extracted from a text file that comes from a camera capture. Thank you in advance for your help. motion simulation MATLAB Answers — New Questions

​

I have an existing build script for a console application that compiles with the following command:
compiler.build.standaloneApplication(opts)
where opts is of the type:
compiler.build.StandaloneApplicationOptions
Is there a way to tell the runtime to exclude the JVM inside the opts argument so that it will load faster? I basically want to do the same thing as:
mcc -m -R -nojvm
Unfortunately, I don’t want to re-write this build script because it is a standard script we use.
Thanks in advance,
DaveI have an existing build script for a console application that compiles with the following command:
compiler.build.standaloneApplication(opts)
where opts is of the type:
compiler.build.StandaloneApplicationOptions
Is there a way to tell the runtime to exclude the JVM inside the opts argument so that it will load faster? I basically want to do the same thing as:
mcc -m -R -nojvm
Unfortunately, I don’t want to re-write this build script because it is a standard script we use.
Thanks in advance,
Dave I have an existing build script for a console application that compiles with the following command:
compiler.build.standaloneApplication(opts)
where opts is of the type:
compiler.build.StandaloneApplicationOptions
Is there a way to tell the runtime to exclude the JVM inside the opts argument so that it will load faster? I basically want to do the same thing as:
mcc -m -R -nojvm
Unfortunately, I don’t want to re-write this build script because it is a standard script we use.
Thanks in advance,
Dave matlab compiler, jvm, build options, standalone application MATLAB Answers — New Questions

​

Hello everybody, i hope you and everyone here can help me in plotting temporal representations of x1(t), x21(t), x22(t) and x23(t)
This code :
N = 1e3; % Number of data points
T0 = pi; % Can be any number
t = linspace(0, T0, N); % time space with N data points
f0=1000; f1=2000; f2=2000;
x1 = sin(2*pi*f0*t)+sin(2*pi*f1*t)+sin(2*pi*f2*t);
T1 =2*pi; % Can be any number
t1 = linspace(0, T1-pi/N, N); % time space with N data points
x21 = sin(2*pi*f0*t1);
T2 = 3*pi; % Can be any number
t2 = linspace(T1, T2-pi/N, N); % time space with N data points
x22 = sin(2*pi*f1*t2);
T3 = inf;
t2 = linspace(T2, T3, N);
x23 = sin(2*pi*f2*t2);Hello everybody, i hope you and everyone here can help me in plotting temporal representations of x1(t), x21(t), x22(t) and x23(t)
This code :
N = 1e3; % Number of data points
T0 = pi; % Can be any number
t = linspace(0, T0, N); % time space with N data points
f0=1000; f1=2000; f2=2000;
x1 = sin(2*pi*f0*t)+sin(2*pi*f1*t)+sin(2*pi*f2*t);
T1 =2*pi; % Can be any number
t1 = linspace(0, T1-pi/N, N); % time space with N data points
x21 = sin(2*pi*f0*t1);
T2 = 3*pi; % Can be any number
t2 = linspace(T1, T2-pi/N, N); % time space with N data points
x22 = sin(2*pi*f1*t2);
T3 = inf;
t2 = linspace(T2, T3, N);
x23 = sin(2*pi*f2*t2); Hello everybody, i hope you and everyone here can help me in plotting temporal representations of x1(t), x21(t), x22(t) and x23(t)
This code :
N = 1e3; % Number of data points
T0 = pi; % Can be any number
t = linspace(0, T0, N); % time space with N data points
f0=1000; f1=2000; f2=2000;
x1 = sin(2*pi*f0*t)+sin(2*pi*f1*t)+sin(2*pi*f2*t);
T1 =2*pi; % Can be any number
t1 = linspace(0, T1-pi/N, N); % time space with N data points
x21 = sin(2*pi*f0*t1);
T2 = 3*pi; % Can be any number
t2 = linspace(T1, T2-pi/N, N); % time space with N data points
x22 = sin(2*pi*f1*t2);
T3 = inf;
t2 = linspace(T2, T3, N);
x23 = sin(2*pi*f2*t2); tomporal domain MATLAB Answers — New Questions

​

I am simulating the projectile motion of an ariel vehicle with its wing deployment in Simscape Multibody. The ariel vehicle has an initial velocity, which is provided by the State Target in 6 DoF joint. The 6DoF joint is connected to the CG location of the whole vehicle at its initial condition. During the projectile motion of the vehicle, the wing deployment occurs.
The CG location of the vehicle changes at every instant due to the wing deployment. I am measuring this distance using an Inertia Sensor. The whole process until this step works well.
My actual aim is to do a co-simulation between Simscape and a CFD solver. I created a Standalone FMU from the simscape model and provided to the CFD solver. From CFD solver I calculate the aerodynamic forces and moments based on new CG location at every time step and give as input to the simscape multibody model (FMU).
Current difficulty:
My current difficulty is on placing a joint on a moving center of mass (CG) in Simscape Multibody. In simscape, when we place a 6-DoF joint, the center of rotation is always on the joint connection point. In my case, this point is fixed at the initial CG location. Instead, I want to have a joint at the moving CG location. If I have a frame on the moving CG location, I can apply the aerodynamic forces from CFD (using 6 DoF joint) directly at that location. I tried to follow some related answers from https://www.mathworks.com/matlabcentral/answers/1576920-placing-a-joint-on-a-moving-center-of-mass-in-simscape-multibody?s_tid=answers_rc1-2_p2_MLT and https://www.mathworks.com/matlabcentral/answers/1851513-simscape-multibody-apply-forces-sense-velocity-of-global-center-of-mass to resolve this issue, but I could not succed.
I appreciate any suggestion/comment in this regard how to place a frame (joint) on a moving center of mass (CG).
Thanks.I am simulating the projectile motion of an ariel vehicle with its wing deployment in Simscape Multibody. The ariel vehicle has an initial velocity, which is provided by the State Target in 6 DoF joint. The 6DoF joint is connected to the CG location of the whole vehicle at its initial condition. During the projectile motion of the vehicle, the wing deployment occurs.
The CG location of the vehicle changes at every instant due to the wing deployment. I am measuring this distance using an Inertia Sensor. The whole process until this step works well.
My actual aim is to do a co-simulation between Simscape and a CFD solver. I created a Standalone FMU from the simscape model and provided to the CFD solver. From CFD solver I calculate the aerodynamic forces and moments based on new CG location at every time step and give as input to the simscape multibody model (FMU).
Current difficulty:
My current difficulty is on placing a joint on a moving center of mass (CG) in Simscape Multibody. In simscape, when we place a 6-DoF joint, the center of rotation is always on the joint connection point. In my case, this point is fixed at the initial CG location. Instead, I want to have a joint at the moving CG location. If I have a frame on the moving CG location, I can apply the aerodynamic forces from CFD (using 6 DoF joint) directly at that location. I tried to follow some related answers from https://www.mathworks.com/matlabcentral/answers/1576920-placing-a-joint-on-a-moving-center-of-mass-in-simscape-multibody?s_tid=answers_rc1-2_p2_MLT and https://www.mathworks.com/matlabcentral/answers/1851513-simscape-multibody-apply-forces-sense-velocity-of-global-center-of-mass to resolve this issue, but I could not succed.
I appreciate any suggestion/comment in this regard how to place a frame (joint) on a moving center of mass (CG).
Thanks. I am simulating the projectile motion of an ariel vehicle with its wing deployment in Simscape Multibody. The ariel vehicle has an initial velocity, which is provided by the State Target in 6 DoF joint. The 6DoF joint is connected to the CG location of the whole vehicle at its initial condition. During the projectile motion of the vehicle, the wing deployment occurs.
The CG location of the vehicle changes at every instant due to the wing deployment. I am measuring this distance using an Inertia Sensor. The whole process until this step works well.
My actual aim is to do a co-simulation between Simscape and a CFD solver. I created a Standalone FMU from the simscape model and provided to the CFD solver. From CFD solver I calculate the aerodynamic forces and moments based on new CG location at every time step and give as input to the simscape multibody model (FMU).
Current difficulty:
My current difficulty is on placing a joint on a moving center of mass (CG) in Simscape Multibody. In simscape, when we place a 6-DoF joint, the center of rotation is always on the joint connection point. In my case, this point is fixed at the initial CG location. Instead, I want to have a joint at the moving CG location. If I have a frame on the moving CG location, I can apply the aerodynamic forces from CFD (using 6 DoF joint) directly at that location. I tried to follow some related answers from https://www.mathworks.com/matlabcentral/answers/1576920-placing-a-joint-on-a-moving-center-of-mass-in-simscape-multibody?s_tid=answers_rc1-2_p2_MLT and https://www.mathworks.com/matlabcentral/answers/1851513-simscape-multibody-apply-forces-sense-velocity-of-global-center-of-mass to resolve this issue, but I could not succed.
I appreciate any suggestion/comment in this regard how to place a frame (joint) on a moving center of mass (CG).
Thanks. simscape, 6-dof joint, moving center of mass, moving frame, co-simulation, simscape multibody, matlab, simulink MATLAB Answers — New Questions

​

hi, i want to implement Takagi Sugeno model in matlab using FIS, how can i set this rule:
"if x is large and y is small, then z=x+y+2"
of course i can add "If (x is large) and (y is small) then (Z is ‘…’)", but i cannot add "x+y+2", how can i do it?
thanks so muchhi, i want to implement Takagi Sugeno model in matlab using FIS, how can i set this rule:
"if x is large and y is small, then z=x+y+2"
of course i can add "If (x is large) and (y is small) then (Z is ‘…’)", but i cannot add "x+y+2", how can i do it?
thanks so much hi, i want to implement Takagi Sugeno model in matlab using FIS, how can i set this rule:
"if x is large and y is small, then z=x+y+2"
of course i can add "If (x is large) and (y is small) then (Z is ‘…’)", but i cannot add "x+y+2", how can i do it?
thanks so much implement sugeno in fis MATLAB Answers — New Questions

​

I’m getting an error while solving this coupled differential equation usually the error is showing issues with vertical concatenation. here’s the equation i’m tring to solve with mu0 = exp(-T0) and Boundary conditions as : U(y = -1) = 0 and U(y= 1) = 0 and T0 = 0 at y = -1 and T0 = 1 at y = 1.

here’s my code:
% Main script to solve the velocity and temperature profile
clear;
clc;

% Define constants
G = 1; % Source term (example value)
Na = 1; % Nusselt number (example value)

% Define the domain for y
y_span = [-1 1];

%% Step 1: Solve Velocity Equation to get U and dU/dy
% The velocity equation is:
% d/dy (mu0 * dU/dy) = G
% Define the velocity equation as a system of two first-order ODEs
function dU = velocity_ode(y, U, mu0, G)
dU = [U(2); (1/mu0) * G]; % U(1) = U, U(2) = dU/dy
end

% Initial conditions for velocity at y = -1
U_initial = [0; 0]; % U = 0 and dU/dy = 0 at y = -1 (you can adjust this)

% Solve the velocity equation using ode45
[y_vel, U_sol] = ode45(@(y, U) velocity_ode(y, U, mu0, G), y_span, U_initial);

% Extract dU/dy from the solution
dU_dy = U_sol(:, 2); % This is the derivative of U with respect to y

% Plot the velocity profile and its derivative
figure;
subplot(2,1,1);
plot(y_vel, U_sol(:, 1), ‘b-‘, ‘LineWidth’, 2); % U(y)
xlabel(‘y’);
ylabel(‘U(y)’);
title(‘Velocity Profile’);

subplot(2,1,2);
plot(y_vel, dU_dy, ‘r-‘, ‘LineWidth’, 2); % dU/dy(y)
xlabel(‘y’);
ylabel(‘dU/dy’);
title(‘Velocity Gradient Profile’);

%% Step 2: Solve Temperature Equation using dU/dy from Step 1
% Temperature equation: d^2T0/dy^2 + Na * mu0 * (dU/dy)^2 = 0
% Define the temperature equation as a system of two first-order ODEs
function dT = temperature_ode(y, T, dU_dy, mu0, Na)
dT = zeros(2,1); % Initialize the output vector

% Interpolate dU/dy from the previously computed solution
dUdy_squared = interp1(y_vel, dU_dy.^2, y, ‘linear’, ‘extrap’);

% First equation: dT0/dy = T(2)
dT(1) = T(2);

% Second equation: d^2T0/dy^2 = -Na * mu0 * (dU/dy)^2
dT(2) = -Na * mu0 * dUdy_squared;
end

% Initial conditions for temperature at y = -1
T0_initial = [0; 0]; % T0 = 0 and dT0/dy = 0 at y = -1 (adjust second value if needed)
mu0 = ex(-T0); % Viscosity (example value)
% Solve the temperature equation using ode45
[y_temp, T_sol] = ode45(@(y, T) temperature_ode(y, T, dU_dy, mu0, Na), y_span, T0_initial);

% Plot the temperature profile
figure;
plot(y_temp, T_sol(:, 1), ‘r-‘, ‘LineWidth’, 2); % T0
xlabel(‘y’);
ylabel(‘T_0(y)’);
title(‘Temperature Profile’);
Please help me here, thanks in advanceI’m getting an error while solving this coupled differential equation usually the error is showing issues with vertical concatenation. here’s the equation i’m tring to solve with mu0 = exp(-T0) and Boundary conditions as : U(y = -1) = 0 and U(y= 1) = 0 and T0 = 0 at y = -1 and T0 = 1 at y = 1.

here’s my code:
% Main script to solve the velocity and temperature profile
clear;
clc;

% Define constants
G = 1; % Source term (example value)
Na = 1; % Nusselt number (example value)

% Define the domain for y
y_span = [-1 1];

%% Step 1: Solve Velocity Equation to get U and dU/dy
% The velocity equation is:
% d/dy (mu0 * dU/dy) = G
% Define the velocity equation as a system of two first-order ODEs
function dU = velocity_ode(y, U, mu0, G)
dU = [U(2); (1/mu0) * G]; % U(1) = U, U(2) = dU/dy
end

% Initial conditions for velocity at y = -1
U_initial = [0; 0]; % U = 0 and dU/dy = 0 at y = -1 (you can adjust this)

% Solve the velocity equation using ode45
[y_vel, U_sol] = ode45(@(y, U) velocity_ode(y, U, mu0, G), y_span, U_initial);

% Extract dU/dy from the solution
dU_dy = U_sol(:, 2); % This is the derivative of U with respect to y

% Plot the velocity profile and its derivative
figure;
subplot(2,1,1);
plot(y_vel, U_sol(:, 1), ‘b-‘, ‘LineWidth’, 2); % U(y)
xlabel(‘y’);
ylabel(‘U(y)’);
title(‘Velocity Profile’);

subplot(2,1,2);
plot(y_vel, dU_dy, ‘r-‘, ‘LineWidth’, 2); % dU/dy(y)
xlabel(‘y’);
ylabel(‘dU/dy’);
title(‘Velocity Gradient Profile’);

%% Step 2: Solve Temperature Equation using dU/dy from Step 1
% Temperature equation: d^2T0/dy^2 + Na * mu0 * (dU/dy)^2 = 0
% Define the temperature equation as a system of two first-order ODEs
function dT = temperature_ode(y, T, dU_dy, mu0, Na)
dT = zeros(2,1); % Initialize the output vector

% Interpolate dU/dy from the previously computed solution
dUdy_squared = interp1(y_vel, dU_dy.^2, y, ‘linear’, ‘extrap’);

% First equation: dT0/dy = T(2)
dT(1) = T(2);

% Second equation: d^2T0/dy^2 = -Na * mu0 * (dU/dy)^2
dT(2) = -Na * mu0 * dUdy_squared;
end

% Initial conditions for temperature at y = -1
T0_initial = [0; 0]; % T0 = 0 and dT0/dy = 0 at y = -1 (adjust second value if needed)
mu0 = ex(-T0); % Viscosity (example value)
% Solve the temperature equation using ode45
[y_temp, T_sol] = ode45(@(y, T) temperature_ode(y, T, dU_dy, mu0, Na), y_span, T0_initial);

% Plot the temperature profile
figure;
plot(y_temp, T_sol(:, 1), ‘r-‘, ‘LineWidth’, 2); % T0
xlabel(‘y’);
ylabel(‘T_0(y)’);
title(‘Temperature Profile’);
Please help me here, thanks in advance I’m getting an error while solving this coupled differential equation usually the error is showing issues with vertical concatenation. here’s the equation i’m tring to solve with mu0 = exp(-T0) and Boundary conditions as : U(y = -1) = 0 and U(y= 1) = 0 and T0 = 0 at y = -1 and T0 = 1 at y = 1.

here’s my code:
% Main script to solve the velocity and temperature profile
clear;
clc;

% Define constants
G = 1; % Source term (example value)
Na = 1; % Nusselt number (example value)

% Define the domain for y
y_span = [-1 1];

%% Step 1: Solve Velocity Equation to get U and dU/dy
% The velocity equation is:
% d/dy (mu0 * dU/dy) = G
% Define the velocity equation as a system of two first-order ODEs
function dU = velocity_ode(y, U, mu0, G)
dU = [U(2); (1/mu0) * G]; % U(1) = U, U(2) = dU/dy
end

% Initial conditions for velocity at y = -1
U_initial = [0; 0]; % U = 0 and dU/dy = 0 at y = -1 (you can adjust this)

% Solve the velocity equation using ode45
[y_vel, U_sol] = ode45(@(y, U) velocity_ode(y, U, mu0, G), y_span, U_initial);

% Extract dU/dy from the solution
dU_dy = U_sol(:, 2); % This is the derivative of U with respect to y

% Plot the velocity profile and its derivative
figure;
subplot(2,1,1);
plot(y_vel, U_sol(:, 1), ‘b-‘, ‘LineWidth’, 2); % U(y)
xlabel(‘y’);
ylabel(‘U(y)’);
title(‘Velocity Profile’);

subplot(2,1,2);
plot(y_vel, dU_dy, ‘r-‘, ‘LineWidth’, 2); % dU/dy(y)
xlabel(‘y’);
ylabel(‘dU/dy’);
title(‘Velocity Gradient Profile’);

%% Step 2: Solve Temperature Equation using dU/dy from Step 1
% Temperature equation: d^2T0/dy^2 + Na * mu0 * (dU/dy)^2 = 0
% Define the temperature equation as a system of two first-order ODEs
function dT = temperature_ode(y, T, dU_dy, mu0, Na)
dT = zeros(2,1); % Initialize the output vector

% Interpolate dU/dy from the previously computed solution
dUdy_squared = interp1(y_vel, dU_dy.^2, y, ‘linear’, ‘extrap’);

% First equation: dT0/dy = T(2)
dT(1) = T(2);

% Second equation: d^2T0/dy^2 = -Na * mu0 * (dU/dy)^2
dT(2) = -Na * mu0 * dUdy_squared;
end

% Initial conditions for temperature at y = -1
T0_initial = [0; 0]; % T0 = 0 and dT0/dy = 0 at y = -1 (adjust second value if needed)
mu0 = ex(-T0); % Viscosity (example value)
% Solve the temperature equation using ode45
[y_temp, T_sol] = ode45(@(y, T) temperature_ode(y, T, dU_dy, mu0, Na), y_span, T0_initial);

% Plot the temperature profile
figure;
plot(y_temp, T_sol(:, 1), ‘r-‘, ‘LineWidth’, 2); % T0
xlabel(‘y’);
ylabel(‘T_0(y)’);
title(‘Temperature Profile’);
Please help me here, thanks in advance ode45, matlab, coupledode MATLAB Answers — New Questions

​

I am getting the following error when I type import matlab.engine in python IDLE

Traceback (most recent call last):
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 45, in <module>
pythonengine = importlib.import_module("matlabengineforpython"+_PYTHONVERSION)
File "C:UsersudvasAppDataLocalProgramsPythonPython37libimportlib__init__.py", line 127, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
File "<frozen importlib._bootstrap>", line 1006, in _gcd_import
File "<frozen importlib._bootstrap>", line 983, in _find_and_load
File "<frozen importlib._bootstrap>", line 965, in _find_and_load_unlocked
ModuleNotFoundError: No module named ‘matlabengineforpython3_7’

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 61, in <module>
pythonengine = importlib.import_module("matlabengineforpython"+_PYTHONVERSION)
File "C:UsersudvasAppDataLocalProgramsPythonPython37libimportlib__init__.py", line 127, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
File "<frozen importlib._bootstrap>", line 1006, in _gcd_import
File "<frozen importlib._bootstrap>", line 983, in _find_and_load
File "<frozen importlib._bootstrap>", line 967, in _find_and_load_unlocked
File "<frozen importlib._bootstrap>", line 670, in _load_unlocked
File "<frozen importlib._bootstrap>", line 583, in module_from_spec
File "<frozen importlib._bootstrap_external>", line 1043, in create_module
File "<frozen importlib._bootstrap>", line 219, in _call_with_frames_removed
ImportError: DLL load failed: The specified procedure could not be found.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "<pyshell#1>", line 1, in <module>
import matlab.engine
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 64, in <module>
‘MathWorks Technical Support for assistance: %s’ % e)
OSError: Please reinstall MATLAB Engine for Python or contact MathWorks Technical Support for assistance: DLL load failed: The specified procedure could not be found.

I am using Matlab 2020a and python 3.7.1.I am getting the following error when I type import matlab.engine in python IDLE

Traceback (most recent call last):
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 45, in <module>
pythonengine = importlib.import_module("matlabengineforpython"+_PYTHONVERSION)
File "C:UsersudvasAppDataLocalProgramsPythonPython37libimportlib__init__.py", line 127, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
File "<frozen importlib._bootstrap>", line 1006, in _gcd_import
File "<frozen importlib._bootstrap>", line 983, in _find_and_load
File "<frozen importlib._bootstrap>", line 965, in _find_and_load_unlocked
ModuleNotFoundError: No module named ‘matlabengineforpython3_7’

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 61, in <module>
pythonengine = importlib.import_module("matlabengineforpython"+_PYTHONVERSION)
File "C:UsersudvasAppDataLocalProgramsPythonPython37libimportlib__init__.py", line 127, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
File "<frozen importlib._bootstrap>", line 1006, in _gcd_import
File "<frozen importlib._bootstrap>", line 983, in _find_and_load
File "<frozen importlib._bootstrap>", line 967, in _find_and_load_unlocked
File "<frozen importlib._bootstrap>", line 670, in _load_unlocked
File "<frozen importlib._bootstrap>", line 583, in module_from_spec
File "<frozen importlib._bootstrap_external>", line 1043, in create_module
File "<frozen importlib._bootstrap>", line 219, in _call_with_frames_removed
ImportError: DLL load failed: The specified procedure could not be found.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "<pyshell#1>", line 1, in <module>
import matlab.engine
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 64, in <module>
‘MathWorks Technical Support for assistance: %s’ % e)
OSError: Please reinstall MATLAB Engine for Python or contact MathWorks Technical Support for assistance: DLL load failed: The specified procedure could not be found.

I am using Matlab 2020a and python 3.7.1. I am getting the following error when I type import matlab.engine in python IDLE

Traceback (most recent call last):
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 45, in <module>
pythonengine = importlib.import_module("matlabengineforpython"+_PYTHONVERSION)
File "C:UsersudvasAppDataLocalProgramsPythonPython37libimportlib__init__.py", line 127, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
File "<frozen importlib._bootstrap>", line 1006, in _gcd_import
File "<frozen importlib._bootstrap>", line 983, in _find_and_load
File "<frozen importlib._bootstrap>", line 965, in _find_and_load_unlocked
ModuleNotFoundError: No module named ‘matlabengineforpython3_7’

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 61, in <module>
pythonengine = importlib.import_module("matlabengineforpython"+_PYTHONVERSION)
File "C:UsersudvasAppDataLocalProgramsPythonPython37libimportlib__init__.py", line 127, in import_module
return _bootstrap._gcd_import(name[level:], package, level)
File "<frozen importlib._bootstrap>", line 1006, in _gcd_import
File "<frozen importlib._bootstrap>", line 983, in _find_and_load
File "<frozen importlib._bootstrap>", line 967, in _find_and_load_unlocked
File "<frozen importlib._bootstrap>", line 670, in _load_unlocked
File "<frozen importlib._bootstrap>", line 583, in module_from_spec
File "<frozen importlib._bootstrap_external>", line 1043, in create_module
File "<frozen importlib._bootstrap>", line 219, in _call_with_frames_removed
ImportError: DLL load failed: The specified procedure could not be found.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "<pyshell#1>", line 1, in <module>
import matlab.engine
File "C:UsersudvasAppDataLocalProgramsPythonPython37libsite-packagesmatlabengine__init__.py", line 64, in <module>
‘MathWorks Technical Support for assistance: %s’ % e)
OSError: Please reinstall MATLAB Engine for Python or contact MathWorks Technical Support for assistance: DLL load failed: The specified procedure could not be found.

I am using Matlab 2020a and python 3.7.1. python, engine MATLAB Answers — New Questions

​

Please how do I setup a DSN-less database connection without using the Database Toolbox GUI on MATLAB R2018b?Please how do I setup a DSN-less database connection without using the Database Toolbox GUI on MATLAB R2018b? Please how do I setup a DSN-less database connection without using the Database Toolbox GUI on MATLAB R2018b? dsn-less, sql server MATLAB Answers — New Questions

​