Answered

Heat transfer differential equation

Call the function this way: [t,y] = ode45(@f,[0,5],1673)

Heat transfer differential equation

Call the function this way: [t,y] = ode45(@f,[0,5],1673)

10 hours ago | 0

Answered

2D information to 1D colorplot

Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

2D information to 1D colorplot

Simple example x = 0:10; y = x.^2; temp = rand(1,10); % temperature cm = jet(10); % colormap ind = 1+rou...

11 hours ago | 0

Answered

Quiver with a system of ODEs.

For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

Quiver with a system of ODEs.

For vector field u = x'; v = y'; To plot the curve [T,Y] = ode45 ... plot(Y(:,1),Y(:,2)) % (x,y) I suppose x and y depen...

23 hours ago | 0

Answered

Solving nonLinear equation and Linearized Eqaution using ODE45

YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

Solving nonLinear equation and Linearized Eqaution using ODE45

YOu have to express , and . Solve the system function du = nonlin(t,y) fi = y(1); psi = y(2); theta = y(3); ...

23 hours ago | 0

Answered

How to plot magnitude and directions between 2 points?

try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

How to plot magnitude and directions between 2 points?

try this (not tested) data = load ('Input.txt'); [row, col] = size (data); Y2 = data (:,1); X2 = data (:,2); Y1 = data (:,3...

1 day ago | 0

| accepted

Answered

surface fixed point value Z

I want to know the value of Z from the surfc function Here is it

surface fixed point value Z

I want to know the value of Z from the surfc function Here is it

4 days ago | 0

Answered

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?

Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

Where are the bugs for this ODE finite difference problem that solve using Newton Raphson method?

Here is my attempt clc,clear N = 10; % 10 Set parameters L = 16*0.3048; %meter b = 19.625 * 0.0254; d = 1.625 * 0.0254; ...

8 days ago | 1

| accepted

Answered

I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.

Try @(V,T)

I am receiving the error message Index exceeds the number of array elements (1) when trying to use ode45.

Try @(V,T)

10 days ago | 0

| accepted

Answered

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?

alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

I have a 14 equation system of second order differential equations which I am trying to solve using dsolve. The program just keeps running without finishing ?

alpha and E should be numerical function main m = 15000; k = 150000; n = 14; A = zeros(n); A(n,n) = -k/m; A(n,n-1) = k/m;...

11 days ago | 0

Answered

Euler 3D rotation between two vectors

Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

Euler 3D rotation between two vectors

Here is rotation using Rodriguez matrix Here is rotation about Y axis and Z. Here is the difference between them See ...

11 days ago | 0

| accepted

Answered

Simultaneous differential equations - derivative of y^3 wrt t

Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

Simultaneous differential equations - derivative of y^3 wrt t

Just get dydt(1)=(-3*D*Cs/rho/r0^2)*y(1)*(1-y(2)); dydt(1)=nthroots(dydt(1),3);

11 days ago | 0

Answered

Numerically derive a continous, non-symbolic function

Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

Numerically derive a continous, non-symbolic function

Derivative is (if you have numerical data) dy = (y(i)-y(i-1)) / (t(i)-t(i-1)); Maybe diff (if you have a function) ? syms x ...

13 days ago | 0

Answered

3D Point Cloud - Gaussian Curvature

one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

3D Point Cloud - Gaussian Curvature

one way clc,clear % generate some data r = 3; t = linspace(0,10)'; x = r*cos(t); y = r*sin(t); z = sin(1*t); t1 = ...

14 days ago | 0

Answered

Runge Kutta for system of eqs

I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

Runge Kutta for system of eqs

I wrote function in the way similar to ode45 y_dot = @(x,dx)[dx; -4*dx-5*x]; %insert function to be solved for i = 1:n k1...

15 days ago | 0

Answered

The problem of finding intersection point of triangle mesh

To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

The problem of finding intersection point of triangle mesh

To check every edge: for j = 1:size(face,2) i1 = 1; for i2 = [2 3 1] % checking edges [1 2], [2 3], [3 1...

16 days ago | 0

| accepted

Answered

How can i solve a system which contains matrices

Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

How can i solve a system which contains matrices

Here is the best way i know: A = [5 3 1 1 1 1 3 1 2]; B = [c1(:) c2(:) c3(:)]'; C = inv(A) * B;

18 days ago | 1

| accepted

Answered

Using a for-loop to score values in different intervals

Correct form if 0 <= x(i) && x(i) < 1

Using a for-loop to score values in different intervals

Correct form if 0 <= x(i) && x(i) < 1

19 days ago | 0

| accepted

Answered

Intersection between line and circle when line ends inside circle (using geom2d)

Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

Intersection between line and circle when line ends inside circle (using geom2d)

Someone should copy their comment to the answers section so it can be accepted. It should be me % Position of line: x1 = [0...

19 days ago | 1

| accepted

Answered

Removing points from meshgrid (mask)

Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

Removing points from meshgrid (mask)

Try this simple script [X,Y] = meshgrid(1:20); n = [1 -1]; n = n/norm(n); % normal of a line dist = -7; ...

20 days ago | 0

Answered

How can I find the surface normal of a 8 by 3 matrix?

Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

How can I find the surface normal of a 8 by 3 matrix?

Use griddata to create a surface Use surfnorm to find normal vectors [nx,ny,nz] = surfnorm(x,y,z); surf(x,y,z) hold on quiv...

20 days ago | 0

Answered

How do I increase the circle size and animate a stick figure?

Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

How do I increase the circle size and animate a stick figure?

Use for loop and pause t = linspace(0,2*pi); x = cos(t); y = sin(t); axis([-10 10 -10 10]) hold on for i = 1:10 h = p...

20 days ago | 0

Answered

Finding specific data point x for given y

I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

Finding specific data point x for given y

I need to find the CB and CS values at t=5 Use interp1 to find i need to find the t and CB values when CS = 0.001 Use fsolve ...

20 days ago | 0

Answered

Surface between 3D lines

You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

Surface between 3D lines

You are concatenating in a wrong way XX = [D1 D1 D1]; YY = [F1 F2 F3]; ZZ = [z1 z2 z3]; surf(XX,ZZ,YY,'EdgeColor','none')

20 days ago | 0

| accepted

Answered

two body problem using ode45

- equation of equilibrium Is this condition fullfield?

two body problem using ode45

- equation of equilibrium Is this condition fullfield?

21 days ago | 0

Answered

Numerically Integrating the differential equation below?

Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

Numerically Integrating the differential equation below?

Any guidance on it is much appreciated! Helpfull page: ODE45 I do not know what I can do if I no tf. Try tf = 5. This should ...

22 days ago | 0

Answered

Streamline plot not complete

Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

Streamline plot not complete

Works ok for me N = 15; xx = linspace(min(x(:)),max(x(:)),10); yy = linspace(min(y(:)),max(y(:)),10); [startx,starty] = mes...

22 days ago | 0

| accepted

Answered

2D Matrix revolution

Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

2D Matrix revolution

Use cylinder for each curve separately: z = linspace(-1,1,20); % axis r = sqrt(1-z.^2); ...

22 days ago | 1

Answered

How do I align the vertices and the adjacency graph of a Buckyball?

Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

How do I align the vertices and the adjacency graph of a Buckyball?

Here is a way ind = full(B); % convert to ordinary matrix [rr,cc] = find(ind); % find non-zero elements ix = ...

24 days ago | 0

Answered

Find the coordinates of a point chosen on a plot

Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

Find the coordinates of a point chosen on a plot

Short example x = linspace(0,2,20)'; % generate some data y = sin(x); plot(x,y,'.-b') xy = ginput(1); ...

26 days ago | 1

| accepted

Answered

How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?

Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

How do I create a sphere without sphere function that will inflate like a balloon utilizing the getframe function?

Use surf to create a 3D surface rr = linspace(0,0.5,20); tt = linspace(0,2*pi,40); [R,T] = meshgrid(rr,tt); [X,Y] = pol2cart...

28 days ago | 0